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PREDICTION  OF  TROPOSPHERIC  RADIO  TRANS  MISS  ION  LOSS 
OVER  IRREGULAR  TERRAIN 

A  COMPUTER  METHOD  -  1968 


by 


A.  G.  Longley  and  P.  L.  Rice 


This  report  describes  a  computer  method  for  predicting  ’"'ng- 
term  median  transmission  loss  over  irregular  terrain.  The 
method  is  applicable  for  radio  frequencies  above  20  MHz  and 
may  be  used  either  with  detailed  terrain  profiles  for  actual 
paths  or  with  profiles  that  are  representative  of  median  terrain 
characteristics  for  a  given  area.  Estimates  of  variability  in 
time  and  with  location,  and  a  method  for  computing  service 
probability,  are  included. 

KEYWORDS:  transmission  loss,  tropospheric  propagation,  irregular 
terrain,  time  availability,  service  probability 

1.  INTRODUCTION 

This  report  describes  a  computer  method  for  predicting  long-t  r:.n 
median  radio  transmission  loss  over  irregular  terrain.  The  method  is 
based  on  well-established  propagation  theory  and  has  been  tested  against 
n  large  number  of  propagation  measurements.  It  is  applicable  for  radio 
frequencies  above  20  MHz  and  may  be  used  either  with  detailed  terrain 
profiles  for  actual  paths  or  with  profiles  that  are  representative  of 
median  terrain  characteristics  for  a  given  area.  Estimates  of  median 
terrain  characteristics  are  based  on  a  large  number  of  terrain  profiles 
for  several  types  of  terrain,  includir  ;  plains,  desert,  rolling  hills,  foot¬ 
hills,  and  rugged  mountains. 

Given  radio  frequency,  antenna  heights,  and  an  estimate  of  terrain 
irregularity,  median  reference  values  of  attenuation  relative  to  die  trans¬ 
mission  loss  in  free  space  are  calculated  as  a  function  of  distance.  For 
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radio  line-of-sight  paths,  the  calculated  reference  is  based  on  two-ray 
theory  and  an  extrapolated  value  of  diffraction  attenuation.  For  trana- 
horizon  paths,  the  reference  value  is  either  diffraction  attenuation  or 
forward  scatter  attenuation,  whichever  is  smaller. 

This  prediction  method  was  developed  for  use  with  a  digital  com¬ 
puter  and  has  been  made  sufficiently  general  to  provide  estimates  of 
transmission  loss  expected  over  a  wide  range  of  frequencies,  path 
lengths,  and  antenna  height  combinations,  over  smooth  to  highly  rugged 
terrain,  and  tor  both  vertical  and  horizontal  polarization.  The  method 
is  described  in  complete  detail  in  annex  3.  Shortcuts  appropriate  for 
limited  applications  are  indicated  throughout  the  body  of  the  report. 
Familiarity  with  other  propagation  models,  such  as  those  described  by 
Rice  et  al,  (1967)  is  not  essential  for  using  the  prediction  method  des¬ 
cribed  here. 

Predictions  have  been  tested  against  data  for  wide  ranges  of  fre¬ 
quency,  antenna  height  and  distance,  and  for  all  types  of  terrain  from 
very  sir^th  p’  .ins  to  extremely  rugged  mountains.  The  data  base  in¬ 
cludes  more  than  500  long-term  recordings  throughout  the  world  in  the 
frequency  range  40  to  10,000  MHz,  and  several  thousand  mobile  re¬ 
cordings  in  the  United  States  at  frequencies  from  20  to  1000  MHz  .  The 
method  is  intended  for  use  within  the  following  ranges: 


Parameter 
frequency 
antenna  heights 
distance 

surface  refractivity 


Range 

20  to  40,  000  MHz 
0.  5  to  3,  000  m 
1  to  2 , 000  km 
250  to  400  N-units 


In  applying  this  prediction  method  to  specific  paths  for  which  de¬ 
tailed  profiles  are  available,  certain  limitations  on  antenna  siting  are 
desirable.  For  example,  the  angle  of  elevation  of  each  horizon  ray 
above  the  horizontal  should  not  exceed  12°,  and  the  distance  from  each 
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antenna  to  its  horizon  should  not  be  less  than  1/  10  ,  or  more  than  three 
times,  the  corresponding  smooth-earth  distance. 

Section  2  discusses  atmospheric  and  terrain  parameters,  and  sec¬ 
tion  3  explains  how  transmission  loss  is  calculated.  The  topics  treated 
in  the  annexes  are  as  fellows: 

Annex  1  shows  how  the  computed  reference  values,  A  ,  are  ad- 

cr 

justed  to  provide  long-term  median  estimates,  A(0.5),  of  attenuation 
relative  to  free  space  for  any  given  set  of  data.  This  annex  also  gives 
estimates  of  the  variability  in  time  and  with  location  and  shows  how  to 
estimate  prediction  errors. 

Annex  2  shows  how  various  path  parameters  have  been  derived 
from  studies  of  terrain  profiles. 

Annex  3  gives  detailed  formulas  and  procedures  required  to  calcu¬ 
late  the  median  reference  attenuation  A  .  A  computer  program  listing, 

cr 

flow  diagram,  and  sample  computations  are  included. 

Annexes  1  and  3  each  contain  a  list  of  symbols  used  in  that  annex, 
with  their  definitions. 

2.  ATMOSPHERIC  AND  TERRAIN  PARAMETERS 
2.  1  Atmospheric  Effects 

Radio  transmission  loss  in  tropospheric  propagation  depends  on 
characveristics  of  the  atmosphere  and  of  terrain.  For  predicting  a  long¬ 
term  median  reference  value  of  transmission  lo3S,  the  refractive  index 
gradient  near  the  earth's  surface  is  the  most  important  atmospheric 
parameter.  This  surface  gradient  largely  determines  the  bending  of  a 
radio  ray  as  it  passes  through  the  atmosphere.  Rays  may  be  repre¬ 
sented  as  straight  lines,  within  the  first  Kilometer  above  the  earth's  sur¬ 
face,  if  an  "effective  earth's  radius",  a  ,  is  defined  as  a  function  of  the 

refractivity  gradient  or  of  the  mean  surface  refractivity,  N  .  In 
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calculating  the  long-term  reference  value,  the  minimum  monthly  mean 
value  of  Ng  is  chosen  to  characterize  average  atmospheric  conditions. 

The  effective  earth's  radius,  which  allows  for  regional  differences 
in  average  atmospheric  conditions,  is  defined  as 

a  =  6370  [l  -  0.04665  exp(0.  005577  N  )]~l  km,  (1) 

where  the  actual  radius  of  the  earth  is  taken  to  be  6370  km. 

The  minimum  monthly  mean  value  of  surface  refractivity  may  be  ob¬ 
tained  from  measurements  or  from  maps  showing  a  related  parameter, 

N  .  The  refractivity,  N  ,  represents  surface  refractivity  reduced  to 
o  o 

sea  level.  Figure  1,  reproduced  from  Bean,  Horn,  and  Ozanich  (I960) 
shows  minimum  monthly  mean  values  of  throughout  the  world.  The 
corresponding  surface  refractivity  is  then 

N  =  N  exp(-  0.  1057  h  ),  (2) 

where  h  is  the  elevation  of  the  earth's  surface  in  kilometers  above 

a 

mean  sea  level.  The  elevation  h  is  determined  at  the  base  of  the  lower 

s 

antenna  for  line-of-sight  paths.  For  a  transhorizon  path,  N  is  taken 

s 

as  the  average  of  two  values  computed  by  substituting  the  heights  of  the 
horizon  obstacles,  h.  and  h  ,  in  (2).  If  an  antenna  is  more  than 

Jkj  1  Lu 

150  m  below  its  horizon,  h  and  N  should  be  determined  at  the  antenna 

s  o 

location.  A  commonlv  used  value  of  N  is  301  with  an  effective  earth's 

s 

radius  a  =  8497  km,  which  corresponds  to  4/  3  of  the  actual  radius. 

Other  atmospheric  effects,  such  as  changes  in  the  refractive  index, 
changes  in  the  amount  of  turbulence  or  stratification,  as  well  as  absorp¬ 
tion  by  oxygen,  water  vapor,  clouds,  and  precipitation,  are  allowed  for 
by  empirical  adjustments.  Adjustments  to  the  median,  allowing  for 
differences  in  climate,  and  estimates  of  variability  relative  to  the  median 
are  described  in  annex  1. 
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MINIMUM  MONTHLY  SURFACE  REFRACTIVITY  VALUES  REFERRED  TO  MEAN  SEA  LEVEL 
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2.2  Description  o£  Terrain 

Transmission  loss  may  be  calculated  for  specific  paths  where  de¬ 
tailed  profiles  are  available,  but  the  prediction  method  is  particularly 
useful  when  little  is  known  of  the  details  of  terrain  for  actual  paths.  To 
characterize  terrain,  profiles  may  be  read  at  regular  intervals  in  both 
N-S  and  E-W  directions,  forming  a  uniform  grid  over  an  area.  Or  an 
actual  or  proposed  deployment  of  propagation  paths,  representing  a  wide 
variety  of  terrain  conditions,  may  be  combined  to  provide  a  single  set 
of  profiles  for  which  an  estimate  of  median  propagation  conditions  is  de¬ 
sired.  The  interdecile  range,  Ah(d)  ,  of  terrain  heights  above  and  be.1  . 
a  straight  line  fitted  to  elevations  above  sea  level,  is  calculated  at  fixed 
distances.  Usually  median  values  of  Ah(d)  increase  with  path  length  to 
an  asymptotic  value,  Ah,  which  is  used  to  characterize  the  terrain. 

When  terrain  profiles  are  not  available,  estimates  of  Ah  may  be 
obtained  from  table  1. 


Table  1.  Estimates  of  Ah. 


Type  of  Terrain  Ah  in  Meters 


Water  or  very  smooth  plains 

0 

-  5 

Smooth  plains 

5 

-  20 

Slightly  rolling  plains 

20 

-  40 

Rolling  plains 

40 

-  80 

Hills 

80 

-  150 

Mountains 

150 

-  300 

Rugged  mountains 

300 

-  700 

Extremely  rugged  mountains 

>700 
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Median  estimates  of  Ah(d)  at  desired  distances  may  be  obtained  from 
the  following  relationship,  which  is  based  on  a  study  of  a  large  number 
of  profiles: 


Ah(d)  =  Ah  [)  -  0.  8  exp(-  0.  02  d)]  m,  (3) 

where  Ah(d)  and  Ah  are  in  meters  and  the  distance  d  is  in  kilometers. 
Studies  of  terrain  are  described  in  annex  2. 

2.  3  Parameters  Required  to  Compute  Transmission  Loss 

For  any  specific  application,  a  minimum  of  four  essential  param¬ 
eters  must  be  supplied  in  order  to  calculate  reference  values  of  trans¬ 
mission  loss.  These  are  the  carrier  frequency  f  in  megahertz,  the 
path  distance  d  in  kilometers,  and  the  transmitting  and  receiving  an¬ 
tenna  heights  above  ground  h  and  h  in  meters.  Other  path  param- 

gi  g  i 

eters  used  in  the  computations,  such  as  horizon  distances  and  elevation 
angles,  may  be  derived  from  these  values  and  available  terrain  infor¬ 
mation  as  described  in  the  next  subsection.  When  detailed  profiles  for 
;ndividual  paths  are  available  one  may  compute  these  additional  path 
parameters  using  the  methods  outlined  in  annex  3. 

In  addition  to  estimates  of  surface  refractivity  N  and  the  terrain 

6 

parameter  Ah,  previously  discussed,  the  ground  constants  applicable  to 
the  intervening  terrain  should  be  considered.  The  conductivity,  cr  ,  of 
the  earth's  surface  and  its  permittivity  or  relative  dielectric  constant,  t, 
enter  into  the  calculations  for  line -of-sight  and  diffraction  attenuation. 
When  these  constants  are  not  known  for  a  given  path  or  area  the  following 
values  may  be  assumed: 
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Table  2.  Typical  Ground  Constants 


Type  of  surface 

a  mho  /  m 

i 

Poor  ground 

0.  001 

4 

Average  ground 

0.  005 

15 

Good  ground 

0.  02 

25 

Sea  water 

5 

81 

Fresh  water 

0.  01 

81 

At  sufficiently  low  frequencies,  the  effect  of  the  conductivity  cr  is  dom¬ 
inant,  while  at  sufficiently  high  frequencies,  the  dielectric  constant  < 
has  the  dominant  effect.  For  oversea  transmission,  this  transition  oc¬ 
curs  between  300  and  3000  MHz,  while  over  "average"  ground,  the 
transition  is  between  5  and  50  MHz  .  For  propagation  over  irregular 
terrain,  at  frequencies  above  100  MHz  ,  and  with  antennas  more  than 
5  m  above  ground,  the  effects  of  the  ground  constants  are  slight,  and  the 
results  for  horizontal  and  vertical  polarization  are  nearly  the  same. 

Under  these  conditions  the  method  may  be  simplified  considerably  by  as¬ 
suming  the  magnitude  of  the  theoretical  reflection  coefficient  R  =  0.  95 

h,  v 

and  the  phase  shift  c  =  0  .  For  many  applications  this  results  in  an  esti¬ 
mate  of  the  effective  reflection  coefficient  R  =0.9,  and  the  attenuation 

e 

A  in  (3.  2)  may  be  calculated  directly,  bypassing  the  calculations  shown 
in  equations  (3.  5)  through  (3.  15)  in  annex  3.  These  approximations  are 
not  applicable  for  transmission  over  the  sea.  For  oversea  transmission 
R^  ^  and  c  may  be  computed  using  the  equations  given  in  annex  3  or 
estimated  from  figure  IH.  1  or  figure  HI.  5  of  annex  IH,  volume  II  of  the 
report  by  Rice  et  al.  (1967). 

The  parameters  required  for  computing  reference  values  of  trans¬ 
mission  loss  L  ^  ,  or  the  corresponding  attenuation  below  free  space 
A  ,  are  then:  frequency  f  in  megahertz,  distance  d  in  kilometers, 

antenna  heights  above  ground  h  and  h  in  meters,  as  well  as 

gl  gi 
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estimates  of  surface  refractivity  N  ,  the  terrain  irregularity  Ah  in 

8 

meters,  the  conductivity  cr  mho/  m  ,  and  the  relative  dielectric  constant 
<  of  the  ground. 

The  next  subsection  shows  how  additional  path  parameters  are  ob¬ 
tained  from  these  basic  parameters. 

2.4  Effective  Antenna  Heights,  Horizon  Distances,  and 

Elevation  Angles 

Additional  path  parameters  that  must  be  known  or  estimated  to  cal¬ 
culate  long-term  reference  values  of  transmission  loas  are  the  effective 

antenna  heights  h  and  h  ,  the  horizon  distances  d  ,  and  d  „  ,  and 
ei  ez  LI  L2 

the  horizon  elevation  angles  8  and  8  .  h  rizon  distances  and  ele- 

ei  ez 

vation  angles  are  shown  in  figure  2.  When  a  detailed  terrain  profile  is 

ava;iable  for  a  given  path,  these  parameters  may  be  computed  by  the 

methods  described  by  Rice  et  al.  (1967)  and  outlined  in  annex  3  of  this 

report.  Otherwise,  estimates  of  effective  antenna  heights  and  of  horizon 

distances  and  angles  must  be  calculated.  Such  estimates  are  based  on 

the  terrain  factor  Ah  ,  on  the  antenna  heights  above  ground  h  and  h  , 

gi  g* 

and  on  the  method  used  for  selecting  antenna  sites. 

When  antennas  are  high  and  the  terrain  is  relatively  smooth,  the 
actual  horizon  distances  d  and  d  are  approximately  equal  to  the 

1j  I 

smooth-earth  horizon  distances  d.  and  d,  .  When  antennas  are  low 

Lsl  Lsz 

and  randomly  located  with  respect  to  hills  or  other  obstructions,  as  with 
many  tactical  comrmmication  nets,  the  actual  horizon  distances  will  vary 
greatly  and  their  median  values  may  be  less  than  the  corresponding 
smooth-earth  values.  When  sites  are  chosen  to  take  advantage  of  hill¬ 
tops,  with  propagation  across  valleys,  as  for  radio  relays,  the  horizon 
distances  and  effective  antenna  heights  may  be  greatly  increased.  Con¬ 
sequently  the  following  estimates  are  used: 
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a)  For  net-type  communications,  with  random  antenna  siting,  the 

effective  antenna  heights  h  and  h  are  assumed  to  be  equal  to  the 

6  et  e2 

structural  heights  above  ground: 

h  =  h  m .  (4a) 

el,  2  gi > 2 

b)  For  radio  relay  links,  with  antennas  located  on  or  near  hilltops, 
the  effective  heights  are  larger  than  the  structural  heights  by  an  amount 
whose  median  value  depends  on  the  strucHral  heights  and  the  terrain  ir¬ 
regularity  Ah: 

h  >  =  h„,  ,  +kexp(-2h  /Ah)  m.  (4b) 

Studies  of  terrain  have  shown  that  the  maximum  difference  to  be  ex¬ 
pected  between  median  values  of  struc  ural  and  effective  height  is  50  m, 
in  which  case  k  would  be  equal  to  50  .  But  in  most  situations  such  a 
difference  would  be  unrealistic,  especially  with  low  antennas  and  limited 
freedom  in  site  selection.  Over  moderately  hilly  to  mountainous  terrain, 
with  structural  antenna  heights  less  than  or  equal  to  10  m,  the  following 
estimates  of  k  may  be  used: 

When  antenna  sites  are  rather  carefully  selected  in  an  area  of 
limited  extent,. 

k  =  /I  +  4  sin  (it  h  /  10)  if  0  -  h  -  5  ,  (4c) 

;  g1*2  g»*2 

(  5  otherwise. 

When  antenna  cites  are  still  more  carefully  selected, 

k  =  {1  +  9  sin  (rr  h  /  10)  if  OH  5  5,  (4d) 

I  g»  »  2  gl  ,  2 

v  10  otherwise  . 
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For  each  application  the  question  of  a  suitable  allowance  for  efi^c- 
tive  antenna  height  should  be  carefully  considered.  The  improved  prop¬ 
agation  conditions  that  can  be  obtained  by  careful  site  selection  may  be 
highly  significant.  Further  study  of  definitions  of  effective  antenna 
height  appears  to  be  the  most  urgent  requirement  for  improving  these 
predictions  for  low  antenna  heights  over  irregular  terrain.  Different 

definitions  of  h  may  be  found  appropriate  for  line-of-sight  and  dif- 
61  «  i 

fraction  formulas. 

When  individual  path  profiles  are  not  available,  median  values  of 
the  horizon  distances  d  and  d  are  estimated  as  functions  of  the 

Xu  X  Lj<- 

median  effective  antenna  heights  h  and  h  determined  above,  the 

6  ei  ez 

terrain  irregularity  factor  Ah,  and  the  smooth-earth  horizon  distances 


d.  and  d„  .  The  distance  from  each  antenna  to  its  horizon  over  a 
Lsi  Lsz 

smooth  earth  is  defined  as 


d  =  */ 0.002a  h  km,  (5a) 

Lsi.z  f  el, l  v 


where  the  effective  antenna  heights  h  are  in  meters  and  the  effec- 

ei  ,  z 

tive  earth’s  radius  a  is  in  kilometers,  as  defined  by  (1).  The  sum  of 


the  smooth-earth  horizon  distances  is 


Median  values  of  horizon  distances  over  irregular  terrain  are  estimated 

as 


=  d 


Lsi ,  * 


exp(-  0.  07  yAh/h^)  km, 


(5c) 


where 
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h  =  /h  for  h  >5  m, 

e  i  ei ,  2  ei ,  2 


!n 


5  otherwise. 


The  total  distance,  d^  ,  between  the  antennas  and  their  horizons  is 

dL=dU  +  dL2  km-  (5d) 

For  paths  whose  antennas  are  within  radio  line  of  sight  of  each 
other,  estimates  of  transmission  loss  depend  on  the  particular  effective 
antenna  heights  that  define  the  dominant  reflecting  plane  between  the  an¬ 
tennas.  Even  for  known  line  -  of- sight  paths  an  estimate  of  the  3um  of 
the  horizon  distances,  d^  ^d,  is  required  to  compute  a  reference  value 

of  attenuation  relative  to  free  space  A  ,  or  of  transmission  loss  L  , 

c  r  c  r 

as  described  in  annex  3,  subject  to  the  restriction 


dLS;=dL2d- 


If  this  condition  is  not  met,a  non-line-of-sight  path  is  implied,  the  esti¬ 
mates  of  h  are  too  low,  and  both  should  be  multiplied  by  the 

ei ,  2 

smallest  factor  that  will  satisfy  (5e). 

The  horizon  elevation  angles  0  and  0  ,  show.  .  figure  2,  are 

ei  e2 

the  angles  by  which  the  horizon  rays  are  elevated,  or  depressed,  rela¬ 
tive  to  the  horizontal  at  each  antenna.  When  detailed  profiles  for  indi¬ 
vidual  paths  are  not  available,  median  values  of  0  may  be  esti- 

ei  ,  z 

mated  as 


C.  1)00  5 


0  1,2  d 


Ls  1  ,  2 


LI,  2 


-  1  Ah  -  4  h 


radians 
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The  sum  of  the  elevation  angles  is 

9=0  +0  or  -  dT  /a  radians  ,  (6b) 

e  ei  e2  L 

whichever  is  larger  algebraically.  In  (6)  all  distances  are  in  kilometers 
and  heights  are  in  meters. 

For  transhorizon  paths  the  path  length  d  is  equal  to  or  greater 
than  the  sum  of  the  horizon  distances  d^.  .  The  angular  distance  for  a 
transhorizon  path  is  always  positive  and  is  defined  as 

9=0  +  d/a  radians,  (7) 

e 

where  d  is  the  path  length  and  a  is  the  effective  earth's  radius,  both  in 
kilometers . 

These  additional  path  parameters,  h  ,  dT  .  ,  0  and  9.  are 

ei ,  2  LI,  2  ei  ,  2  * 

used  in  computing  reference  values  of  attenuation  relative  to  free  space 

A  ,  or  transmission  loss  L  .  When  only  the  basic  parameters  are 
c  r  c  r 

supplied,  estimates  of  these  additional  parameters  are  calculated  using 
equations  (4)  thi  ough  (7).  When  detailed  profiles  are  available  for  de¬ 
sired  paths,  these  additional  parameters  are  obtained  as  described  by 
Rice  et  al.  (  1967)  and  outlined  in  annex  3. 


3.  transmission  loss  calculations 


This  section  describes  how  the  various  parameters  discussed  in 

section  2  are  used  to  compute  transmission  loss.  Median  reference 

values  A  of  attenuation  below  free  space  are  computed  first.  The 
cr 

reference  values  L  of  transmission  loss  are  then  the  sum  of  the  free 

cr 

space  basic  transmission  loss,  L^.,  and  the  reference  attenuation  relative 

to  fre^  space,  A  : 

cr 


L  =  L.  +  A  dB  . 
cr  bf  cr 


(8) 
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The  free-space  basic  transmission  loss  is 


Lb£  =  32.45  +  20  log1()  f  +  20  log1Q  d  dB  .  (9) 

where  the  radio  frequency  f  is  in  megahertz,  and  the  distance  d  is  in 
kilometers . 

The  reference  attenuation  A  is  computed  using  methods  based  on 

different  propagation  mechanisms  for  three  distance  ranges.  Well 

within  radio  line  of  sight,  the  formulas  of  two-ray  optics  are  used  to 

compute  attenuation  relative  to  free  space.  Just  beyond  line  of  sight, 

diffraction  is  the  dominant  mechanism.  The  prediction  method  computes 

a  weighted  average,  ,  of  estimates  of  diffraction  attenuation  over  a 

double  knife  edge,  and  over  irregular  terrain.  At  greater  distances, 

well  beyond  the  radio  horizon,  the  dominant  propagation  mechanism  is 

usually  forward  scatter.  The  prediction  method  for  this  distance  range 

is  a  modification  of  the  scatter  computations  described  by  Rice  et  al. 

(  1967).  The  reference  attenuation  A  for  transhorizon  paths  is  either 

cr 

the  diffraction  attenuation  A  ,  or  the  scatter  attenuation  A  ,  whichever 

d  s 

is  smaller.  The  distance  at  which  diffraction  and  scatter  losses  are 
equal  is  defined  as  d 

To  provide  a  continuous  curve  of  the  computed  reference  attenuation 

A  as  a  function  of  distance  the  following  methods  are  used, 
cr 

3.  1  To  Compute  A  Within  Radio  Line  of  Sight 

When  transmitting  and  receiving  antennas  are  within  radio  line  of 
sight,  the  two-ray  optics  formulas  described  in  annex  3  are  used  to  com¬ 
pute  attenuations  A  and  A,  at  specified  distances  d  and  d,  that  are 

o  l  o  1 

well  within  the  horizon.  The  distance  d^  is  chosen  to  approximate  the 
greatest  distance  at  wh-  h  the  attenuation  below  free  space  is  zero.  The 
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distance  is  greater  than  d  but  well  within  thrt  range  for  which  two- 
ray  optics  formulas  are  valid.  The  methods  described  in  the  next  sub¬ 
section  are  used  to  compute  the  diffraction  attenuation  A^g  at  the 

smooth  earth  horizon  distance  dT  .  These  three  values  of  attenuation, 

Ls 

A  ,  A,  ,  and  A,  ,  computed  at  the  distances  d  ,  d,  ,  and  d_  ,  res- 
o  1  Ls  o  1  Ls 

pectively,  are  used  to  determine  the  slopes  k  and  k  of  a  smooth 

J.  U 

curve  of  A  versus  distance  for  the  range  1  -  d  ^  d,  : 

cr  6  Ls 

Acr=Ac+k1(d-do>+k2l°810<d/do>  dB'  (I°> 

Note  that  the  smooth-earth  horizon  distance  d  may  be  greater  than 

.Ls 

the  actual  horizon  distance  d  . 

Equation  (10)  may  be  simplified  by  defining  a  term  A  : 


A  =  A  -  k  d  -  k  log  d  dB  . 
e  o  1  o  2  °10  o 


Then  for  all  distances  greater  than  one  and  less  than  d  ,  reference  at- 

La  S 

tenuation  is  defined  as 


A_  =  A  +  k  d  +  k  log  d  dB  .  (12) 

c  r  c  i  w  i  v 

Detailed  methods  for  computing  k^  and  k7  are  given  in  annex  3. 

Whenever  the  value  of  A  computed  from  (12)  is  less  than  zero,  it  is 

cr 

assumed  that  A  =  0  .  This  prediction  method  does  not  attempt  to  des- 
c  r 

cribe  the  lobery  within  line  of  sight  over  irregular  terrain. 

3.2  To  Compute  Diffraction  Attenuation  A, 

d 

The  diffraction  attenuation  is  computed  by  combining  estimates  of 
knife-edge  diffraction,  based  on  Fresnel-Kirchhoff  theory,  with  a  modi¬ 
fication  of  the  method  for  computing  diffraction  over  smooth  terrain 
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developed  by  Vogler  (1964).  Vogier's  method  estimates  the  diffraction 
attenuation,  A^_  ,  over  the  bulge  of  the  earth  in  the  far  diffraction  region, 
and  is  applicable  for  smooth  terrain.  Knife-edge  diffraction  theory  is 
used  to  estimate  attenuation  over  an  isolated  hill  or  ridge.  In  this  ap¬ 
plication  the  knife-edge  attenuation,  A  ,  is  computed  as  though  the  radio 
path  crossed  two  sharp,  isolated  ridges.  In  general,  for  irregular  ter¬ 
rain,  the  diffraction  attenuation,  A^  ,  is  computed  as  a  weighted  average 

of  the  two  estimates  A  and  A,  : 

r  k 

A  ,  =  (1  -  w)  A.  +  w  A  dB  ,  ( 13) 

d  k  r 

where  the  weighting  factor,  w,  is  determined  empirically  as  a  function 

of  radio  frequency  and  terrain  parameters  and  is  defined  in  annex  3. 

The  diffraction  attenuation  A,  is  calculated  at  distances  d  and 

d  3 

d^  in  the  far  diffraction  region  using  the  formulas  given  in  annex  3.  A 

straight  line  through  these  points  (A^,  d^)  and  (A^,  d^)  is  then  defined 

bv  the  intercept,  A  ,  ,  and  slope,  m,  ,  as  follows: 

ed  d 

m,  =  (A  -  AJ/(d  -  dj  dB/km,  (14) 

d  4  3  4  3 

and 

A  =  A  +  A  m  d  dB  ,  (15) 

ed  fo  4  d  4 

where  A  is  a  clutter  factor  defined  in  annex  3. 

-O 
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The  reference  attenuation  A  ^  at  any  distance  d  greater  than 

the  smooth  earth  horizon  distance  dT  and  Less  than  the  distance  d 

Ls  x 

where  diffraction  and  scatter  attenuation  are  equal  is 

A  =  A  =  A  +  m,  d  dB,  for  dT  —  d  "S  d  .  (16) 

cr  d  ed  d  Ls  x 

3.  3  To  Compute  Forward  Scatter  Attenuation  A 

s 

When  the  path  Length  d  or  the  angular  distance  9  is  large,  the 

forward  scatter  attenuation  A  may  be  less  than  the  diffraction  atten- 

s 

uation  A  Therefore,  when  the  product  of  the  distance  in  kilometers 
d 

and  the  angular  distance  in  radians  exceeds  0.  5  (0d  >0.5),  forward 

scatter  attenuation  A  is  computed  for  comparison  with  A  . 

s  d 

For  Large  values  of  0d,  the  scatter  attenuation  A  is  assumed  to 
have  a  linear  dependence  on  distance;  therefore,  A^  is  computed  at  two 
large  distances  d  and  d  .  A  straight  line  through  the  points  (A.,  d  ) 
and  (A^.  d^)  is  then  defined  by  the  intercept  A  and  the  slope  m  as 
follows: 


A  =  A  -  m  d  dB  >  (17  a) 

es  5  s  5 

and 

m  -  (A  -  A  )  / (d  -  d_)  dB-  (17b) 

s  6  5  6  o 

The  reference  attenuation  A  at  any  distance  d  greater  than  the 

cr 

distance  d  where  diffraction  and  scatter  attenuations  are  equal  is 
x 

given  by 


A  -A  =  A  +md  dB  ,  for  did.  (  18) 

c  r  s  e  s  s  x 
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The  detailed  equations  nd  techniques  requiied  to  compute  the  ref¬ 
erence  attenuation  A  as  a  function  of  distance  are  given  ir.  annex  3. 

cr 

Figures  3  through  8  show  computed  and  measured  values  -'f  attenua¬ 
tion  below  free  space  plotted  as  a  function  of  distance.  The  data  were 
obtained  in  a  measurement  program  described  by  Johnson  et  al.  (1967). 
Measurements  were  made  at  reque  cies  of  20,  50,  and  100  MHz  ir. 
three  areas.  Figures  3  and  4  show  data  from  northeastern  Ohio,  with 
corresponding  predicted  values.  Figu:  e  3  shows  medians  of  data  re¬ 
corded  at  distances  of  10,  20,  30,  a->d  50  km  from  the  transmitter, 
using  both  vertical  and  horizontal  polarization  at  100  MHz  ,  and  with  re¬ 
ceiver  heights  of  "  6,  and  9  m  .  The  number  of  measurements  at 

each  distance  is  tabulated,  and  corresponding  medians  of  point-to-point 
predictions  based  on  individual  profiles  for  each  of  the  measurement 
f.aths  r  -e  shown.  The  lines  represent  computed  reference  values  A  , 
based  on  frequency,  antenna  heights,  and  the  asymptotic  value  of  terrain 
irregularity  Ah  .  Figure  4  shows  data  and  predicted  values  at  frequen¬ 
cies  of  2  0  and  50  MHz  . 

Figures  5  and  6  show  measured  and  predicted  values  of  attenuation 
below  free  space  for  an  area  in  the  plains  near  Boulder,  Colorado. 
Measurements  were  made  at  distances  of  5  to  80  km  from  the  trans¬ 
mitter.  Medians  of  data  and  the  number  of  measurements  made  at  each 
distance  are  shown,  with  corresponding  point-to-point  predictions.  The 
lines  repre v -it  computed  reference  values  A  ,  based  on  the  asymptotic 
value  of  te  rain  irregularity  Ah. 

Figur-  s  7  and  8  show  measured  and  predicted  values  of  attenuation 
below  free  space  in  the  Colorado  mountains  west  of  Boulder.  The  com¬ 
puted  value  '  of  A  represented  by  lines  assume  Ah  =  650  m.  Actually 
cr 

the  shorter  paths  extend  only  into  the  foothills.  For  the  5-km  paths  the 
median  value  of  Ah(d)  is  100  m.  The  triangles  on  figures  7  and  8  show 
values  computed  using  these  estimates  of  Ah. 
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COMPUTED  AND  MEASURED  VALUES  OF  ATTENUATION  BELOW  FREE  SPACE 

NORTHEASTERN,  OHIO 
f-IOOMHz,  Ah=90m,  H«  =  3I2 


Distance  in  Kilometers 


ED  ANfO  MEASURED  VALUES  OF  ATTENUATION  BELOW  FREE  SPACE 

COLORADO  PLAINS 

f=IOO  MHz,  Ah* 90m,  N,*290 


p 


-23- 


Distance  In  Kilometers 


Figure 
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ANNEX  1 

VARIABILITY  IN  TIME  AND  WITH  LOCATION,  PREDICTION 
ERROR,  AND  SERVICE  PROBABILITY 

This  annex  deals  with  definitions  of  satisfactory  service  for  any 
radio  communication  system,  the  variability  of  service  with  time  and 
location, and  ways  of  allowing  for  random  errors  of  prediction.  A  speci¬ 
fied  grade  of  service,  g  ,  from  a  wanted  signal  in  the  presence  of  un¬ 
wanted  signals,  and  rapid  or  fine-grained  variations  of  these  signal 
levels  with  time  or  receiving  location,  is  assumed  to  depend  only  upon 
whether  a  required  wanted-to-unwanted  signal  ratio  R^g^)  i®  exceeded. 

If  increasing  values  of  available  grades  rep  '  sent  better  gradeB  of  ser¬ 
vice,  R  (g)  must  increase  with  g.  For  any  reasonable  selection  of  a 
receiving  system,  the  probability  of  power-independent  distortion  is  as¬ 
sumed  negligible,  or  at  least  minimized. 

Let  R(q^)  represent  the  available  wanted-to-unwanted  signal  ex¬ 
ceeded  for  a  fraction  q„  of  a  specified  oeriod  of  time.  With  R  and  qm 

T  r  T 

fixed,  satisfactory  service  exists  if  R(q^)  >R  .  Data  normalized  to 
correspond  with  the  given  conditions  can  be  used  to  predict  R{q^,)  •  The 
service  probability,  Q,  is  defined  as  the  expected  fraction  of  normalized 
data  for  which  R(q^,)  >  R  - 

To  avoid  accepting  an  unsatisfactory  system,  a  high  service  pro!  n- 
bility  is  required,  while  there  should  be  a  low  service  probability  for 
any  system  that  is  rejected,  other  things  being  equal.  Service  proba¬ 
bility  can  also  be  used  as  a  weighting  factor  in  cost-benefit  studies  or 
in  studying  complex  netting  and  multiple  access  problems.  For  some 
applications,  such  as  broadcasting  or  relay  networks  with  many  receivers 
and/  or  transmitters,  it  may  be  worthwhile  to  describe  variations  from 
location  to  location  statistically  instead  of  examining  every  possible 

propagation  link.  Then  the  quantity  R(q  ,  q  ,  ''})  is  defined  as  the 

X  Xi 
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wanted-to-unwanted  signal  ratio  exceeded  for  at  least  a  fraction  q  of  a 
specified  period  of  time,  for  a  fraction  q  of  all  receiving  antenna  loca- 
tions  or  propagation  paths  in  a  statistically  homogeneous  area  or  group, 
and  with  a  probability  Q  . 

Interference  probability  is  the  complement,  P  *  1  -  Q,  of  the  ser¬ 
vice  probability.  Harmful  interference  may  be  said  to  exist  at  the  an¬ 
tenna  terminals  of  a  particular  receiving  system  if  the  service  proba¬ 
bility  is  less  than  0.95,  i.  e.,  the  interference  probability  is  greater 
than  0. 05  . 

The  first  section  of  this  annex  discusses  the  time  variability  of  propa¬ 
gation  attenuation,  the  next  section  describes  variability  with  location, 
and  the  last  section  deals  with  the  statistics  of  R  . 

1-1.  Variability  in  Time 

1-1.1  Adjustment  Function  V{0.  5,  d  ) 

c 

Adjustments  to  the  calculated  long-term  reference  value  A  of  at¬ 
tenuation  relative  to  free  space  may  be  required  to  provide  long-term 
median  estimates  of  A(0.  5)  for  given  sets  of  data.  Adjustments  des¬ 
cribed  here  provide  estimates  of  A{0.  5)  for  specific  periods  of  time  and 
for  various  climatic  regions.  Climatic  regions  are  distinguished  on  the 
basis  of  meteorological  data  as  described  by  the  CCIR  (1966)  and  by 
Rice  et  al.  (1967). 

The  calculated  long-term  reference  value  A  represents  the  me¬ 
dian  attenuation  to  be  expected,  assuming  minimum  monthly  mean  values 
of  surface  refractivity  (see  sec.  2.1).  Therefore,  differences 
are  to  be  expected  between  the  reference  value  A  and  the  long-term 
median  A(0.  5)  that  represents  all  hours  of  the  year.  In  the  northern 
temperate  zone,  for  example,  minimum  monthly  mean  values  of  Ng 
occur  during  the  winter  months. 
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The  all-year  median  attenuation  A(0.  5)  differs  from  the  computed 
value  A  by  an  amount  V{0.  5)  dB  : 

A(0.  5)  =  Acr  -  V(0.  5)  dB.  (1.1) 

This  difference  V(0.5)  between  the  all-year  median  attenuation  and  A 

cr 

is  shown  in  figure  1.  1  for  several  climatic  regions  as  a  function  of  an 

effective  distance,  d  ,  expressed  in  kilometers. 

e 

The  effective  distance  d^  depends  on  the  distance  at  which  diffrac¬ 
tion  and  forward  scatter  losses  are  approximately  equal  over  a  smooth 

earth,  and  on  d.  ,  which  is  the  sum  of  the  smooth-earth  horizon  dis- 
Lo 

tances  for  an  effective  earth's  radius  a  =  9000  km.  Define  8  as  the 

si 

angular  distance  where  diffraction  and  scatter  losses  are  approximately 
equal  over  a  smooth  earth  of  effective  radius  a  =  9000  km.  Then, 

d  =a  0  =  65(100/  f)^  km,  and  (1.2a) 

81  81 


where  f  is  the  frequency  in  megahertz  and  h  are  the  effective  an- 

C 1  »  • 

tenna  heights  in  meters,  defined  in  section  2.4  (4).  The  effective  dis¬ 
tance  d  is  then  defined  as  follows: 
e 

For  d  5  (d,  +  d  )  km ,  d=130d/(dT  +d  )  km  ( 1.  3a) 

Lo  si  e  Lo  si 

For  d  >  (d  +  a  )  km ,  d  =  130  +  d-(d,  +d  )  km  .(1.  3b) 

Lo  si  e  '  Lo  si  ' 

In  each  climatic  region  the  calculated  reference  value  should  be  ad¬ 
justed  by  the  amount  V(0.  5,  d^)  to  obtain  the  predicted  all-year  median 
value  of  attenuation. 

The  curves  shown  in  figure  1.  1  represent  V(0.  5,  d  )  for  all  hours 

© 

of  the  day  throughout  the  entire  year.  For  some  applications  it  is  im¬ 
portant  to  know  something  about  the  diurnal  ar.d  seasonal  changes  that 
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ADJUSTMENT  FACTOR  V(0.5,  d.)  FOR  8  CLIMATIC  REGIONS 
A(05)  =  Acr  -  V10.5,  d«)  dB 
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may  be  expected.  Such  changes  have  been  studied  in  the  continental 
United  States  where  a  large  amount  of  data,  recorded  for  periods  of  at 
least  a  year,  is  available.  The  variation  of  long-term  median  levels 
with  season  usually  shows  maximum  attenuation  in  midwinter,  especially 
on  winter  afternoons,  and  minimum  attenuation  during  the  summer 
months,  particularly  during  the  morning  hours. 

The  data  were  divided  into  summer  and  winter  periods,  May  through 
October,  and  November  through  April,  and  the  hours  of  the  day  were  di¬ 
vided  int .  four  groups  providing  the  following  eight  time  blocks: 


No. 

Months 

Hour  b 

1 

Nov  -  Apr 

0600  -  1300 

2 

Nov  -  Apr 

1300  -  1800 

3 

Nov  -  Apr 

1800  -  2400 

4 

May  -  Oct 

0600  -  1300 

5 

May  -  Oct 

1300  -  1800 

6 

May  -  Oct 

1800  -  2400 

7 

May  -  Oct 

0000  -  0600 

8 

Nov  -  Apr 

0000  -  0600 

The  data  for  time  blocks  1,  2,  3,  and  8  were  combined  to  form  a 
winter  all-hours  group  and  the  data  for  time  blocks  4,  5,  6,  and  7  were 
combined  to  form  a  summer  all-hours  group.  The  adjustment  factor 
V(C.5)  for  each  of  these  periods  of  time  is  shown  in  figure  1.2  for  the 
United  States,  which  is  typical  of  a  continental  temperate  climate. 

For  other  climatic  regions  an  indication  of  the  seasonal  variation  to 
be  expected  may  be  obtained  from  the  annual  range  of  monthly  mean  N 
shown  in  figure  1.  3.  Much  of  the  data  in  the  United  States  was  recorded 
where  the  annual  range  of  monthly  mean  N^  ia  40  to  50  N-units  .  In  re¬ 
gions  where  the  annual  range  is  less  than  20  N-units,  seasonal  variations 
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are  expected  to  be  negligible.  One  would  also  expect  less  diurnal  change, 
for  example,  in  a  maritime  temperate  climate  where  changes  in  temper¬ 
ature  during  the  day  are  less  extreme.  In  climates  where  the  surface 
refractivity  N  changes  considerably  throughout  the  year,  the  consecu¬ 
tive  4- to  o-month  period  when  the  monthly  mean  value  of  is  lowest 
is  assumed  to  correspond  to  "winter",  whatever  months  are  involved. 

1-1.  Z  Long-Term  Variability  About  the  Median 

Long-term  variability,  identified  with  the  variability  of  hourly  me¬ 
dian  values  of  attenuation,  usually  results  from  slow  changes  in  average 
atmospheric  refraction,  in  the  defc~ee  of  atmospheric  stratification  or 
in  the  intensity  of  atmospheric  turbulence.  Estimates  of  long-ter;  .*  var¬ 
iability  to  be  expected  are  important  to  insure  adequate  service  and  to 
avoid  possible  interference  between  services  operating  on  the  same  or 
adjacent  frequencies. 

Estimates  of  variability  in  time  about  the  long-term  median  attenua¬ 
tion  are  ba.-ed  on  measurements.  Hourly  median  values  of  transmission 
loss  recorded  for  a  long  period  of  time  over  a  single  path  may  show  wide 
variations,  especially  in  areas  where  marked  seasonal  changes  occur  in 
the  surface  refractivity  and  in  the  refractive  index  gradient. 

Figure  1.4  shows  the  variability  expected  at  100  MHz  in  a  contin¬ 
ental  temperate  climate  for  all  hours  of  the  year,  for  the  summer  hours, 
and  for  the  winter  hours.  The  curves  show  the  interdecile  range  of  at¬ 
tenuation  as  a  function  of  the  effective  distance  d  .  They  were  drawn 

e 

through  medians  of  a  large  amount  of  data  recorded  in  the  frequency 

range  88  to  108  MHz  .  The  curve  Y  (0.  1)  is  the  difference  between  the 

o 

median  attenuation  A(0.  5)  and  the  attenuation  A(0.  1)  not  exceeied  10 
percent  of  the  time.  Similarly,  the  curve  Y^(0.9)  represents  the  dif¬ 
ference  between  the  median  attenuation  and  the  attenuation  A(0.9)  not 
exceeded  90  percent  of  the  time: 
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The  variability  in  time  of  attenuation  over  a  given  path  also  depends 
on  frequency-related  effects.  The  frequency  factors  g(0.  1,  f )  and 
g(0.9,  f)  shown  in  figure  1.5  adjust  the  predicted  variability  for  100  MHz, 
shown  in  figure  1.4  for  use  at  other  frequencies: 

Y(0.  1)  =  Y  (0.  1)  g(0.  1,  f)  dB  (1.5a) 

Y(0.9)  =  Y  (0.9)  g(0.9,  f)  dB  .  (1.5b) 

The  empirical  curves  g(q,  f )  should  not  be  regarded  as  an  estimate  of 
the  dependence  of  long-term  variability  on  frequency;  they  represent 
only  an  average  of  many  effects,  some  of  which  are  frequency-sensitive. 

The  apparent  frequency  dependence  is  a  function  of  the  relative  domin¬ 
ance  of  various  propagation  mechanisms,  which  in  turn  depends  on  cli¬ 
mate,  time  of  day,  season  and  the  particular  types  of  terrain  profiles 
for  which  data  are  available.  An  allowance  for  year-to-year  variability 
is  included  in  g(q,  f). 

The  estimates  of  A  ,  V(Q.  5),  Y(0.  1)  and  Y(0.  9)  may  be  used  to 

cr 

predict  an  entire  cumulative  distribution  of  attenuation  relative  to  free 
space  using  the  following  ratios: 

Y(0.0001)  =  3.33  Y(0.1)  Y(C.9999)  =  2.90  Y(0.  9) 

Y(0.001)  =2.  ?3Y(0.1)  Y (0 .  999)  =  2.  41  Y(0.  9) 

Y(0.01)  =  2.00  Y(0, 1)  Y (0. 99)  =  1.82  Y(0.9) 


YJO.  1)  =  A(0.  5)  -  A(0.  1)  H  1^(0.  5)  -  1^(0.  1)  dB  (1.4a; 

Yo(0.9)  =  A(0.5)  -  A(0.9)  =  1^(0.  5)  -  1^(0.  9)  dB  (1.4b) 
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Figure  1.6  shows  such  a  predicted  distribution  compared  with  values 
measured  over  a  path  for  more  than  2  years.  The  figure  shows  the  bi¬ 
variate  distribution  that  is  typical  of  a  continental  temperate  climate. 

The  preceding  discussion  and  figures  1.2,  1. 4, and  1.  5  are  for  a  con¬ 
tinental  temperate  climate  and  are  based  largely  on  data  from  measure¬ 
ments  made  in  the  United  States  and  in  West  Germany.  Data  from  o.her 
geographical  regions  were  used  to  develop  curves  for  other  climates, 
notably  data  from  the  British  Isles  and  the  west  coastal  areas  of  Europe 
that  characterize  the  maritime  temperate  overland  and  maritime  tem¬ 
perate  oversea  curves.  The  curves  presented  by  Rice  et  al.(1967)  for 
other  climatic  regions  are  based  on  relatively  few  measurements. 


1-1.  3  Computer  Method  for  Estimating  Long-Term  Variability 


Long-term  variability  about  the  median  attenuation  may  be  estimated 

in  terms  of  a  standard  deviation,  cr  ,  ind  a  standard  normal  deviate, 

Ta 

ZQ(q)  .  The  symbol  q,  representing  any  fraction  between  0  and  1,  and 

the  standard  normal  deviate  z  (q)  may  be  expressed  in  terms  of  the  er- 

°  .-1 

ror  function,  erf  x,  and  its  inverse,  erf  x. 


q  =  0.5+0.5erf(z  / \/""2)  (1.6a) 

o  * 

zQ(q)  =  er*_1  (^q-i)  •  (l.fab) 


Then  the  long-term  variability  about  the  median  exceeded  at  least  a  frac¬ 
tion  q  of  the  time  may  be  expressed  as 


-V  *o(qT>  dB' 


(1.  7) 


where  q^  is  any  desired  fraction  of  time  and  Y(q  )  is  the  difference 
between  the  median  attenuation  A(0.  5)  and  the  attenuation  not  exceeded 


a  fraction  q  of  the  time, 
i 
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COMPARISON  OF  MEASURED  AND  PREDICTED  ATTENUATION, 
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For  a  continental  temperate  climate.Y  (0.  1)  is  not  equal  to  Y  (0.9). 

The  cumulative  distribution  may  be  considered  as  consisting  of  two  parts, 

represented  by  standard  deviations  cr  (0.  1)  and  cr  (0.9),  where 

I  d  X  <1 

Y(0.  1)  =  1.282  <r_  (0.1),  and  (1.8a) 

I  a 

Y(0.9)*-1.282<r  (0.9).  (1.8b) 

Define  x  =  d  /  100 . 
e 

For  d  s  200 : 
e 

aTa(0.  1)  =  8  x2g(0.  1,  f)  exp  (  -0 .  36  x2).  (1. 9a) 

For  d  >  200; 

e 

«rT  (0.  1)  =  g(0.  1,  f)  [4.2  +  16.  5  exp ( -0. 77  x)  ]  .  (1.9b) 

For  d  s  250: 

e 

trTa(0. 9)  =  4.  6  x2  g  (0 . 9,  0  exp  (-0. 26  x^) .  (1.10a) 

Ford  >  2  50 : 
e 

<rTa (°- 9)  =  g(0.9,  f)  [2.3  +  15  exp(-0.6x)  ]  .  (1.  10b) 

For  60  s  f  s  1 600  MHz: 

g(0.  1.  f)  =  0. 21  sin  1*5.22  log10(f/ 200)1  +1.  28,  (1.1  1a) 
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g(0.9,  f)  =  o.  18  sin  [5.22  log10  <f/200)J  +1.23.  (1.11b) 

For  f  >  1600: 

g(0.  1,  f)  =  g(0.9,  £)  =  1.05.  (1.  11c) 

Equations  (l.lla-c)  represent  approximations  to  the  empirical  estimates 
of  g(0.  1)  and  g(0.9)  plotted  in  figure  1.  5  for  all  year. 

1-2.  Variability  with  Location 

There  is  a  path-to-path  variation  in  the  available  wanted  Bignal 
power  that  may  be  referred  to  as  location  variability.  Such  random 
variations  from  location  to  location  for  any  given  may  be  assumed 
to  be  normally  distributed  with  a  standard  deviation  cr  dB  .  Then 
A(q^)  is  the  attenuation  below  free  space  not  exceeded  for  at  least  a 
fraction  q^  of  all  randomly  chosen  paths  for  which  other  parameters, 
such  as  frequency,  antenna  heights,  and  path  length,  are  fixed.  A  value 
q^  =  0.5  would  represent  median  conditions,  while  q^  =  0.  1  would 
mean  that  the  antennas  are  assumed  to  be  located  at  sites  selected  to  be 
among  the  best  10  percent  of  all  possible  locations  within  a  given  area. 
For  low  antennas,  over  irregular  terrain,  the  best  estimate  from  pre¬ 
sently  available  data  indicates  a  value  of  a  =  10  dB,  as  shown  in 

J><  C 

annex  2.  The  location  variability  Y(q^)  for  any  fraction  q^  may  be  ex¬ 
pressed  in  terms  of  the  standard  deviation  a,  and  a  standard  normal 

La 

deviate  z  (qT  )  , 
o  J- 

Yt  *  Y(q.  )  =  -  c r  s  (q.  )  dB,  (1.  12a) 

-1— <  L.  A  O  Xj 

where 

Zo(qL)  =  V^erf"1  (2qL  -  D  .  (1.12b) 
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1-3.  Prediction  Error  and  Service  Probability 

A  communication  link  is  assumed  to  provide  satisfactory  service  of 
a  given  grade  g  if  the  available  signal-to-noise  ratio  R  exceeds  the 
required  protection  ratio  R  (g  )  for  at  least  a  fraction  q^.  of  the  time 
and  a  fraction  q^  of  all  similar  links.  The  noise-limited  service  re¬ 
quirement  may  be  written  as 

R(q„,  qT,  Q)  =  R  +  Y_+  YT  +  Y  dB  ,  (1.13a) 

i  J— *  O  1  i— »  c 

where  R  is  the  median  available  signal-to-noise  ratio  and  the  symbols 
o 

Y  ,  Y  ,  and  Y  represent  biases  or  safety  factors  required  to  achieve 
1  -1— <  c 

protection  for  at  least  the  desired  fraction  q^,  of  the  time,  the  desired 
fraction  q^  of  locations,  and  with  a  probability  Q.  These  biases  are 
considered  in  the  order  given,  allowing  first  for  time  availability  q  , 
then  for  location  variability  q^,  with  q^.  fixed,  and  finally  for  service 
probability  Q,  with  both  q^.  and  q^  fixed.  That  is,  we  may  write 
( 1.  13a)  as 

R(qT,  qL.  Q)  =  Rq  +  YT(qT)  +  YJq^  qj  +  Y^qj.,  q^  Q)  dB  .  (1.  13b) 

In  this  discussion  we  consider  interference  due  to  noise  at  the  term¬ 
inals  of  the  receiving  antenna,  including  both  internal  system  noise  and 
externally  generated  unwanted  signals  whose  effect  can  be  represented 
by  an  equivalent  noise  power  whenever  the  probability  of  power- 
independent  distortion  is  negligible.  Further  discussion  of  interference 
from  unwanted  signals  may  be  found  in  annex  V  of  the  report  by  Rice 
et  al.  (1967). 

The  biases  Y_,  Y  ,  and  Y  expressed  in  decibels  are  assumed  to 
T  L  c 

be  normally  distributed  in  time  and  with  location,  with  errors  of  pre¬ 
diction  also  normally  distributed.  The  biases  may  then  be  expressed  in 
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terms  of  a  standard  normal  deviate  z  (q)  ,  and  the  variance  of  the  avail- 

o 

able  signal-to-noise  ratio  in  time,  a^.,  with  location,  a and  with  pre- 
diction  error,  <r  .  The  variance  a“  is  defined  in  terms  of  variances  in 

2  2  T 

time  cr_  and  cr-,  associated  with  the  available  wanted  signal  power  W 
Ta  Tn  a 

and  the  noise  power  W  ,  respectively,  and  of  p^,  the  coefficient  of  cor¬ 
relation  between  them: 


2  2  2 

=  o„  +  at  -  2p 


Ta 


Tn 


T  Ta  Tn 


dB 


(1.  14a) 


where  a  is  defined  by  (1.9)  and  { 1 .  10) . 

2 

Similarly  the  path-to-path  variance  cr^  of  the  available  signal-to- 
noise  ratio  may  be  expressed  as: 


o 


2 

L 


2,2 
crT  +  cr 
La  Ln 


ZpL  ^La  °Ln 


dB2  , 


(1.  14b) 


where  <r.  and  <t  are  the  location-vo-location  standard  deviations  of 
La  Ln 

the  available  wanted  signal  power  W  and  the  noise  power  W  ,  respec- 

&  D 

tively,  and  is  the  normalized  coefficient  of  correlation  between 

them.  For  service  limited  by  a  background  of  man-made  noise,  a  study 

of  available  data  indicates  that  <r_  4  dB  and  cr  ^4  dB,  while  <rT  =* 

1 n  Ln  La 

10  dB»as  previously  stated.  The  correlation  coefficients  p^.  and  p^ 
are  usually  assumed  to  be  positive  because  when  the  received  signal  is 
hi  eh  as  a  result  of  good  propagation  conditions  the  noise  level  is  likely 
to  be  high  also. 

The  estimated  standard  error  of  prediction  n  may  be  defined  in 

2  C 

terms  of  the  variance  a  associated  with  the  received  power  exceeded 

ca 

at  least  a  fraction  q  of  the  time  at  a  fraction  q  of  the  locations,  the 

2 

variance  a  associated  with  the  corresponding  noise  power,  and  the 
c  n 

correlation  p  between  them,  with  additional  allowances  for  time  and 

c 

location  variability: 
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p 


where  c r  =  5fl  +  0.6exp(-d/100)]  dB.  (1.16) 

ca  e 

The  estimated  prediction  error  CTC  includes  an  allowance  for  errors 
in  predicting  the  required  signal-to-noise  ratio  Rr  (gr).  The  coefficients 
0.12  and  4  in  (1.15)  and  5  and  0.6  in  (1.16)  represent  the  best  empirical 
estimates  presently  available. 

Estimates  of  cr  =*  4  dB ,  a  2,1  5  dB  are  probably  adequate,  and  the 
c  n  x 

standard  deviation  cr  is  computed  using  (116)  as  a  function  of  the  ef- 

C  21 

fective  distance  d  ,  defined  by  (1.  3),  Little  information  is  available 
e 

regarding  the  correlation  p  between  prediction  errors  for  wanted 

c 

signals  and  noise.  The  correlation  is  expected  to  be  positive  and  less 
than  unity. 

The  safety  factors  or  biases  shown  in  (1.  13)  may  be  allowed  for 
separately  in  the  order  given.  The  bias  Y,^,  assuming  q^.  =  0.99,  is 
the  difference  between  an  available  signal-to-noise  ratio  R(0.99,  q^)  dB 
exceeded  at  least  99  percent  of  the  time  at  any  location  and  the  median¬ 
time  value  R(0.  5,  q  )  .  Then  the  bias  Y  ,  assuming  q  =0.99  and 

JLj  j— -  Jl 

q  =0.1,  is  the  difference  between  the  ratio  R(0. 99,  0.  1)  exceeded  at 

JLa 

least  99  percent  of  the  time  at  10  percent  of  the  locations  and  the  median¬ 
time,  median-location  value  R(0.5,  0.5): 


yt  =  -  <rT(pT)  *0(q.r) 

dB, 

(1.  17a) 

YL  =  -  'lV  ‘oV 

dB, 

(1.  17b) 

where  p^.  is  the  time  correlation  and  is 

the  path- 

.v,-path  cor rela  - 

tion  between  wanted  signals  and  noise. 

The  bias  or  safety  margin, 

Y  *  -  cr  (p  )  z  (Q) 
c  c  c  o 

dB, 

(1. 17c) 

1-18 
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required  to  achieve  a  given  signai-to-noise  ratio  with  a  service 

probability  Q  depends  in  turn  on  the  correlation  p  assumed  between 

c 

prediction  errors  for  wanted  signals  and  noise  in  (1.  15). 

The  available  wanted  signal  power  is 

W  =  W  +  G  -  L  dBW,  (1.  18) 

a  t  p  b 

where  Wt  is  the  total  radiated  power  in  dBW,  G  is  the  path  antenna 

gain  in  dB,  and  L  is  the  basic  transmission  loss  in  dB.  The  path  an- 
b 

tenna  gain  may  be  expressed  as 

G  =  G  +  G  -  L  dB,  (1.  19) 

p  1  2  gp 

where  G^  and  G^  are  the  free  space  antenna  gains  in  decibels  relative 
to  an  isotropic  radiator.  The  loss  in  path  antenna  gain,  ,  may  be 

approximated  as  follows: 

For  G,  +  G  s  50  dB  L  ±  \  dB ,  (1.20a) 

1  2  8P 

For  50  <  G  +  G^  «  100  dB 

L  *  0.07  exp[0. 055  (G  +  Gj]  dB  .  (1.20b) 

8P  L  1  2  J 

In  most  situations  we  assume  that  free-space  gains  are  realired.  Where 

the  loss  in  gain  may  be  appreciable,  other  methods  for  estimating  G  are 

P 

described  in  the  report  by  Rice  et  al.  (1967)  .  The  approximation  in 

(1.20b)  tends  to  give  substantially  larger  values  of  L  for  large  G  + 

g  P  1 

G,  than  the  methods  reported  by  Rice  et  al.  Experience  with  actual 
communication  links  at  4500  MHz  using  10  m  parabolic  reflectors  sug¬ 
gests  that  (1.20b)  may  provide  more  realistic  estimates  of  link  perform¬ 
ance.  However,  more  Uu-oretical  and  experimental  studies  are  needed 
to  resolve  this  problem. 
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The  available  signal  to  noise  ratio  R  is 


R  =  W  -  W  dB, 
a  n 


(1.21) 


where  W  is  the  total  equivalent  r-f  noise  power  in  dBW,  It  is  often  con- 

n 

venient  to  express  W  and  W  as  power  spectral  densities,  in  dB(W/hHz), 

a  n  ' 

rather  than  expressing  the  total  power  in  the  r-f  passband  of  the  receiver 
in  dBW. 

The  powers,  or  power  spectral  densities,  W  and  W  and  the  basic 

an 

transmission  loss  L  are  assumed  to  be  normally  distributed  in  time 

b 

and  from  path  to  path.  The  median  values  of  W  ,  W  ,  and  L  may 

an  b 

be  denoted  as  W  ,  W  ,  and  L,  .  Then  the  median  available  signal- 
o  no  do 

to-noise  ratio  R  is 
o 


R  =  W  -  W  =WiG-L,-W  dB, 
o  o  no  t  p  bo  no 


(1.22) 


and  the  signal-to-noise  ratio  R  may  be  expressed  in  terms  of  its  median 


value  R  and  the  biases  Y  .  Y,  and  Y  : 
o  T  L  c 


R  =  R  +  Y  +  Y  +  Y  dB, 
o  T  L  c 


(1.23) 


where 


Y  =  -  ir  z  (Q)  dB  . 
c  co 

The  grade  of  service  at  any  receiving  location  is  satisfactory  if  the 

available  wanted  signal-to-noise  ratio  R  exceeds  the  ratio  R  required 

r 

for  satisfactory  service  in  the  presence  of  fine-grained  time  and  space 
variations  of  signals  and  noise,  that  is,  if  R  -  R  >  0  ,  where 


R  -  R  =  R  +  Y_  +  Y  -  R  -  j  z  (Q)  dB  . 
r  o  T  L  r  c  o 


1. 24) 
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Substituting  (1.22)  in  (1.24)  and  defining 

S  =  W  +  G  -R  -  W  +  Y  +  Yt  dB,  (1.25) 

c  t  p  r  no  T  L  '  ' 

we  can  write 

R-R  =  S  -  L,  ~  a  *  (Q)  .  (1.26) 

r  o  bo  co 

The  service  probability  Q  is  calculated  bv  setting  R  -  R^  equal  to 
zero  and  solving  for  Q  in  terms  of  the  error  function: 


a  z  (Q)  =  S  -  L.  dB, 
co  o  bo 


(1.27a) 


Q  =  0. 5  +  0. 5  erf[(So  -  V^)]  ♦ 


(1.27b) 


Or,  for  a 
solved  for  the 


given  value  of  Q,  such  as  =  0.95,  (1.27a)  may  be 
value  of  S  required  to  achieve  this  service  probability. 
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1-4  List  of  Symbols  and  Abbreviations 
In  the  following  list  the  English  alphabet  precedes  the  Greek  alpha¬ 
bet  and  lower  case  letters  precede  upper  case  letters.  In  general, 
upper  case  letters  are  used  for  quantit?.es  expressed  in  decibels. 

a  an  effective  earth's  radius  that  allows  for  average  refraction 

of  radio  rays  near  the  surface  of  the  earth,  (1.2a). 

A  a  predicted  reference  value  of  attenuation  below  free  space, 
c  r 

expressed  in  decibels,  (1.  1). 

A(q^  )  attenuation  below  free  space  not  exceeded  for  at  least  a  frac¬ 
tion  of  all  randomly  chosen  paths  for  which  other  param¬ 
eters, such  as  frequency,  antenna  heights,  and  path  length, 
are  fixed. 

A(0.  5)  long-term  median  value  of  attenuation  below  free  space  for 
a  specified  period  of  time,  climatic  region,  etc.  (1.  1). 
A(0.1),  A(0.9)  attenuation  below  free  space  not  exceeded  for  10 
percent  and  90  percent  of  the  time,  respectively,  (1.4). 
d  great  circle  distance  in  kilometers 

dB  abbreviation  for  decibel. 

d  an  effective  distance  in  kilometers,  defined  by  (1.3). 
e 

d,  the  sum  of  the  smooth-earth  horizon  distances  for  an  effec- 
Lo 

tive  earth's  radius  a  =  9000  km,  (1.2). 

d  the  distance  between  horizons  for  which  diffraction  and  for- 
si 

ward  scatter  losses  are  approximately  equal  over  a  smooth 
earth,  (1.2a). 

erf  x  the  error  function  of  x ,  (1.6),  defined  as  H(x)  * 
erf  *  x  the  inverse  error  function  of  x.  (1.6). 
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f  radio  wave  frequency,  expressed  in  megahertz,  (1.2). 
g(q,  f)  frequency  factors  used  to  adjust  the  predicted  variability  in 
time  at  100  MHz  for  use  at  other  frequencies,  (1.5)  and 
figure  1.4. 

g(0.  1,  f),  g(0. 9,  f)  frequency  factors  used  to  adjust  Y^O.  1)  and 

Y  (0.9).  respectively,  (1.5)  and  figure  1.4. 

G  path  antenna  gain  in  decibels  above  the  unit  gain  of  an  iso- 
P 

tropic  radiator,  (1.19). 

G  ,G  free-space  antenna  gains,  in  decibels  relative  to  an  iso- 
tropic  radiator,  of  the  transmitting  and  receiving  antennas, 
respectively,  (1.19). 

h  ,  h  effective  heights  in  meters  of  the  transmitting  and  re¬ 
el  ei 

ceiving  antennas,  respectively,  (1.2b). 

L  basic  transmission  loss  in  decibels,  (1.  18). 
b 

L  ,  L  (0.5)  long-term  median  value  of  basic  transmission  loss 
bo  b 

in  decibels,  (1.4)  and  (1.22). 

L  (0.  1),  L  (0.9)  basic  transmission  loss  not  exceeded  for  frac- 
b  b 

tions  0.1  and  0.9  of  hourly  medians,  (1.4). 

L  loss  in  path  antenna  gain,  ( 1 .  19)  and  ( 1 .  20) . 
gP 

N  atmospheric  refractivity  at  the  surface  of  the  earth, 
q  a  symbol  that  represents  any  fraction  between  0  and  1,  (1.  6). 
q^  any  desired  fraction  of  all  randomly  chosen  paths  for  which 
such  parameters  as  frequency,  antenna  heights,  and  path 
length  are  fixed,  (l.  12). 
q^  any  desired  fraction  of  time,  (1.7). 

Q  a  symbol  that  represents  service  probability,  (1.17)  and 
(1.27). 
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R  the  available  signal-to-noise  ratio  expressed  in  decibels, 

(1.21). 

R  the  signal-to-noise  ratio  required  to  provide  satisfactory  ser- 
r 

vice  in  the  presence  of  fine-grained  time  and  space  variations 
of  signals  and  noise,  (1.24). 

R  median  value  of  the  available  signal-to-noise  ratio  expressed 
in  decibels,  (1.22). 

R(q^,,  q^,  Q)  the  sign^l-to-noise  ratio  available  for  at  least  a  de¬ 
sired  fraction  q^,  of  the  time  at  a  desired  fraction  q  of 
locations,  with  a  pi  obability  Q,  (1.13). 

R(0.5,  q^)  the  median  time  value  of  the  available  signal-to-noise 
ratio  at  any  location. 

R(0.99,  q^)  the  available  signal-to-noise  ratio  exceeded  at  least 

99  percent  of  the  time  at  any  location. 

R(Q.99,  0.1)  the  available  signal-to-noise  ratio  exceeded  at  least 

99  percent  of  the  time  at  10  percent  of  the  locations. 

S  a  term  defined  by  (1,2  5), 
o 

V(0.  5)  the  difference  in  decibels  between  a  computed  reference 

value  A  and  the  median  attenuation  A(0.  5)  expected  for  a 
c  r 

specified  climate,  season,  time  of  day,  or  desired  group  of 
paths,  ( 1 .  1) . 

V(0.  5,  d  )  the  adjustment  factor  V(0.  5)  as  a  function  of  the  eifec- 

e 

tive  distance  d  ,  figure  1.1. 
e 

W  radio  frequency  signal  power  or  power  spectral  density  that 
would  be  available  from  an  equivalent  loss-free  receiving 
antenna,  ( 1 .  18) . 

W  the  total  equivalent  r-f  noise  power  or  power  spectral  density 
at  the  terminals  of  a  loss -free  receiving  antenna,  including 
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both  internal  system  noise  referred  to  these  terminals  and 
externally  generated  .unwanted  signals  whose  effect  can  be 
repr es ented  by  an  equivalent  noise  power,  ( 1 .  21). 

W  the  total  power  radiated  from  a  transmitting  antenna  in  a 
given  band  of  radio  frequencies,  (1.  18). 

W  the  median  value  of  W  ,  (1.22). 
no  n 

W  the  median  value  of  W  ,  (1.22). 
o  a 

x  a  parameter  defined  as  x  =  d  / 100  used  in  (1.9)  and  ( 1 .  10) . 

e 

Y^  an  allowance  in  decibels  for  prediction  error,  (1.  13)  and 
(1.  17). 

Y  an  allowance  in  decibels  for  random  variations  in  transmis¬ 
sion  loss  from  location  to  location,  (1.  13)  and  (1.  17). 

Y  an  allowance  in  decibels  for  long-term  variability  in  time, 

(1.  13)  and  (1.  17). 

Y(q^)  the  location  variability  for  any  fraction  of  all  randomly 
chosen  paths  for  which  other  parameters, such  as  frequency, 
path  length  and  antenna  heights, are  fixed,  (1.  12a). 

Y(q^)  the  variability  in  time  about  the  long-term  median  exceeded 
at  least  a  fraction  q^.  of  the  time,  (1.  7). 

Y(0.  1),  Y(0.01)  the  difference  between  the  long-term  median  at- 

Y(0.9),  Y(0.99)  tenuation  and  that  not  exceeded  for  fractions  0.  1, 
0.01,  0.9,  and  0 . 99  of  the  time,  respectively,  (1.5)  and  fol¬ 
lowing  discussion. 

Y  (0.1),  Y  (0.9)  values  of  Y(0.1)  and  Y(0.9),  respectively,  ex- 

o  o 

pected  at  100  MHz  in  a  continental  temperate  climate, 
plotted  versus  an  effective  distance  d  ,  figure  1.4  and  (1.4). 
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z  a  stand  : -d  normal  deviate,  (1.6). 
o 

z  (q)  a  standard  normal  deviate  where  the  symbol  q  represents 
o 

any  fraction  between  zero  and  unity,  (1.6). 

z  (qT  )  a  standard  normal  deviate  where  the  symbol  q  represents 
o  L. 

any  desired  fraction  of  all  randomly  chosen  paths  for  which 
such  parameters  as  frequency,  distance, and  antenna  heights 
are  fixed,  (1.  12)  and  (1.  17). 

z  (q^,)  a  standard  normal  deviate  where  the  symbol  q  represents 
any  desired  fraction  of  time,  (1.7)  and  ( 1 .  17). 
z  (Q)  a  standard  normal  deviate  where  the  symbol  Q  represents 
the  desired  service  probability,  (1.  17). 

9  the  angular  distance  at  which  scatter  and  diffraction  trans- 
si 

mission  losses  are  approximately  equal  over  a  smooth  earth 
c  effective  radius  a  =  9000  km,  (1.2). 

p  the  cross -cor relation  coefficient  between  the  variations  of 

c 

the  received  signal  power  exceeded  at  least  a  fraction  q_ 
of  the  time  at  a  fraction  q^  of  locations  and  variations  of 
the  corresponding  noise  power,  (1.  15). 
p  the  coefficient  of  correlation  between  the  location- to- location 

Xmt 

variations  of  the  available  wanted  signal  power  and  the  noise 
power,  ( 1 .  14b) . 

p^.  the  coefficient  of  correlation  in  time  between  the  available 

wanted  signal  power  and  the  noise  power,  (1.  14a). 

a  estimated  prediction  error,  defined  by  (1.  15). 
c 

j  ,  u  the  standard  deviations  of  the  received  wanted  signal 
c  a  c  n 

power  and  noise  power,  respectively,  exceeded  at  least  a 
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fraction  q^,  of  the  time  at  a  fraction  q,  of  the  locations, 

(1,15);  a  is  defined  by  ( 1 .  16),  a  "4  dB  . 
ca  cn 

cr  (p  )a  value  of  cr  computed  using  a  specified  value  of  ,  (1.  17c). 
cr  ,  <r  ^  the  standard  deviation  and  variance,  respectively,  from 
location  to  location  of  the  available  signal-to-noise  ratio, 

(1.  14b), 

<rT  ,  a,  the  location-to-locatior  standard  deviations  of  the  avail- 
L-a  Ln 

able  wanted  signal  power  and  noise  power  respectively, 

(1.14b);  available  data  indicate  crT  =*  10  dB  and  cr  a 

X  3.  X  n 

4  dB  , 

o'  (p  )  the  location-to-location  standard  deviation  of  the  available 
signal-to-noise  ratio  assuming  a  specified  value  of  p  , 

(1.  17b). 

2 

cr^  the  standard  deviation  and  variance,  respectively,  of  the 

available  wanted  signal-to-noise  ratio  in  time,  (1.  14a). 

<r_  ,  cr  the  standard  deviation  in  time  of  the  available  wanted 
Ta  Tn 

signal  power  and  noise  power,  respectively,  (1.  14a);  cr 

i  <1 

is  computed  using  ( 1 . 9)  and  (1.10);  available  data  indicate 

cr  ~  4  dB  . 

T  n 

cr  (0.  1),  c  (0.9)  standard  deviations  representing  the  bi-normal 
id  id 

time  distribution  of  attenuation  relative  to  free  spice,  illus¬ 


trated  in  figure  1.6,  ( 1 .  <t)f  and  ( 1 .  10) . 

o^.(p  )  the  standard  deviation  in  time  of  the  available  signal-to- 

noise  ratio  assuming  a  cpecified  value  of  p ^ ,  (1.  17a). 

a  a  term  used  in  estimating  prediction  error  that  allows  for 
x 

errors  in  predicting  the  required  signal-to-noise  ratio,  a  =» 
5  dB,  (1.15). 
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ANNEX  2 


STUDIES  OF  TERRAIN  PROFILES 
2-1  Introduction 

Radio  transmission  loss  over  irregular  terrain,  for  the  frequencies 
and  distances  considered  in  this  report,  depends  mainly  on  the  profile 
characteristics  of  a  great  circle  path  between  transmitting  and  receiving 
antennas.  Some  allowance  is  made  for  vegetation  and  man-made  clutter, 
while  large  buildings  and  dense  vegetation  are  treated  in  the  same  way 
as  features  of  the  terrain  itself. 

For  point-to-point  transmission  loss  calculations  for  a  given  ter¬ 
rain  profile  and  antenna  locations,  the  parameters  of  interest,  in  order 
of  their  usual  importance,  are  the  sum  of  horizon  ray  elevation  angles 
0  ,  the  effective  antenna  heights  h  ,  the  path  distance  d  ,  the  in- 

terdecile  range  of  terrain  elevations  Ah(d)  ,  the  horizon  distances 
d  ,  and  the  effective  dielectric  constant  t  ,  and  conductivity  a, 

JLj  Ay  Lm 

for  a  terrain  or  sea  surface.  The  definition,  use,  and  importance  of 
these  parameters  are  explained  in  the  body  of  the  report. 

To  obtain  calculated  reference  values  of  propagation  attenuation  for 
specified  sets  of  terrain  profiles  and  antenna  locations,  or  to  obtain  es¬ 
timates  of  variability  from  location  to  location  with  several  or  all  of  the 
above  parameters  fixed,  within  limits,  and  to  estimate  prediction  er¬ 
rors,  statistical  descriptions  of  these  parameters  are  needed,  especially 
estimates  of  median  values  for  commonly  occurring  situations.  Ter¬ 
rain  statistics  were  developed  for  selected  areas  by  reading  a  large 
number  of  terrain  profiles.  Each  profile  is  represented  by  dis¬ 
crete  elevations  at  uniform  distances  of  half  a  kilometer.  The  areas 
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selected  for  terrain  study  include  one  in  the  tree-covered  rolling  hills 
of  northeastern  Ohio,  and  one  in  the  plains  and  foothills  of  Colorado, 
where  an  exhaustive  transmission  loss  measurement  program  was  car¬ 
ried  out  (Miles  and  Barsis,  1966).  About  100  paths,  each  60  km  in 
length,  were  selected  at  random  throughout  the  continental  United  States 
to  provide  an  estimate  of  "average"  terrain  statistics.  Two  limited  re¬ 
gions,  one  entirely  in  the  mountains  crossing  the  Continental  Divide  and 
the  other  entirely  in  the  plains,  were  chosen  for  an  intensive  study  of 
correlation  along  and  between  profiles. 

Within  each  region  selected  for  intensive  study,  36  profiles,  60  km 
in  length,  were  read  in  each  of  six  directions  separated  in  azimuth  by 
30°,  providing  a  total  of  216  profiles  that  form  a  rather  closely  spaced 
"grid"  over  a  100  km  square  area.  The  101  "random"  paths  throughout 
the  continental  United  States  were  so  chosen  that  they  would  not  ap¬ 
proach  or  cross  each  other.  The  separation  between  adjacent  paths 
ranges  from  60  to  320  km  with  a  median  separation  of  about  200  km  . 

These  "random"  paths  lie  chiefly  in  four  directions,  N-S,  E-W,  NE-SW, 
and  NW-SE.  None  are  located  in  Colorado  which  is  well  represented  by 
the  "plains"  and  "mountains"  grids  and  the  area  over  which  radio  propa¬ 
gation  measurements  were  made. 

In  the  measurement  program  antenna  sites  were  located  randomly 
on  or  near  roads  without  regard  to  the  proximity  of  natural  or  man-made 
obstacles.  To  correspond  with  these  measurements  the  path  profiles 
were  also  selected  arbitrarily,  without  regard  to  the  location  of  hill.)  or 
other  obstacles.  Each  of  the  60  km  profiles  in  the  mountains,  plains, 
and  random  U.  S.  samples  was  considered  in  lengths  of  5,  10,  20,  30,  40, 
50,  and  60  km,  starting  from  one  end  of  each  profile,  to  study  the  effect 
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of  path  length  on  various  terrain  parameters.  A  further  statistical  ter¬ 
rain  study,  assuming  antenna  sites  placed  advantageously  on  hilltops, 
should  be  carried  out  to  show  the  improvement  in  radio  transmission  to 
be  expected  by  increasing  the  effective  height  of  one  or  both  antennas. 

2-2  The  Terrain  Parameter  Ah 

The  path  profiles  described  above  were  used  to  obtain  estimates  of 
several  terrain  parameters  for  paths  of  a  given  length.  Of  these  the 
interdecile  range  Ah(d)  of  terrain  heights,  above  and  below  a  straight 
line  fitted  by  least  squares  to  elevations  above  Bea  level,  was  calculated 
at  fixed  distances.  Usually  the  median  values  of  Ah(d)  ,  for  a  specified 
group  of  profiles,  increase  with  path  length  to  an  asymptotic  value,  Ah  . 
As  explained  in  detail  in  the  body  of  the  report,  this  parameter  Ah  is 
used  to  characterize  nondeterministic  aspects  of  terrain  irregularity. 

Figures  2.1  and  2.2  show  cumulative  distributions  of  Ah(d)  for  the 
random  U.S.  mountain,and  plains  paths.  Each  distribution  represents 
101  profiles  for  random  terrain  and  216  profiles  each  for  mountain  and 
plains  terrain.  These  distributions  show  a  consistent  increase  in  the 
median  value  of  A  h(d)  with  increasing  path  length.  As  one  would  expect, 
the  variance  of  A  h(d)  for  the  randomly  chosen  paths  is  much  greater 
than  that  for  the  more  homogeneous  terrain  included  in  the  plains  and 
mountain  areas. 

Distributions  of  Ah(d)  for  the  area  in  Colorado  where  radio  meas¬ 
urements  were  made  are  shown  in  figure  2.3.  These  paths  are  locat  'd 
in  an  area  that  includes  plains,  foot. nils,  and  mountains,  but  were  con¬ 
sidered  in  only  two  groups.  For  these  much  smaller  groups  of  profiles, 
the  same  trends  are  observed,  the  median  value  of  Ah(d)  increasing  con¬ 
sistently  with  increasing  path  length.  This  trend  was  not  observed  in  the 
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CUMULATIVE  DISTRIBUTIONS  OF  Ah  ;) 
U.  S.  Random  Paths 


Ah(dj  In  AActers 


CUMULATIVE  DISTRIBUTIONS  OF  Ah(d) 


Percent  of  Profiles 
Figure  2.2 
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Ah(d)  in  Meters 
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CUMULATIVE  DISTRIBUTIONS  OF  Ah(d) 
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CUMULATIVE  DISTRIBUTIONS  OF  Ah (d) 

Nortiaasfsrn  Ohio 


FVcent  of  Profiles 
Figure  2  4b 
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distribution  of  Ah(d)  for  northeastern  Ohio,  as  shown  in  figure  2.  4a, 
where  for  path  lengths  from  20  to  50  km  no  consistent  ase  in  Ah(d) 

is  shown.  These  data  include  a  group  of  profiles  that  form  radials  from 
the  location  of  the  "central  transmitter",  and  a  smaller  group  of  pro¬ 
files  that  radiate  from  each  of  5  peripheral  transmitter  locations.  When 
the  profiles  from  each  of  the  peripheral  transmitter  locations  are 
grouped  together,  the  usual  trend  is  observed,with  the  median  Ah(d)  in¬ 
creasing  consistently  with  path  length,  as  shown  in  figure  2.4b.  For  the 
groups  of  paths  from  the  location  of  the  central  transmitter  the  median 
Ah(d)  is  higher  at  all  distances  but  particularly  for  the  10-and20-km  pro¬ 
files.  This  indicates  that  the  terrain  is  more  rugged  in  the  vicinity  of 
the  central  transmitter  than  in  the  remainder  of  the  area. 

Median  values  of  Ah{d)  are  plotted  versus  distance  in  figures  2.  5, 

2.6,  and  2.7  for  random,  plains,  and  mountain  paths.  For  each  median 
Ah(d)  ,  an  estimate  of  the  asymptotic  value  Ah  was  calculated  using  (3) 
of  the  main  body  of  this  report: 

Ah(d)  =  Ah  [l  -  0.  8  exp(-  0.  02  d)]  m,  (3) 

where  Ah(d)  and  Ah  are  in  meters  and  the  distance  d  is  in  kilometers. 
Choosing  appropriate  values  of  Ah  from  these  calculations,  the  smooth 
curves  on  the  figures  were  plotted.  Figure  2.  5  shows  for  the  101  ran¬ 
domly  chosen  paths  how  the  medians  of  data  at  each  distance  compare 
with  a  curve  computed  assuming  the  asymptotic  value,  Ah  =  90  m  .  In 
the  plains  and  mountain  areas  36  profiles,  60  km  in  length,  were  read  in 
each  of  6  directions.  To  determine  whether  the  terrain  changes  in  a 
predictable  way  from  one  direction  to  another  the  paths  in  each  direc- 
tion  were  considered  separately.  In  figures  2.6  and  2.7  each  symbol 
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THE  PARAMETER  Ah(d>  VERSUS  DISTANCE 

U  S.  Plains  Paths 


Distance  in  Kilometers 
Figure  2  6 


PARAMETER  Ah(d)  VERSUS  Di STANCE 

U  S.  Mortar  Raths 


Distance  in  Kilometers 
Figure  2.7 


p 


represents  the  median  value  of  Ah(d)  for  36  paths  in  a  specified  direc¬ 
tion.  Figure  2,6  shows  that  paths  in  an  east-west  direction  consistently 
show  a  larger  interdecile  range  of  terrain  heights  than  those  in  a  north- 
south  direction.  Data  from  the  E-VV  paths  suggest  an  asymptotic  value 
Ah  =  80  m  ,  while  those  from  the  N-S  paths  suggest  Ah  =  50  m.  The 
median  value  for  paths  in  all  six  directions  is  Ah  =  62  m.  In  the  figure 
smooth  curves  calculated  from  (3)  for  Ah  -  50,  62  and  80  m  are  com¬ 
pared  with  medians  from  profiles  in  each  direction.  The  data  for  the 
mountain  \  ofiles  do  not  show  as  consistent  a  directional  trend, but  in 
this  case  the  east-west  paths  show  the  smallest  values  and  the  east- 
northeast  paths  tend  to  yield  the  largest  values  of  Ah(d)  as  shown  on 
figure  2.7.  The  data  for  the  most  part  fall  between  the  limits  Ah  =  800m 
and  Ah  =  1000  m,  nd  the  terrain  is  represented  by  the  asymptotic  value 
Ah  =  900  m. 

Figure  2.8  shows  median  values  of  Ah(d)  tor  the  few  Colorado 
mountain  paths  over  which  measurements  were  made,  and  a  curve  of 
Ah(d)  for  Ah  =  6  50  m  .  The  value  Ah  =  650  m  was  assumed  in  calculat¬ 
ing  transmission  loss  to  be  expected  over  these  paths.  For  the  shorter 
paths  this  value  is  too  large,as  indicated  in  the.  figure.  For  the  paths  in 
Ohio  and  those  in  the  Colorado  plains,  the  value  Ah  =  90  m  was  as¬ 
sumed  in  calculating  expected  transmission  losses.  Figure  2.9  shows 
the  corresponding  smooth  curve  of  Ah(d)  versus  d  compared  with 
median  values  of  Ah(d)  from  pat''  profiles.  The  symbol  x,  represen. ing 
median  values  of  Ah(d)  for  Colorado  plains  paths,  shows  good  a  greement 
with  the  curve  especially  at  the  greater  distances  where  most  of  the  radio 
data  were  recorded.  The  square  symbol  represents  median  values  of 
Ah{d)  for  all  Ohio  paths  at  each  distance,  the  circle  represents  values 
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THE  PARAMETER  Ah(d)  VERSUS  DISTANCE 

N.  E.  Ohio  and  Colorado  Plains 
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Figure  2.9 


for  the  central  transmitter,  and  the  triangle  those  for  the  peripheral 
transmitters  at  each  distance.  Although  the  terrain  in  this  Ohio  area  is 
obviously  not  homogeneous,  the  value  Ah  =  90  m  was  chosen  as  being 
■  r esentative  of  the  terrain, 

he  following  values  of  Ah  were  chosen  as  representative  of  the 
terrain  in  each  of  the  areas  for  which  terrain  statistics  were  obtained: 


Area 

Ah  in  m 

U .  S.  random 

90 

Plains  grid 

62 

Mountain  grid 

900 

Colorado  plains  (meas.  paths) 

90 

Colorado  mountains  (meas.  paths) 

650 

NE  Ohio  (meas.  paths) 

90 

These  figures  provided  most  of  the  basis  for  the  values  of  Ah  listed 
in  table  1  of  the  main  body  of  the  report. 

2-3  The  Horizon  Distance,  d^ 

When  a  detailed  terrain  profile  is  available  for  a  given  path,  the 
horizon  distances  d^  ^  and  d^  an<3  their  sum  d^  may  be  obtained  di¬ 
rectly  from  the  profile  information  and  an  estimate  of  the  effective 
earth'f  radius.  When  individual  path  profiles  are  not  available,  median 
values  of  d^  and  are  estimated  as  functions  of  the  median  effec¬ 

tive  antenna  heigh'. s,  the  terrain  irregularity  factor  Ah,  and  the  corres¬ 
ponding  distances  to  the  horizon  over  a  smooth  earth. 
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The  estimates  d  and  d  computed  using  equation  (5c)  of  the 

Jm  i  -L  L. 

main  body  of  the  report, 

d  =  d  exp(-  0.07  \/Ah/he)  km,  (5c) 

LI,  2  Lsi ,  2  ’ 

approach  the  smooth-earth  values  dT  and  d.  as  the  terrain  factor 
Ah  approaches  zero  or  as  the  antenna  heights  become  very  large.  To 
determine  the  constant  in  (5c),  horizon  distances  were  obtained  for  the 
sets  of  U.  S.  random,  plains,  and  mountains  profiles  with  pairs  of  an¬ 
tenna  heights  above  ground,  h  and  h  ,  chosen  as  follows: 

gi  gi 


h 

h 

h 

h 

h 

gl 

g* 

gi 

gi 

gi 

gi 

1 

1 

3 

3 

10 

10 

1 

10 

3 

10 

30 

30 

3  30 


For  paths  of  lengths  30,  40,  50,  and  60  km,  the  horizon  distances  d^  ^  , 

d^7,  and  the  sum  of  the  horizon  distances  dT  were  obtained  from  each 

profile.  Figures  Z.  10,  Z.  11,  and  2.  12  show  cumulative  distributions  of 

dj^,  d  and  d^  .  As  shown  in  table  2.  1,  these  horizon  distances  are 

independent  of  path  length,  provided  that  the  path  length  chosen  is  greater 

than  d^  .  Figure  2.10  shows  cumulative  distributions  of  d  ,  d^7,  anc* 

d  for  10 1  random  paths  with  h  =  h  =1  m.  Figures  2.  11  and  2.  12 
L  gi  gi 

show  similar  distributions  for  216  paths  in  the  plains  and  in  the  rugged 

mountains  with  h  =  h  =  10  m  .  Note  that  the  median  value  of  d,  is 
gi  gi  L 

always  greater  than  the  sum  of  the  medians  of  d  an  !  d  . 

i-  1  Lb. 
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Horizon  01  stance  in  Kilometers 
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CUMULATIVE  DISTRIBUTIONS  OF  HORIZON  DISTANCES 
U.  S.  Random  Paths 


^rcsnt  of  Profiles 
Figure  2.10 
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CUMULATIVE  DISTRIBUTIONS  OF  HORIZON  DISTANCES 
U.  S.  Ploln»  Path* 
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CUMULATIVE  DISTRIBUTIONS  OF  HOR.ZON  DISTANCES 
U.  S.  Mountain  Pott)* 
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Median  Values  of  Horizon  Distances 
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The  constant  in  (5c)  was  estimated  by  computing  the  smooth  earth 

distances  d  and  d  that  correspond  to  these  median  values  of 
L<Sl  Lj  82 

d  .  and  d  and  rewriting  (5c)  as 

Li  1  Lj£ 

dLl,2  =  dLB,,2  (2.1a) 

kl.  2  =  ^  (dLsi  ,  ,  /  ' dul.  2>  ^  h^'UZ-  ^  lb> 

The  median  of  all  computed  values  was  k  =  0.  07  . 

Table  2.2  shows  calculated  values  of  d  and  d  ,  for  four  an- 

Lsi  LI 

tenna  heights  for  the  random,  plains,  and  mountains  paths.  Table  2.  3 
shows  a)  median  values  of  d  from  profiles  and  b)  corresponding  cal- 
culated  values  for  seven  antenna  height  combinations. 

The  values  show  rather  good  agreement  in  all  areas  for  low  to 
medium  antenna  heights,  with  a  decided  tendency  to  overestimate  d 

L 

when  both  antennas  are  as  much  as  30  m  above  ground,  especially  in 
the  rr  mntains.  These  comparisons  are  all  made  assuming  effective  an¬ 
tenna  heights  equal  to  structural  heights. 

2-4.  The  Elevation  Angle  8 

For  each  of  the  large  number  of  terrain  profiles  available,  the  ele¬ 
vation  angles  9  and  9^  were  computed  using  (3.  1)  of  annex  3 , and  for 
each  profile  the  sum  of  these  angles  9  was  also  computed.  Figures 

2.  13,  2.  14,  and  2.  15  Bhow  cumulative  distributions  of  8  ,8  ,  and 

ei  ez 

8^  given  in  milliradians .  These  are  plotted  on  logarithmic  probability 
paper,  but  it  should  be  kept  in  mind  that  a  small  percentage  of  these 
angles  are  negative,  especially  in  the  group  of  plains  paths.  As  pre¬ 
viously  observed  for  the  horizon  distances  d  ,  d  ,  and  d  ,  these 
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Calculated  values  of  d  „  and  d,  , 

S  i  i-  1 
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el 
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"  r  1 1,1  ’?■"  7  - 

10 

30  m. 

Random  Paths,  Ah  = 

90  m,  N 

=  310,  a 
s 

=  8640  km. 

d 

4.  16 

7.  20 

13.15 

22.  77 

Lsi 

d 

3.  1 

5.  6 

10.  7 

20. 2 

LI 

Plains  Paths,  Ah  =  62 

=  290,  a  =• 

8 330  km  . 

d 

4.08 

7.07 

12.  91 

22.  36 

Lsi 

d 

3.2 

5.  5 

10.8 

20.2 

LI 

Mountain  Paths,  Ah 

=  900  m, 

N  =  2  50  , 
s 

a  =  7  8  50  km  . 

d 

3.  96 

6.  86 

12.  53 

21. 70 

Lsi 

d 

1 . 6 

2.  7 

6.  4 

14.8 

LI 

*d  =  dr  exp{  -  0.07  *V"Ah/  he) 
LI  Lsi 
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Table  2.  3 


The  Sum  of  the  Horizon  Distances  dT 


Random  Paths 


a) 

11.  2 

13.  5 

15.  2 

17.  5 

19.  7 

24.  2 

35.  2 

b) 

6.  2 

11.  2 

13.  8 

16.  3 

21.4 

25.  8 

40.  4 

Plains  Paths 

a) 

12.  5 

15.2 

17.  2 

18.  7 

21.  4 

25.  5 

35.  4 

b) 

6.  4 

11.0 

14.  0 

16.  3 

21.  6 

25.  7 

40.  4 

Mountain  Paths 

a) 

9.2 

9.2 

10.  2 

10.  2 

10.  7 

11.2 

13.  5 

b) 

3.  2 

5.  4 

8.  0 

9.  1 

12.  8 

17.  5 

29.  6 

a)  Median  value  from  profiles 

b)  Calculated  using  (5e) 
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elevation  angles  are  independent  of  path  length,  provided  that  the  path 

length  chosen  is  greater  than  the  sum  of  the  horizon  distances. 

Figure  2.13  shows  cumulative  distributions  8,0  ,  and  0  for 

ei  ei  e 

101  random  paths  with  h  =  1  m,  and  figures  2.  14  and  2,15  show 

gi,* 

similar  distributions  for  216  paths  in  the  plains  and  in  the  mountains, 

with  h  =  10  m.  As  with  the  horizon  distances,  the  median  value  of 
gii  2 

0  is  always  greater  than  the  sum  of  the  median  values  of  8  and  9 
e  ei  e2 

Table  2.4  shows  median  values  of  8  ,8  and  8  for  several  antenna 

ei  ei  e 

height  comomations . 

Estimates  of  the  elevation  angles  were  made  usin^  he  expression 
(6a)  given  in  the  main  body  of  this  report. 


0 

ei ,  2 


0.0005  T  „  (dLsi  i 

'  d.  1,3  l  d.  * 

Lai,j  L  \  -i-l,2 


■) 


^h 


radians  .  (6a) 


Table  2.5  shows  median  values  of  0  ,0  ,  and  8  calculated  using 

ei  ei  e 

this  expression  and  actual  values  of  d  and  d  f  as  read  from  the  pro- 

i-  i.  i-  «— 

files.  Since  no  consistent  distance  dependence  is  observed,  median 
values  of  8  for  each  antenn3  height  combination  are  compared  with  the 
calculated  values.  For  the  random  and  plains  paths,  the  median  esti¬ 
mates  of  8  are  approximately  equal  to  the  median  values  given  in 

Table  2.4.  For  the  mountain  paths  the  eetimated  values  of  0  are  al- 

e 

most  twice  as  large  as  those  obtained  directly  from  profiles. 

Considering  the  situation  where  predictions  are  based  solely  on  es¬ 
timates  of  terrain  characteristics,  without  actual  profiles  for  individual 

paths,  values  of  0  ,0  ,  and  0  were  calculated  using  equation  (6a) 

r  ei  ei  e 
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CUMULATIVE  DISTRIBUTIONS  OF  ELEVATION  ANGLES 
U,  S.  Plains  Palhs 


Table  2.  4 


Median  Values  of  the  Elevation  Angles 


0,8,  and  9  in  Milliradians 
ei  e  2  e 


h  1,1  3,  3  m 

gii  * 


a,  0  0  9  0  0  9 

km  ei  e2  e  el  ti  e 


Random 


30 

2.95 

10. 72 

19.  73 

2.  59 

9.  33 

17.  63 

40 

3.84 

12. 64 

25.  10 

2.85 

11.00 

21.  38 

50 

3.  98 

13.  11 

22.  52 

3.25 

12.  78 

19.  34 

60 

3.  96 

1  5.  36 

30.  38 

3.  30 

13.  71 

29.  73 

Plains 

30 

3.  91 

7.  30 

13.  51 

2.  62 

6.  20 

11.  72 

40 

3.  67 

8.  58 

1  5.  24 

2.  68 

8.06 

13.  04 

50 

3.  51 

10.46 

15.  50 

2.  59 

9.  18 

13.  68 

60 

3.  69 

11.07 

1  5.  94 

2.71 

10.  17 

14.  19 

Mountains 

30 

108.4 

124.  8 

268.  1 

106.  1 

122.  6 

263.  1 

40 

105.  0 

81.8 

212.  2 

104. 4 

81.  3 

207.  3 

50 

99.  5 

75.  7 

206.  5 

99.  1 

73.  9 

203.  4 

oO 

100. 0 

8  5.  0 

227.  9 

97.  1 

84.  4 

225.  9 
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Table  2.4  (continued) 


h  3,  10  10,  10  m 

g‘.  * 


d,  e  e  e  e  e  e 

km  ei  ej  e  ei  ez  e 


Rand*.  m 


30 

2. 

68 

6. 

Oi 

15.  C 4 

1. 

59 

7. 

26 

13. 

22 

40 

3. 

42 

7. 

54 

20.04 

1. 

67 

7. 

70 

19. 

72 

50 

3 . 

28 

10. 

53 

16.44 

1. 

83 

10. 

,  86 

14. 

79 

60 

3. 

40 

11. 

26 

25.  22 

1 

33 

i  X  . 

26 

23. 

16 

Plains 


30 

2.  70 

4. 40 

9.  86 

1.  37 

4  44 

8.47 

4C 

2.  70 

6.25 

10. 90 

1.47 

6.  26 

9.  75 

50 

2.  59 

7.  78 

1L.-8 

1.  58 

7.  78 

9.  80 

60 

2.71 

8.  62 

12.22 

1.  26 

8. 62 

10.  62 

Mountains 


30 

106.  1 

117.  1 

255.3 

100. 9 

117. 

249.  0 

40 

104.4 

77.  1 

201.  7 

10-  .  2 

77.  1 

197.  9 

50 

99.  1 

68.  3 

199.  3 

94.  9 

68.  3 

192.  1 

60 

97.  3 

79.  0 

220.  1 

9  3.  2 

79.  0 

21^.4 
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Table  2.  5 

Median  Values  of  the  Elevation  Angles  Calculated 
Using  (6a),  With  d^  ,  ^rom  Pr°file3 

9,9  ,  9  in  milliradians 

ei  e2  e 


h  1,1  3,  3  m 

gi.  2 


d,  0  9  9  9  0  0 

km  ei  e2  e  ei  ez  e 


Random 


30 

0.  11 

8.88 

24.  17 

0.81 

7.  78 

25.  46 

40 

0 

14.  73 

33.  71 

0.06 

10.  56 

25.  63 

50 

0 

4.98 

33.  36 

0.06 

5.  69 

23.  79 

60 

-2.  17 

8.88 

30.01 

0.06 

7.  78 

26.  54 

Plains 

30 

1.  15 

0 

11.  51 

1.  51 

1.  51 

10.  10 

40 

1.  15 

-0.  14 

20.09 

0.  78 

i.  51 

8.  78 

50 

1.  15 

0 

17.  65 

0.17 

0.  78 

9.  93 

60 

1.  15 

-1.11 

14,  28 

0.  17 

0 

9.  65 

Mountains 

30 

46.  8 

241.8 

678.6 

108. 9 

303.  9 

754.0 

40 

46.8 

85.  8 

359.  6 

108. 9 

147.  9 

479.  6 

50 

46.  8 

35.8 

418.  6 

108. 9 

147. 9 

510.  2 

o0 

46.  8 

144.  3 

581.  1 

108. 9 

206.  4 

705.  2 
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Table  2.5  (continued) 


h  3,  10  10,  10  m 

g».  2 


d,  0  0  0  0  9  9 

km  el  ez  e  ei  e2  e 


Random 


30 

1.  26 

1.  36 

15.  20 

0.  54 

1.  36 

17.  99 

40 

0.  44 

3. 04 

23.  36 

0.20 

3.  04 

20.68 

50 

0.06 

3.  79 

15.  85 

0 

3.  79 

14.  55 

60 

0.06 

4.  68 

15.  59 

0 

4 .  6b 

14.  74 

Plains 

30 

1.  96 

0 

5.  96 

0.  37 

0 

4.47 

40 

0.  78 

0 

4.  99 

0.07 

0 

2.  68 

50 

0.  17 

-0.  19 

4.  77 

-0.  19 

-0.  19 

3.  12 

60 

0.  17 

0.07 

5.  72 

-0.  19 

0.07 

2.03 

Mountains 

30 

108.9 

244.  2 

651.  1 

118.8 

244.2 

618. 4 

40 

106.  9 

146.  7 

504.  7 

118.8 

146.  7 

526.  8 

50 

108.9 

146.  7 

495.  6 

118.8 

146.  7 

520. 9 

60 

108.  9 

146.  7 

580.  6 

118.8 

146.  7 

544.  1 
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with  values  of  d  and  d  calculated  as  functions  of  Ah  (5c),  As  one 

Lj  -L 

would  expect,  these  estimates  of  6  do  not  correspond  closely  to  me¬ 
dian  values  from  individual  profiles. 

Table  2.6  shows  the  comparison  between  (a)  median  values  of  9 
obtained  directly  from  terrain  profiles,  (b)  those  estimated  using  actual 
horizon  distances  in  (6a),  and  (c)  values  calculated  using  (6a)  and  esti¬ 
mates  of  and  dT  y  computed  using  (5c).  For  the  random  and  plains 


Table  <1. 6 


Elevation 

Angles  9 

in  Milliradians 

h 

1,  1 

i,  3 

3.  10 

10,  10 

gl  .  2 

Random 

a) 

23.  8 

20.  5 

18.  2 

17.  2 

b) 

31.7 

2  5.6 

13.  7 

16.  4 

0 

8.  7 

3.9 

1.  5 

-1.0 

Plains 

a) 

15.4 

13.  4 

10.  9 

9.  8 

b) 

16.  0 

9.  8 

5.  4 

2.  9 

0 

4.  5 

1.  5 

-0.  2 

-1.  9 

Mountains 

a) 

220.  0 

218.  6 

210.  9 

20o.  6 

b) 

49'"1.  8 

607.7 

542.  6 

5  3  5.4 

c) 

459.  2 

264.  0 

174.  4 

84.  8 

a) 

b) 


c) 


median  values  of  9  from  profiles 

e 

median  values  of  0e  calculated  using  (ba)  and  values  of 
d  ,  d  from  profiles 

l_j  1  Lit, 


values  of  0e  calculated  using  (6a)  with  values  of 
calculated  using  (A  ) 


2-  V 


paths,  method  (c)  uuderestimates  9  for  all  antenna  height  combinations 
tested,  but  for  the  mountain  paths  it  provides  a  better  estimate  than 
method  (b),  which  overestimates  the  elevation  ai  by  a  wide  margin, 
yielding  values  more  than  twice  as  large  as  the  median  values  from  ter¬ 
rain  profiles.  Since  the  elevation  angle  0  is  probably  the  single  most 

e 

important  terrain  parameter  for  predicting  transmission  loss,  this  wide 
range  in  estimates  introduces  considerable  prediction  error. 


2-5  Terrain  Parameters  for  Colorado  Plains, 

Mountains, and  Northeastern  Ohio 

Profiles  of  the  paths  in  the  foothills  and  plains  of  Colorado  and  in 
northeastern  Ohio  were  used  to  obtain  horizon  distances  and  elevation 
angles.  These  are  compared  with  values  calculated  as  previously  des¬ 
cribed  in  sections  2-3  and  2-4  for  the  U.  S.  random,  plains,  and  moun¬ 
tain  paths. 

These  paths  represent  much  smaller  samples  than  those  previously 
discussed,  so  one  would  expect  the  terrain  statistics  to  be  less  consist¬ 
ent.  Table  2.7  shows  median  values  of  the  horizon  distances  d^  , 
and  d^  for  the  antenna  height  combinations  used  in  the  measurement 
program,  at  distances  of  20,  30,  50  and  80  km.  The  total  number  of 
paths,  N,  in  each  group  is  also  tabulated.  The  horizon  distances  do  not 
appear  to  be  independent  of  path  length  as  they  were  for  the  U.  S.  random, 
plains,  and  mountain  paths.  This  may  be  attributed  in  part  to  the  small 
sample  size,  and  in  Colorado  to  the  fact  that  all  paths  radiate  from  a 
single  transmitter  located  in  the  plains,  resulting  ir.  median  values  of 

d  much  greater  than  those  of  d  ,  even  with  comparable  antenna 
LI  L  l. 
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Table  2.  7 

The  Horizon  Distances,  Colorado  Plains,  Mountains, and  Ohio 
Median  Values  front  Profiles 


h 

g» 

4 

h  l 
g* 

0.6 

l 

.  7 

3 

6 

9  m 

dkm 

N 

dLl 

dL2 

dL 

Tz 

dL 

dL2 

dL 

dL2 

dL 

dL2 

dL 

Colorado  Plains, 

Ah  = 

:90m, 

N  = 
s 

290, 

a  =  8  327  km 

20 

14 

11.2 

2.0 

19.  1 

2.  7 

18.  1 

4.  7 

19.2 

7.  2 

19.4 

7.  2 

19.4 

30 

33 

12.  6 

2.0 

21.6 

2.0 

21.6 

2.  0 

21.6 

2.0 

20.4 

4.  4 

2  3.4 

50 

43 

14.  4 

5.0 

o 

oo 

^  J 

6.0 

29. 0 

6. 0 

29.  6 

9.  2 

32.  3 

18.0 

33.  7 

80 

52 

21.0 

4.  5 

28.  5 

10.8 

33.  3 

16.  5 

35.  8 

Colorado  Mountains, 

Ah  =  650 

m , 

N  - 
s 

290  , 

* 

II 

OO 

ru 

7  km 

20 

10 

7.  6 

1.2 

10.  5 

1.  2 

10. 9 

1.2 

10. 9 

2.  0 

10. 0 

2.0 

10. 9 

30 

14 

9.  8 

1.  5 

11.  5 

1.  5 

11.  5 

1.  5 

11.  5 

1.  5 

11.  5 

1  5 

11.5 

50 

16 

10.  2 

2.8 

15.0 

2.8 

15.0 

2.8 

15.0 

4.  3 

15.  3 

4.  3 

15.  3 

NE  Ohio, 

Ah  =  90  m, 

N  =  3i2 

3 

,  a 

=  867  6  km 

20 

42 

5.0 

?  ) 

4_  ,  w 

14.  i 

2.8 

12.  3 

3.  0 

12.  6 

5.  0 

15.  0 

30 

o2 

7.0 

2.  5 

18.0 

2.8 

18.0 

5.  t; 

23.  2 

6.  0 

2  3.:. 

50 

92 

8. 6 

3.  8 

22.  3 

4.  4 

23.  1 

5.  3 

25.  3 

7.  9 

26.  3 
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heights.  The  median  value  of  d^  is  not  always  much  greater  than  the 
sum  of  the  medians  of  d  and  d  especially  with  the  higher  receiver 
heights.  The  paths  in  Ohio  radiate  from  six  different  antenna  locations 
and  show  median  values  of  d  much  greater  than  the  sum  of  the  medians 

■L< 

of  dL1  aad  dL2  . 

The  smooth-earth  distances  dT  and  dT  were  calculated  for 

L  s  1  Lsz 

each  antenna  height  and  used  in  (5c)  to  calculate  values  of  dT  ,  d^  , 
and  d^ : 


d.  ,  ,  =  d  ex 

L 1,  2  I. si ,  z 


p(-  0.07  h/he)  . 


(5c) 


Table  2.8  shows  these  calculated  values  compared  with  median  values 
from  the  profiles.  The  calculated  values  of  d^  are  consistently  less 
than  the  median  values  from  terrain  profiles,  es pecially  in  the  Colorado 
plains.  The  calculated  values  of  d^  *  however,  correspond  quite  well 
with  those  from  the  profiles.  This  better  agreement  with  d^  is  to  be 
expected  as  the  receivers  are  much  more  randomly  located  th«*n  the 
transmitters.  The  calculated  values  of  d^  and  d^»  using  (5c), never 
exceed  the  corresponding  smooth-earth  values,  but  for  the  Colorado 
plains  paths  the  median  d^  is  always  much  greater  than  the  corres¬ 
ponding  value  of  d  .  This  results  xrom  the  fact  that  the  transmitter 

X-  6 

is  located  in  a  bowl  or  depression  with  rising  ground  in  all  directions  so 
the  horizon  is  much  farther  away  than  it  would  be  due  to  the  normal  fall- 
off  of  the  earth. 

For  each  of  the  terrain  profiles  in  the  Colorado  plains  and  foothills 

and  in  northeastern  Ohio,  the  elevation  angles  0,0  ,  and  their  sum 

e  l  ez 

0  were  computed  using  (3.  1)  from  annex  3  of  this  *  eport.  Table  2. 9 
shows  median  values  at  each  distance  for  the  antenna  heights  used  in  the 
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Table  2.8 

The  Horizon  Distances,  Colorado  Plains,  Mountains, and  Ohio 
Median  and  Calculated  Values 


h 

gi 

-4m,  h 

0.6 

g2 

1.7 

3 

6 

9  m 

dn 

dL2 

dL 

dL2 

dL 

dL2 

dL 

dL2 

dL 

dL2 

dL  km 

Colorado 

Plains,  Ah 

=  90  m,  N 

3 

O 

O' 

fxj 

II 

a  = 

8  327  km 

a) 

13.  8 

3.  5 

2  5.7 

4.0 

25.  3 

4.0 

28.  5 

5.  6 

32.  3 

11.  2 

33.  7 

b) 

6.  1 

2.4 

8.  5 

4.0 

10.  1 

5.  2 

11.3 

7.  6 

13.  7 

9.8 

15.  9 

c) 

8.  2 

3.  2 

11.4 

5.  3 

13.  5 

7.  1 

15.  2 

10. 0 

18.  2 

12.  2 

20.4 

Colorado 

Mountains, 

Ah  =  i 

6  50  m , 

N  = 
s 

290, 

a  =  8  32  7  km 

a) 

9.8 

1.  5 

11.  5 

1.  5 

11.  5 

1.  5 

11.  5 

2.0 

11.  5 

2.  0 

11.  5 

b) 

3.  7 

1.4 

5.  1 

2.4 

6.  1 

3.  2 

6.  9 

5.  0 

8.  7 

6.  8 

10.  5 

c) 

8.  2 

3.  2 

ll.i 

5.  3 

13.  5 

7.  1 

15.  2 

10.0 

18.  2 

12.  2 

20.4 

NE  Ohio, 

Ah  = 

90  m , 

N  = 
s 

312,  a 

=  8676  km 

a) 

7.0 

2.  5 

18.0 

2.  8 

18. 0 

5.  3 

23.  2 

6.  C 

23.  5 

b) 

6.  2 

3.  2 

9.  4 

5.  4 

11.6 

7.8 

14.  0 

10.0 

16.  2 

c) 

8.  3 

4.  2 

12.  5 

7.  2 

15.  5 

10.  2 

18.  5 

12.  5 

20.8 

a)  Median  values  from  profiles 

b)  Calculated  values  using  (5c) 

c)  Smooth  earth  values,  dT  and  d 

Lsi ,  J 
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Table  2.  9 

The  Elevation  Angles,  Colorado  Plains,  Mountains  and  Ohio 


Median  Values  in  Milliradians 


h  4 

Pi 

h 

g2 

0.  6 

1.7 

3 

6 

9  m 

d,  0 

9 

0 

9 

0 

0  0 

0 

0 

0  0 

km  ei 

e2 

e 

e2 

e 

e2  e 

e2 

e 

e2  e 

Colorado  Plains,  Ah  =  90  m,  N  =  290  ,  a  =  8327  km 

s 


20 

1.  3 

5.  5 

10.4 

6.  5 

10.  3 

6.  1 

10.  1 

5.  1 

9.  3 

4.  3 

8.2 

30 

0.  5 

16.8 

16.  9 

14.  5 

14.  6 

11.9 

12.  1 

8.8 

7.4 

7.  1 

6.9 

50 

0.  9 

2.  5 

5.  7 

2.0 

5.  6 

1.8 

5.0 

1.0 

4.  6 

0.4 

3.4 

80 

-2.  1 

1.9 

0.8 

1.  3 

O 

O 

1 

0.4 

-1.4 

Colorado  Mountains,  Ah  =  650  m,  N  =  290,  a 

s 

=  8327  km 

20 

25.  4 

77.  4 

99.8 

76.  6 

98.6  75.7 

97.  3 

73.  7 

94.  2 

71.  5  91.2 

30 

33.  7 

112.  9 

135.  2 

111.2 

134. 3  109. 2 

133.  3 

104.  7 

129.  4 

100.  1  124.9 

50 

47.  2 

105.4 

15''.  « 

104.  2 

158. 6  102.8 

158.  5 

99.8 

158.  3 

98.  6  158.0 

NE 

Ohio,  Ah  =  90  m,  Ns 

=  312, 

a=  8676  km 

20 

3.  5 

6.  0 

11.2 

5.  6 

11.4 

4.8 

11.0 

4.  0 

9.9 

30 

1.0 

6.  1 

10. 6 

4.8 

9.  6 

3.  6 

9.0 

3.  0 

6.  6 

50 

0.  3 

5.  3 

7.  6 

4.  3 

7.  3 

3.  1 

6.  6 

2.  8 

6.0 
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„surements.  These  values  show  no  consistent  dependence  on  distance. 


and  9  decreases  only  slightly  with  increasing  height  of  the  receiving 
e 

antenna.  Table  2.  10  shows  values  of  9  ,9 

ei  ei  '  e 


,  and  9  calculated  using 
e 

(6a),  compai  ed  with  median  values  of  0  for  each  antenna  height. 


0 


.«i22«  A  Ah_4h  1 

dL„  L  \  LI  /  elJ 


radians  . 


loa) 


2-6  Location  Variability 

The  path-to-path  variation  in  available  wanted  signal  power  is  dis¬ 
cussed  in  annex  1.  Such  random  variations  from  location  to  location 

are  assumed  to  be  normally  distributed  with  a  standard  deviation  <rT  dB. 

La 

An  estimate  of  cr  is  required  to  calculate  the  service  p  bability  Q  . 

l~d  cl 

Analysis  of  path-to-path  variability  of  radio  transmission  loss  for  a 
given  frequency  and  terrain  variance  assumes  statistical  homogeneity  of 
the  terrain.  It  has  been  noted  that  the  plains  and  mountain  areas  show  a 
predictable  change  in  the  variance  of  terrain  from  one  direction  to 
another,  and  that  in  the  area  studied  in  NE  Ohio  the  greatest  terrain  ir¬ 
regularity  occurs  in  the  vicinity  of  the  central  transmitter. 

Transmission  loss  data  from  the  measurement  program  reported  by 
Miles  and  Barsis  (1966)  were  used  to  obtain  an  estimate  of  cr  The 

interdecile  range,  AL,  of  values  of  transmission  loss  recorded  for  each 
frequency,  polarization,  antenna  height  combination, and  distance  was 
tabulated  for  the  Colorado  plains  and  mountain  areas,  and  the  area 
studied  in  NE  Ohio.  These  interdecile  ranges,  given  in  table  2.  11, 
show  no  consistent  dependence  on  antenna  height  combinations  or  on  path 
length,  but  do  increase  quite  consistently  with  frequency  and  terrain  ir¬ 
regularity.  The  interdecile  ranges  of  transmission  loss  AL  were 
plotted  versus  the  parameter  Ah(d)M  and  a  smooth  curve  was  drawn 
through  overlapping  median  values.  The  analytic  expression. 
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Table  2.  10 

The  Elevation  Angles,  Colorado  Plains,  Mountains, and  Ohio 
Median  and  Calculated  Values  in  Milliradians 


96000000000 
ei  e2  e  e2  e  ei  e  ei  e  e  2  e 


Colorado  Plains 


a)  0.9  5.5  10.4  6.5  10.3  6.1  10.1  5.1  7.4  4.3  6.9 

b)  1.5  5.7  7.2  3.0  4.5  2.2  3.7  0.6  2.1  -0.3  1.2 


Colorado  Mountains 


a)  33.7  105.4  135.2  104.2  134.3  102.8  133.3  99.8  129.4  98.6  124.9 

b)  61.7  169- 4  231.1  9  5.  7  157. 4  71.7  133. 441.0  10  2. 7  26.0  87.7 


NE  Ohio 


a)  1.0  6.0  10.6  4.8  9.6  3.6  9.0  3.0  6.6 

b)  1.4  3.5  4.9  1.9  3.3  0.6  2.0  -0.3  1.1 


a)  Median  values  from  proj.'i’"s 

LI,  2 


b)  Calculated  using  (6a), with  d 


calculated  using  (5e) 


0 


Table  Z.  11 

Interdecile  Ranges  of  Transmission 


Frequency  100  MHz _ 

Polarization  Vertical  Horizontal 


hsi  =  4,  hg^  = 

3 

6 

9 

3 

6 

9 

Colorado  Plains 

d  =  10  km.  Z4.  5 

Zl.  9 

Z0.  9 

Z3.Z 

Z5.  5 

Z4.  7 

ZO 

Zl.  Z 

Zl.  8 

Z6  6 

Zl.  8 

Z5.9 

Z6.  1 

30 

3Z .  3 

31.  3 

Z9 .  8 

31.0 

33.  3 

34.  5 

50 

17.  4 

ZZ.  0 

ZZ.  0 

ZO.  0 

Z0.  6 

ZZ.  5 

30 

17.  0 

Zl.  8 

Zl.  0 

19.  7 

18.  4 

17.  8 

Colorado  Mountains 

d  =  1 0  km.  Z5.  4 

31.6 

Z7.  0 

39.  1 

67.  3 

ZO 

Z6.  3 

Z8.  3 

Z9.  Z 

40.  6 

46.9 

43.  4 

30 

30.  5 

Z6.  1 

Z7.  3 

35.  1 

36.  1 

50 

19.6 

Z6.  6 

Z8.  3 

NE  Ohio 

d  =  10  km. 

34.6 

Z6.  7 

Z7.  3 

ZZ.  5 

Z7.  3 

Z7.  5 

ZO 

Z5.  6 

Z3.  4 

Z  3.  5 

Z4.  8 

Z4.6 

Z4.  Z 

30 

33.  5 

Z7.  5 

ZZ.  0 

Z6.  6 

Zl.  1 

Z0.  3 

Z  -40 


AL  in  dB 


Z  3.  4 

ZZ.  1 

9.  4 

ZZ  .6 

16.  0 

9.  0 

17.  5 

15.  8 

19.  6 

16.  3 

17.  5 

8.  6 

Z0.  0  36. 0  38.  3 

Z7.Z  Z7.3  Z9-6 

ZZ.3  Z4.Z  Z  1 . 6 

18.8  1Z.7  16.7 

Zl.Z  Z0.  1  19.9 

16.3  19.4  15.3 

16.6  Z3.1  lb. 7 


(2.2) 


<rLa  =  [0,  i  +  0.  2  K/AM^]"1  dB, 

was  then  fitted  to  these  values,  where  a.  =  0.  AL.  This  function  in- 

La 

creases  rapidly  to  about  9  dB  for  Ah(d)/\  =  20  and  then  slowly  in¬ 
creases  further  to  a  maximum  value  of  10  dB  . 

The  presently  available  data  indicate  larger  values  of  cr  with  hori- 

La 

zontal  than  with  vertical  polarization  at  100  MHz  in  the  mountains,  but 
no  significant  polarization  effect  is  observed  in  the  Colorado  plains  or  in 
Ohio.  Further  studies  of  location  variability  should  be  made,  especially 
at  higher  frequency  ranges  and  with  higher  antennas.  The  estimate  of 
given  by  (2.2)  depends  entirely  on  the  examination  of  data  at  20,  50, 

3. 

and  100  MHz  in  Colorado  and  Ohio. 

2-7  The  Terrain  Roughness  Factor  cr^ 

The  terrain  roughness  factor,  in  (3.  5)  annex  3,  for  line -of- sight 
calculations  represents  the  rms  deviation  of  terrain  and  terrain  clutter 
within  the  limits  of  the  first  Fresnel  zo  .e  in  the  dominant  reflecting 
plane.  For  this  report  the  factor  is  defined  by  (3.6)  as 

cr  ( d)  -  0.  78  Ah(d)  exp  j-  0.5  [Ah(d)]  ^  ** |  m  ,  for  Ah(d)  >4  m  ,(  3.  6a) 

e  (d)  =  0 .  39  Ah(d)  m ,  for  Ah(d)  -  4  m  .  (3.6b) 

h 

These  analytic  expressions  were  developed  from  a  study  of  about  70 
line-of-sight  radio  paths  where  detailed  terrain  profiles  were  available. 
For  each  of  these  paths  the  interdecile  range  of  terrain  heights  Ah(d) 
was  calculated  and  cr  was  computed  using  the  formulas  given  in  section 
5,  volume  1,  and  annex  III,  volume  2,  of  the  report  by  Rice  et  al.  (1967). 
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p 


These  formulas  define  the  points  at  which  the  first  Fr  snel  ellipse  cuts 

the  great  circle  plane.  The  factor  a  was  then  calculated  as  the  rms 

h 

deviation  of  modified  terrain  elevations  relative  to  a  smooth  curve  within 
these  limits 

The  computed  value  of  cr^  was  plotted  versus  the  corresponding 
value  of  Ah(d)  for  each  path.  Equation  (3.6)  defines  a  smooth  curve 
fitted  to  these  computed  values. 
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ANNEX  3 

EQUATIONS  AND  METHODS  FOR  COMPUTING  THE 
REFERENCE  ATTENUATION  A 


cr 

The  minimum  input  parameters  required  to  compute  the  reference 
attenuation  relative  to  free  space  are  the  rudio  frequency  f  in  megahertz, 
the  path  distance  d  in  kilometers,  and  antenna  heights  above  ground  h 

gl 

and  h  in  meters.  Estimates  of  surface  refractivity  N  ,  terrain  ir- 

gl  8 

regularity  Ah,  and  the  ground  constants  o-  and  t  may  be  selected  as 
described  in  section  2,  when  measured  values  are  not  available. 

When  detailed  proiiles  of  individual  paths  are  not  available-,  equa¬ 
tions  (3)  through  (6),  section  2,  are  used  to  estimate  median  values  of 

the  additional  parameters  h  ,  d,  ,  ,  and  0 

r  el  ,2  LI,  2  el  ,2 

When  detailed  profile  information  is  available  for  a  specific  path. 

the  actual  horizon  distances  d  anc  d  ,  horizon  elevation  angles 

L 1  Lj£ 

9  and  8  ,  and  effective  antenna  heights  h  and  h  above  the  dom- 

ei  ei  el  ei 

inant  reflecting  piar  ;  are  used  in  computing  A  .  The  location  of  a 
horizon  obstacle  may  be  determined  by  testing  all  possible  horizon  ele¬ 
vations  and  selecting  the  one  for  which  the  horizon  elevation-  angle  0 


ej 


or  0  is  a  maximum: 
ei 


0'001<hL.-V  dLl 


el 


LI 


za 


radians 


(3.  la) 


0.  001  (h 


L2 


h«  > 

82 


ez 


L2 


_L 2 

2a 


radians  , 


(3.  lb) 


where  h  are  the  antenna  heights  above  sea  level  in  meters,  a  is  the 

81,2 

effective  earth's  radius  in  kilometers,  h  .  are  the  heights  in  meters 
above  sea  level  of  the  horizon  obstacles,  and  d  are  the  great  circle 

i-J  1  ,  £ 

distances  in  kilometers  from  each  antenna  to  its  horizon.  The  picdic- 

tio.n  method  is  limited  to  values  of  0  -  0.  2  radians.  For  larger 

cl  ,  2 
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elevation  angles  the  assumption  of  an  effective  earth's  radius  a,  based 

on  the  surface  refractivity  N  ,  is  not  applicable. 

s 

An  alternative  procedure  is  first  to  compute  a  least-squares  fit  of 

a  straight  line  to  terrain  eleva'ions  above  sea  level.  The  heights  h(  ^  s 

and  h  are  then  defined  relative  to  this  curve  fit,  rather  than  rela- 

LI,  c 

tive  to  sea  level.  This  amounts  to  replacing  sea  level  by  an  arc  of 
radius  "a"  that  is  a  least-squares  fit  to  the  great  circle  path  terrain 
profile . 

For  line-of-sight  paths,  the  effective  antenna  heights  h  are  de- 

6 1  i  i 

fined  as  the  height  of  each  antenna  above  the  dominant  reflecting  plane 
between  the  antennrs,  or  the  structural  height,  whiche%fer  is  greater. 

The  effective  heights  may  be  calculated  as  heights  above  a  smooth  curve 
fitted  to  great  circle  profile  terrain  elevations  that  are  intervisible  to 
both  antennas.  A  straight  line  is  first  fitted  by  least  squares  to  equi- 
distant  heights  h.  ,  and  an  amount  d.  /  2a  is  then  subtracted  at  each  dis- 

l  l 

tance  d.  to  allow  for  the  path  curvature  1/a  .  When  terrain  is  so  ir¬ 
regular  it  ^nnot  be  reasonably  well  approximated  by  one  or  more 
such  reflecting  planes,  the  effective  heights  are  estimated  using  (4a)  or 
(4b)  in  the  main  body  of  this  report. 

The  total  input  required  to  compute  A  is  then:  f,  d,  h  ,  h  , 

cr  gl  gi 

oolarization,  and  actual  or  estimated  values  of  N  ,  Ah,  o  ,  and  t  . 

*  s 

When  available  for  specific  paths,  the  parameters  Ah(d)  ,  d  ,  d  , 

1j  X  JL/4- 


0  ,9  ,  h  ,  and  h  are  also  included  as  input, 

e  i  e  i  e  i  o  z 


3-1.  Two-Ray  Optic6  Formulas  for  Computing  A  and  A 

o  1 

At  distances  d  and  d,  ,  which  are  '  -11  within  radio  line  of  sight, 
o  1 

but  are  so  chosen  that  the  difference  between  the  direct  and  ground- 
reflected  rays  never  exceeds  one  fourth  of  the  wavelength,  the  following 
formula  is  used  to  compute  the  attenuation  relative  to  free  space: 
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p 


f 


A=  -101°glO  H  *  Re  ~  ^  Re  COS(~*T~"  ~  C)j  +Gp‘  101og10(gox  go2)dB- 

(3.2) 

Here  g^  and  g  represent  the  directive  gain  for  each  antenna  in  the 
direction  of  the  other,  v/hile  ZttAt/x.  is  the  path  length  difference  be¬ 
tween  direct  and  ground-reflected  rays,  expressed  in  electrical  radians 
and  in  degrees  as 


- —  =  4.  19  17  x  10  ^  f  h  h  /  d  radians  , 

X  el  ez 

=  2.  4017  x  10  f  h  h^  /  d  degrees  , 


(3.3a) 

(.3b) 


with  f  in  MHz,  h  in  meters,  and  d  in  kilometers.  R  is  the  mag- 

ei  ,  2  e 

nitude  of  an  effective  reflection  coefficient  and  c  is  its  phase  relative  to 
rr  radians.  As suming  matched  polarizations ,  the  median  path  antenna 
gain  may  be  approximated  as 


G  =  10  log,,  (g  g  )  dB  , 
p  6  10  'sol  6oz' 


(3.4) 


and  these  terms  in  (3.2)  then  cancel  each  other. 

No  divergence  factor  is  included  in  the  definition  of  R  since  its  use 

e 

will  not  add  significantlv  to  the  accuracy  of  the  method  described  for  ir¬ 
regular  terrain.  (See  the  smooth-earth  formulas  for  D,  h 
Rice  et  .  (1967). 

Let 

tion  of 


te 


h  in 
re 


^  and  represent  the  directive  antenna  gains  in  the  direc- 

oint  of  ground  reflection.  Then, 

i/ 


R  -  R, 
e  h 


(3.5) 
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Usually,  and  g  =  g ^  ,  unless  beams  are  very  narrow  or  are 

directed  away  from  the  earth's  surface  to  minimize  reflection  from  the 
surface.  ^  is  the  magnitude  of  the  theoretical  plane  earth  reflection 

coefficient,  the  subscripts  h  and  v  referring  to  horizontal  and  vertical 
polarization  respectively,  and  the  factor  ir  ^  in  the  exponent  is  the  rms 
deviation  of  terrain  and  terrain  clutter  within  the  limits  of  the  first 
Fresnel  zone  in  the  dominant  reflecting  plane.  For  this  report  the 
factor  a-  and  the  grazing  angle  4,  are  defined  as  follows: 

1/ 

<r  Ad)  =  0,  78  Ah(d)  exp  |-  0.  5  [Ah(d)]  /  \  m,  for  Ah(d)  >  4m,  (3.6a) 


<j  (d)  =  0.  39  Zih(d),  for  Ah(d)  -4  m, 
h 


(3.6b) 


4,  =  tan"1  [(h  +  hfi2)/  (1000  d)]  . 


(3.  7) 


R  expf-  ( 2 tt  <r  sin  4>)A]  >0.5  and  >  vsin  y  ,  (3.  8a) 

h,  v  L  h  ’ 


Otherwise, 


R  =  R  . 
e  e 


n  rg« g”  xA 

R  =  I -  sin  v 

e  Is  8 
1*01  oz 


(3.8b) 


The  theoretical  plane  earth  reflection  coefficients  R,  ,  R  and  the 

h  v 

phase  angle  c  are  functions  of  the  radio  frequency  f,  grazing  angle  4j  , 
and  the  ground  constants  o-  and  t  .  Their  magnitudes  may  be  read  from 
figures  III.  1  through  HI.  8,  volume  Z  of  the  report  by  Rice  et  al.  (1967), 
or  computed  as  follows: 


x  =  18000  cr /  f  ,  q  =  x/  (Zp) 

1/ 

Z  L  Z  L  Z 

Z  p~  =  [(<  -  cos  v)  +  X  ]  +  (t  -  COS  v) 


(3.9a) 

(3.9b) 
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2  2 

bv  '  2  Z  ’  bh  "  Z  Z 
p  +  q  p  +  q 


radians 


_  2(pt  +  ox)  2  p 

v  Z  Z  ’  mh  -  Z  2  * 

p  +  q  P  +  q 


(3.10) 


(3.11) 


Then 


R  =  fl  +  b  sin  ^  -  m  sin  4>]/  [l  +  b  sin^  y  +  m  sin  (3.  IZa) 

V  V  V  V  yj\/ 

=  [j.  +  b^  s-'  d  -  sin  y]/  [l  +  b^  sin^  y  +  sin  .  (3.  IZb) 

The  phase  angle  c  in  (3.  2)  is  defined  below  for  both  horizontal  and 

ve..i.'  1  polarization,  c,  and  c  .  The  angle  c,  defined  as 
^  h  v  h 


c  =  tan  I — — - - 

h  V  p  +  6in  y 


]  -  tan  |  - - )  radians  (3.13) 

/  \P  -  sm  y/ 


is  always  negative  and  ranges  in  value  from  0  ^  c  ^  -  0.  1  radians.  The 

h 

angle  c  changes  suddenly  from  near  zero  to  tt  /  2  at  the  pseudo-Brewster 

'  - 1  i 

angle,  sin  \l/b  .  To  define  c  ,  let 

v  v 


y^  -  (x  sin  c<  +  q ) /  («  sin  +  p)  ,  y  =  (x  sin  y  -  q)/ (t  sin  y  -  p) 


t  sin  y  >  p  : 


(3.  14) 


-1  -1 

c  =  tan  v  -  tan  v,  +  radians  . 


(3. 15a) 


If  *  sin  -  v  p  and  p  sin  v  0.  5  : 


-1  -1 

c  -  tan  y  +  tan  y  radians  . 

V  1  4- 


It  *  sin  v-  i  p  and  p  sin  v  -  0.  5  : 


(3.  15b) 


-1  -1 

c  =  tan  v,  -  tan  v  radians  . 
v  1  '  Z 


(3.  15c) 


3-5 
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In  the  abo%'e  formulas,  tan  y  is  in  the  first  quadrant  if  y  is  positive 

and  in  the  fourth  quadrant  if  y  is  negative. 

The  two-ray  optics  formulas  (3.2)  to  (3.  15)  are  used  to  compute 

values  of  attenuation  A  and  A  at  distances  d  and  d  ,  respectively. 

ot  l t  o  l  7 

For  A  ,  -  0  ,  define 
ed 

-5 

d  =4x10  h  h  f  km,  or  0.5  d,  ,  whichever  is  smaller,  (3.16a) 
o  ei  e2  L 

For  A  ,  <■  0  ,  define 
ed 

d  =  -  A  /  m  km  ,  or  (d  -  2)  km  ,  whichever  is  smaller  ,  ( j.  lfeb) 

oi  ed  d  L 


d  for  d  2  0.  5  d 
o:  ol  L 


0.  5  d  otherwise 

1  Li 


(3.  16c) 


d  =  d  +  0.  25  (d  -  d  )  km.  (3.  16d) 

i  o  L  o 

In  (3.  16)  the  radio  frequency  i  is  in  MHz,  the  effective  antenna  heights 

h  are  in  meters,  d.  is  the  sum  of  the  horizon  distances  in  kilo- 

e  l,  l  L 

meters  and  the  attenuation  A  ,  and  slope  m  ,  are  defined  in  tiie  next 

ed  d 

subsection  ( 3.  38). 

In  addition  to  the  two-ray -theory  estimates  A  and  A  of  attenua- 

ot  1 1 

tion  at  the  distances  d  and  d  ,  estimates  of  diffraction  attenuation  A  ,  , 

o  l  od 

A  ,  ,  and  A  are  also  computed  at  d  ,  d  ,  and  d  ; 
id  Ls  o  i  Ls 


A  =  A  .  +  m  d 
od  ed  do 

A  =  A  +  m  d 
id  ed  d  i 

A  "A  +  m  d  , 

Ls  ed  d  Ls 

where  A  and  m  are  defined  in  the  next  subsection  bv  (3.  38). 
ed  d 


(3. 17a) 
(3.  17b) 
(3. 17c) 


The  estimates  of  attenuation  A  and  A,  at  the  distances  d  and  d 

o  1  o  1 

are  then  computed  as  weighted  averages  of  the  two-ray  theory  and  the 
diffraction  estimates 


A  =  w  A  +  (1  -  w  )  A  or  A  ,,  whichever  is  smaller,  (3.  18a) 

o  o  ot  o  od  od 

A  =  w  A  +  (1  -  w  )  A  or  A  ,  whichever  is  smaller,  (3.  18b) 

1  o  It  o  Id  Id 


w  =  (1  +  f  Ah  10  )  , 

o 


(3.  18c) 


For  distances  less  than  the  smooth-earth  horizon  distance  d,  ,  the 

L  s 

calculated  reference  value  A  is  defined  by  a  smooth  curve  fitted  to 

c  r 

the  three  values  of  attenuation  below  free  space,  A  ,  A  ,  and  At  ,  at 

o  1  Ls 

the  distances  d  ,  d  .and  d,  . 

o  1  Ls 


For  0  <  d  -  dT  : 

Ls 


A  =  A  +  k,  (d  -  d  )  +  k,  log  fd/d  )  dB. 
cr  o  1  o  .110  o 


(3.  19) 


The  constants  k  and  k?  in  (3.  19)  are  evaluated  as  follows.  First  esti- 

mates  k  ,  k  of  the  slopes  k  ,  k7  in  (3,  19)  are  computed  as 
1  Z  1 


(A  -  A  )  (d  -  d  )  -  (A  -  A  )  (dT  -  d  ) 

.  _  _ Ls  o  1  o _ 1  £  Ls  o 

Z  (d  -  d  )  log  (d  /d  )  -  (d  -  d  )  log  (d  Id  ) 
1  o  1C  Ls  o  Ls  o  10  1  o 

or  0  ,  whichever  is  larger  algebraically, 


( 3 .  Z  0) 


k  =  [(A  -A  )  -  k  log  (d  / d  )]/(d.  -d)  dB/km.  (3.Z1) 

1  *-  I.s  o  Z  10  Ls  o  J  Ls  o 
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If  k^  <  0  set  kj  =  0  and 


k  =  (At  -  A  )  /  log .  (dT  /  d  ). 
2  '  Ls  o'  6 10  '  Ls  o' 


(3.22) 


If  the  reference  attenuation  A  ^  computed  from  (3.  19)  is  less  than 

zero  at  any  distance  0  -  d  S  d  ,  let  A  =0  for  that  distance. 

Ls  cr 


3-2.  Formulas  for  Computing  Diffraction  Attenuation  A, 

d 

In  the  far  diffraction  region,  the  attenuation  A^  is  computed  as  a 

weighted  average  of  two  estimates,  A  for  smooth  terrain  and  A,  for 

r  k 

highly  irregular  terrain.  In  general,  A^  is  defined  by  (13)  as 


A  =  (1  -  w)  A,  +  w  A  dB  , 

Cl  K.  I* 


where  the  empirically  determined  weighting  factor  w  is  defined  as 

I  -1 


w  =  <1  +  0 


L  gl  g« 


a  8  +  d 

e 


(3.2  3) 


with  — —  —  i  1000.  In  the  accompanying  computer  program  and  output 

C  =  0  .  For  low  antennas  with  known  path  parameters  C  -  10  . 

In  (3.  23)  the  radio  w-avelength  \,  te rrain  irregular ity  Ah(d)  ,  and 

effective  and  structural  antenna  heights  h  ,  h  are  in  meters: 

el  ,  2  g  i,  z 

the  effective  earth's  radius  a,  the  horizon  distance  d  ,  and  the  dis- 

L 

tance  d ,  at  which  \  and  A^_  are  computed  are  in  kilometers;  and  the 
sum  of  the  elevation  angles  8  is  in  radians.  For  very  smooth  terrain, 
the  weight  w  =  1  and  A^  =  A  ,  and  for  highly  irregular  terrain,  the 
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p 


weight  w  =  0  and  =  A^  .  The  prediction  approaches  when  either 
the  frequency  or  the  terrain  irregularity  are  very  large;  therefore,  a 
limit  is  placed  on  this  ratio. 

The  diffraction  attenuation  is  computed  at  distances  d.  and  d  , 

*  3  4 

chosen  well  beyond  the  horizon: 

2  t  2  I 

d3  =  d^  +  0.  5  (a  /  f  )T  km  ,  d^  =  d3  +  (a  /  f  )T  km  .  (3.  24) 


If  d  <  d  set  d  =  d  . 

3  Ls  3  Ls 

At  these  distances,  d  and  d  ,  the  attenuations  A_  and  A  are 

3  4  3  4 

computed  using  the  following  formulas,  substituting  d3  and  d^  for  d  in 
(3.23)  to  obtain  and  w  : 


A  =  (1  -  w  )  A  +  w  A 
3  3  kj  3  r  j 


(3.25a) 


0 


3 


A4  =  (1  -  w4>  Ak. 

=  9e  +  V‘*  94 


+  W  A 
3  r  4 


+  d4/a  . 


(3.25b) 

(3.25c) 


The  estimates  A^  and  A  for  highly  irregular  terrain  are  com¬ 
puted  as  though  the  horizon  obstacles  were  sharp  ridges  or  hills,  and  the 
attenuation  is  computed  for  a  double  knife-edge  path. 


v  1 .  3  1Z915  03  Vf  dH  <d3-dL)/(d3-dL2> 


"1.4  *  !'2<m  °4 

"2.4  =  ,'i,1S  °4  V,dL2  <V^IVV 


(3.  26a) 

(3.  26b) 

(3.26c) 

(3.  26d) 
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k.3 


1.3 


2.  3 


rj ,  4 


3,  4 


1. 27 

v  for 

0iv<2.4 

(3.  27a) 

10  v 

for  v  >  2 . 4 

(3.  27b) 

Ak4 

=  A(vx 

,)  +  A(v  )  . 

4  2.4 

(3.27c) 

A 

and  A 

are  defined  as 

r  3 

r< 

<v  - 

F(x2,  - 

20  dB, 

(3.28) 

'h,4> 

depend 

on  the  radio  frequ 

ency. 

polarization,  and  ground  constants  a  and  t,  the  distances  d  ?,  d  , 

and  the  effective  earth's  radii  a  for  the  terrain  between  the  antennas 

I  t  *1- 

and  their  horizons  and  a  for  the  terrain  between  horizons.  The 

3,  4 

latter  are  defined  as 

a  =  d^  /  (0 . 002  h  )  km ,  a  =  df  /(0. 002  h  )  km  ( 3 .  -9a) 

1  LI  ei  2  L2  ez 

a  =  (d  -  d  )  /  9  km,  a  =  (d  -  d  )/0  km.  (  3.  29  b) 

33L3  3  3  i_  4 

Then  the  distances  x  ,  are  defined  as 

1  i  »  ■  »  * 


h  =  Bi  h 


km  , 

X> 

=  B,  a,’  1  dLZ  km 

(3. 

30) 

<V 

"l1 

+  x  +  x  >  km 

(3. 

31a) 

‘V 

dL» 

+  +  x>  km  , 

(3. 

31b) 

where  the  parameter  B  is  defined  for  both  vertical  and  horizontal 

*  i  »  »  ^  i  * 

polarization  as 
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(3.  32) 


l,2,3,4  =  416-“fI  t‘-607  -Kh.v<a1.2,3.4>3 


The  parameters  K.  (a)  for  horizontal  and  K  (a)  for  vertical  polarization 

h  v 

are  defined  as 


I  2  2  1 

K  (a)  =  0.  36278  (a  f  p[(€  -  1)  +  x  ]"u 


(3.  33a) 


K  (a)  =  K  (a)  [tZ  +  xZ]2  , 


(3.  33b) 


where  \  is  defined  by  (3.9a)  as  x  =  18000  u /  f  ,  and  the  ground  constants 
cr  and  <  are  included  in  the  input. 

The  functions  F(x^)  and  F(x^)  may  be  read  from  figures  8.  5  or  8.  6 
of  the  report  by  Rice  et  al.  (1967)  or  may  be  computed  using  the  following 
formulas. 


1.  For  0  <  x,  200  and  0  -  K.  (a,  ,)  5  10' 
1,2  -  h,  v  1,2' 


F(x12)  =  <0  1og10xli2  -  117,  or 


(3.  34a) 


F(x  i  2 )  =  "  117  dB  * 


(3.  34b) 


whichever  yields  the  smaller  absolute  value. 


2.  To v  0  <  x  5  200  and  10  t  K.  (a  )  <  1  , 

i ,  2.  -  h ,  v  1,2 

and  X  >-  450/|lo8l0[KhiV(a12)]|3, 

F(x  5)  is  calculated  using  (3.34a).  Otherwise, 

*  » 


Fh(  2>  *  L°g10  Kh.  vK.  2>  +  ^  5  *  10'  XI,  v^l,  2)  "  15  dD  • 


(3.  34c) 
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Note  that  when  K.  (a ,  )  >  0.  1  no  test  on  x  is  required  and  ( 3.  34c  ,  is 

h,  v  1, 2  '  ' 

always  used. 


3.  For  200  <x  S  2000,  define 
1 ,  2 


w  =  0. 0134  x  exp(-  0.  005  x  ) 
i  #  i  )  ^  1  > 


(3.35a) 


Then 


F(xl,  i)  =  wliJ(«l°g10*1>  2-  "■)+(! “5751  *it2-‘01og10xli  2)dB. 

(3.35b) 


4.  For  x  >  2000  , 
I ,  L 


F<xl,2>  =  05751  *1,2  •  101°SlOXl,2  dB' 


(3.36) 


The  parameter  G(x  )  is  defined  as 

j,  4 


G(x3  4)  =  0.05751  -  10  1og,0x3  4  dB. 


(3.  37) 


Values  of  A  as  given  by  (3.27)  and  of  A  as  given  by  (3.  28) 

K J  ,  4  rj,4 

are  substituted  in  (3.24a,b)  to  obtain  A^  and  A  .  These  computed 

values  cf  A^  at  d^  ,  and  A^  at  d^  are  used  to  compute  the  slope  m^ 

and  intercept  A  ^  that  define  a  straight  line.  The  reference  attenuation 

A  at  any  distance  d,  •  d  -  d  is  then 
cr  Ls  x 


A  =  A  ,  =  A  ,  +  m  ,  d  dB  , 
cr  d  ed  d 


(3. 38a) 


Aed  =  Afo  +  A4  -mdd4’  and  md  =  (A4  -  A3)/{d4  “  d3)  ’  <3'3Sb) 


where  A.  is  a  "clutter  factor",  defined  as 
to 
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A.  =  5  log .  ,  [l  +  h  h  fa  (d  )10_5]dB, 
fo  6 10  L  gl  g2  h  '  Ls'  J 


or  15  dB  ,  whichever  is  smaller, 


(3. 38c) 


and  the  terrain  roughness  term  a  (d  )  is  obtained  by  substituting  d 

n  J—  8  X- •  8 

for  d  in  (3. 6)  and  (3) . 


3-3.  Formulas  for  Computing  Scatter  Attenuation  A 

s 

At  distances  d  and  d  ,  defined  below,  the  following  formulas  are 
5  6 

used  to  obtain  initial  estimates  A  and  A  of  forward  scatter  at- 

5  6 


tenuation  relative  to  free  space: 


d  =  d  +  200  km,  d  =  d  +  400  km  (3.39) 

5  1-  6  L 

0-  =  0  +  d  /a,  9.  =  0  +  d  /  a  radians  (3.40) 

o  e  5  o  e  b 


H5.6-'(r_+'ir_)a95.6f|0'007  -  0-058  \6 

\  ei  ei' 

or  15  dB  ,  whichever  is  smaller. 


dB 


(3.41) 


S  =  H  +  10  log  (f  9.)  -  0.  1  (N  -  301)  exp( -  9C  d _/ 40)  dB  (3.42a) 

5  5  10  T>  s  5  5 

3,  =  H  .  +  10  log  (f  0^)  -0.1  (N  -  301)  exp(-  l  .  /  40)  dB.  (3.42b) 

6  6  iu  6  s  t  6 

Substitute  d  ,  0.  ,  S  ,  and  d  ,  8  ,  S  in  the  following  expressions  to 

6  6  6  6  6  6 

a  a  \  A 

obtain  A  =  A  at  dr,  and  A,  =  A  at  d.  . 

5  s  5  6  s  6 

For  0d  i  10  : 


A  =  S  +  103.4  +  C  .  332  0d  -  10  log,,  (0d)  dB. 
s  10 


(3. 43a) 


For  10  i  Odi  70: 


A  =  S  -t  97.  1  +  0. 212  0d  -  2.  5  log  (0d)  dB. 

8  I U 


(3.43b) 
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For  0d  > 70 


As  =  S  +  86.  8  +  0.  157  6d  +  5  log  (0d)  dB 


(3.43c) 


3-3.  1.  For  H  s  10  dB 
5 

When  the  frequency  gain  function,  H_  ,  computed  at  d  is  less  than 

5  5 

or  equal  to  10  dB  ,  formulas  (3.  39)  through  (3.  4  3)  give  the  actual  pre¬ 
dicted  scatter  loss  at  the  distances  d  and  d  ,  and 

5  6 

A  =  Ar  dB  ,  and  A.  =  A,  dB  . 
o  5  do 

The  scatter  attenuation  A  ,  at  any  distance  d,  is  then  given  by  (17) 
and  (18)  as 


A  =  A  +  m  d  dB  , 
s  e  s  s 


where 


A  =  Ac  -  m  d  ,  and  m  =  (A  -  Ac)/  (d  -  d  )  . 
es  5  s  5  s  6  5  6  5 

The  distance  d  ,  where  diffraction  and  scatter  attenuations  are 
x 


equal,  is 


d  =  (A  -  A  ,)/  (m  -  m  )  km 
x  es  ed  d  s 


(3.44a) 


or  d  t  0.Z5  (a"/  f)!  log  I,  whichever  is  greater  ,  (3.44b) 

Lj  1U 

where  A  .  and  m,  are  defined  bv  (15),  and  (3.38b).  When  (3.44b)  is 
ed  d 

used  to  define  d  ,  redefine  A  as 
x  e  s 


A  =  A  ,  +  ( m  -  m  )  d 
es  ed  d  a  x 


(3.44c) 


The  reference  attenuation  A  for  transhorizon  paths  is  then 

c  r 


for  d  -  d  -  d  ,  A  =  A,  =  A  ,  +  m  d  dB 
Ls  x  cr  d  ed  d 


for  d-d  -  1500  km,  A  =  A  =  A  +  m  d  dB 
x  c  r  s  e  s  s  J 


(3.45) 
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3-3.  2.  For  H  >  10  dB  and  5  15  dB 
5 

When  the  frequency  gain  function  H  computed  at  d  is  greater 

.  .  5  5 

A  A 

than  10  dB,  the  estimates  A  and  A  are  modified  by  comparison  with 

5  6 

the  scatter  loss  expected  over  a  smooth  earth,  Ah  =  0  ,  To  determine 

the  distance  d  ,  where  diffraction  and  scatter  losses  would  be  equal 
xo 

over  a  smooth  earth,  the  diffraction  loss,  with  Ah  =  0  ,  is  also  computed. 

A 

F  or  the  special  case,  Ah  =  0,  let  A  =  A  ,  m,  =m  ,  and  A 

do  ed  do  d  so 

be  the  preliminary  estimate  of  scatter  attenuation  at  d  .  Assume  that 

5 

the  slope  m  is  not  changed.  Then  one  estimate  of  d  is  obtained  by 
s  xo 

substituting  in  (3.44): 


d  =d  =  (A  -  m  dr  -  A  ,  )/(m  ,  -m)  km  . 

xi  xo  so  s  5  do  do  s 


(3. 46a) 


When  H,  is  large,  a  good  estimate  of  d  is 
r  xo 


•\  ?  - 

d  =d  =  d  +0,25  (a  /f)J  logln  f  krn  . 
X2  xo  L  10 


(3.46b) 


B  or  smaller  values  of  H  ,  d  is  the  better  estimate  of  d  ,  and 

5  xi  xo 

for  larger  values  d  is  the  better  estimate.  Therefore,  a  weighted 

function  is  used  to  compute  d  as  follows: 

xo 


d  =d  (3-0.2  H  )  t  d  (0.2  H  -2)  km  (3.46c) 

xo  xi  5  X2  5 


For  Ah  =  0,  scatter  and  diffraction  losses  are  equal  at  d  .  The 

xo 

diffraction  attenuation  A  at  d  is 

xo  xo 


A  =A  +m  d  dB.  (3.-17) 

xo  do  do  xo 


It  is  assumed  that,  in  general,  the  forward  scatter  attenuation  A  at 

sx 

d  =  d  for  any  value  of  Ah  is 
xo 

A  -  A  +  (A  -  A  )  dB  .  (3.48) 

sx  xo  5  'O 
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The  intercept  at  d  =  0  would  then  be 


A  =  A  -  m  d  dB  . 
es  sx  s  xo 


(3.49) 


Substituting  this  value  of  A  in  (3.44a  or  b) 


es 


tance  d  ,  and  for  any  distance  did  , 
x  x 


determines  the  dis- 


A  =A=A  +md  dB  . 
cr  s  es  s 
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f 


3-4.  List  of  Symbols  and  Abbreviations 
In  the  following  list  the  English  alphabet  precedes  the  Greek  alphabet 
and  lower  case  letters  precede  upper  case  letters.  In  general,  upper 
case  letters  are  used  for  quantities  expressed  in  decibels. 


a  an  effective  earth's  radius  that  allows  for  average  refraction 

of  radio  rays  near  the  surface  of  the  earth,  (1)  and  (3.  1). 

a  ,a  effective  earth's  radii  for  the  terrain  between  the  transmit- 
1  £ 

ting  or  receiving  antennas,  respectively,  and  the  corres¬ 
ponding  horizon,  (3.29a). 

a  ,a  effective  earth's  radii  for  the  terrain  between  horizons  at 
3  4 

distances  d^  and  d^  respectively,  (3.29b). 

A  attenuation  relative  to  free  space,  expressed  in  decibels.  At¬ 

tenuation  below  free  space  is  written  as  positive  values  of  A  . 


A  a  predicted  reference  value  of  attenuation  below  free  space, 

cr 


A 

A 


d 

do 


A 

e 


A 

is 


expressed  in  decibels. 

the  diffraction  attenuation  in  dB,  (13),  (16),  and  section  3-2. 
the  diffraction  attenuation  in  dB  equivalent  to  A  ,  but  com¬ 
puted  assuming  a  smooth  spherical  earth,  (3.47). 
attenuation  below  free  space  in  dB  defined  by  (11),  to  sim¬ 
plify  (10). 

estimated  diffraction  attenuation  below  free  space  in  dB, 
extrapolated  to  zero  distance,  (15)  and  (3.  38). 
estimated  scatter  attenuation  below  free  space  in  dB,  ex¬ 
trapolated  to  zero  distance,  (17)  ;•  d  (3.49). 
an  estimate  of  attenuation  due  to  surface  clutter,  (3.  38c). 
an  estimate  of  knife-edge  diffraction  attenuation,  (13). 

A  an  estimate  of  knife-edge  attenuation  computed  at  distances 
k4 

d^  and  d^  ,  respectively,  (3.27). 
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A ^  diffraction  attenuation  in  dB,  computed  at  the  distance  d^  . 

A  diffraction  attenuation  in  dB,  computed  at  the  distance  dT  , 

Ls  L  s 

(3. 17c). 

A  an  estimate  of  the  diffraction  attenuation  over  the  bulge  of 

r 

the  earth,  (13),  and  subsection  3-2, 

A  ,,  A  the  diffraction  attenuation  A  computed  at  distances  d, 
r3  r4  r  3 

and  d^  ,  respectively,  (3.28). 

A  forward  scatter  attenuation  in  dB,  (18). 

s 

A  forward  scatter  attenuation  computed  at  d  =  d  ,  (3.48). 

sx  xo 

A  an  estimate  of  the  forward  scatter  attenuation,  (3.43). 

s 

A  an  estimate  of  the  diffraction  attenuation  over  a  smooth  earth, 


computed  at  d  ,  (3.47). 
xo 

A(v)  an  estimate  of  knife-edge  diffraction  as  a  function  of  the 
parameter  v,  (3.26). 

A  ,  A  attenuation  below  free  space  computed  at  the  distances  d 
o  1  o 


and  d^,  respectively,  (10)  and  (3.  18) 


A  t>  A  ^  estimates  of  attenuation  below  free  space  computed  at  the 
distances  d  and  d^,  respectively,  using  two-ray  optics, 
(3.  18). 

A  ,  A  predicted  diffraction  attenuation  computed  at  distances  d 
3  4  3 

and  d^  ,  respectively,  (3.25). 

A  ,  A  predicted  scatter  attenuation  computed  at  distances  d  and 
5  6  5 


d,  ,  respectively, 
o 


A 

A,,  A  estimates  of  scatter  attenuation  computed  at  distances  d 
5  u  5 

and  d  ,  respectively,  (3.43). 

6 

A 

A  a  preliminary  estimate  of  scatter  attenuation  computed  over 

a  smooth  earth  at  the  distance  d  ,  (3.46). 

5 
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A(0.  5)  a  long-term  me,';an  estimate  of  attenuation  relative  to  free 


space  for  any  particular  set  of  data, 
b  ,  b  parameters  used  in  computing  the  theoretical  plane  earth  .■ re¬ 
flection  coefficients  for  horizontal  and  vertical  polarization, 
respectively,  (3.10). 

B  parameters  us  i  in  c  mputing  the  modified  distances 

1,  2,  3,  4 

*1.2.  3.  4’  ,3-32)' 

c  phase  angle  relative  to  it  radians  of  an  effective  reflection 

c oeffi cient,  (3.2). 

c  ,  c  the  phase  angle  c  for  horizontal  and  vertical  polarization, 
h  v 

re  spec  ely,  (  3.  1  3)  to  (  3 .  1  5) . 
d  great  circle  path  distance  in  kilometers. 

dB  decibels,  10  log  (power  ratio). 

d^  one  of  a  secies  of  equal  distances  at  which  terrain  heights  h. 

are  read,  p.  3-2, 

d^  sum  of  the  distances  d  and  d^  from  each  antenna  to  the 
corresponding  horizon,  (  5d) . 

d  ,  d^  the  distances  from  the  transmitting  and  receiving  antenna, 

respectively,  to  their  corresponding  horizons,  (5c). 

d,  the  sum  of  the  smooth-earth  horizon  distances  d,  ,  and 

Ls  1  s  1 


I .  s2  * 


,  d  distances  from  the  transmitting  and  receiving  antennas, 

1  1.  s  2 

r  spectively,  to  their  corresponding  horizons  over  a  smooth 
earth,  (5a). 

the  distance  at  which  diffraction  and  scatter  attenuations  art- 
equal,  (3.44). 


d  the  distance  at  which  diffraction  and  scatter  attenuation  would 

xo 

be  equal  over  a  smooth  earth,  (3.46c). 
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d  , ,  u  ,  estimates  of  the  distance  d  ,  defined  by  (3.46a)  and 
xl  x2  xo 

(3.46b). 

A 

d  a  preliminary  estimate  of  the  distance  d  ,  (3.46a)  and 

xo  xo 

(3. 46b). 

d^  a  distance  chosen  to  approximate  the  greatest  distance  at 
which  the  attenuation  below  free  space  is  zero  dB,  (3.  16). 


one  estimate  of  the  distance  d  ,  (3.  16), 

o 

a  distance  greater  than  d  but  less  than  dT  , 

o  L 


(3.  1 6d)  - 


defined  by 


distances  defined  by  (3.  24)  at  which  diffraction  attenuation  is 
calculated. 


d 


5 


distances  defined  by  (3.  39)  at  which  scatter  attenuation  is 
calculated. 


e  the  base  for  natural  or  Naperian  logarithms,  e  ^2.7183, 

(3.21) 

f  radio  wave  frequency,  expressed  in  megahertz  (MHz)  in  this 

report. 

F(x  ),  F(x  )  a  function  used  in  computing  diffraction  attenuation, 

1  *— 

(3.  34). 

g  , .  g  directive  gain  of  each  antenna  in  the  direction  of  the  other, 
ol  o2 

(3.2)  and  (3.  5). 

g  , ,  g  directive  gain  of  each  antenna  in  the  direction  of  a  point  of 
r  1  r  2 

ground  reflection,  (3.  5). 

Gp  path  antenna  gain  expressed  in  decibels  above  the  unit  gain  of 
an  isotropic  radiator,  (3.  3). 

G(x  ,  x  )  a  function  used  in  computing  the  diffraction  attenuation  at 

the  distances  d  and  d  ,  (3.  37). 

3  4 
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h  subscript  referring  to  horizontal  polarization, 

h.  any  one  of  a  series  of  equidistant  heights  of  terrain  above 

sea  level. 


h  a  height  in  meters  used  in  computing  the  horizon  distances, 

e 

(5c). 

h  ,  h  effective  antenna  heights  of  the  transmitting  and  receiving 
el  e2 

antennas,  respectively,  (4)  and  (3.  3). 

h  ,  h  structural  antenna  heights  above  ground,  (4). 

8^8^ 

h  ,  h  height  above  sea  level  of  the  horizon  obstacle  for  the  trans- 
mitter  and  receiver,  respectively,  (3.  1). 

h  height  of  the  surface  of  the  ground  above  sea  level,  (2). 

s 

h  ,  h  height  above  sea  level  of  the  transmitting  and  receiving  an¬ 
si  s2 

tennas,  respectively,  (3.1). 

H^,  frequency  gain  function  computed  at  the  distances  d^  and  d^  , 
respectively,  (3.41). 

k  a  coefficient  used  in  defining  effective  antenna  heights,  (4). 

k  ,  k  coefficients  that  define  the  slope  of  a  smooth  curve  of  A 
12  cr 

versus  distance  for  distances  0  i  d  ^  d,  ,  (10)  and  (3.  19). 

Ls 

k  estimates  of  the  coefficients  k^  and  k^  ,  (3.20). 

K  (a),K  (a)  parameters  for  horizontal  and  vertical  polarization. 


log 


10 


L. 


bf 


L 


cr 


m 


respectively,  used  in  computing  diffraction  attenuation,  (3.  3  3). 
logarithm  to  the  base  10  . 

basic  transmission  ’oas  in  free  soace,  (9). 
median  reference  value  of  transmission  loss  (8). 

ope  of  the  curve  of  diffraction  attenuation  A  versus  dis- 

d 

tance,  (loM  ~nd(3.38). 
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m,  slope  of  diffraction  attenuation  in  dB/ km  for  Ah  =  0,  (3.46a). 
do 

m  slope  of  the  curve  of  scatter  attenuation  A  versus  distance, 

s  ,  s 

(17b)  and  (3.45). 

m^,  parameters  used  in  computing  the  magnitudes  of  the  theoret¬ 
ical  plane  earth  reflection  coefficients  R,  and  R  ,  (3.11). 

h  v 

N  the  surface  r  ef.ractivity,  (2). 

surface  refractivity  reduced  to  sea  level,  fig.  1. 
p  parameter  used  in  computing  the  theoretical  plane  earth  re¬ 

flection  coefficient,  (3.9b). 

q  pa  rameter  used  in  computing  the  theoretical  plane  earth  re¬ 

flection  coefficient,  (3.9a). 

R  the  magnitude  of  an  "effective"  reflection  coefficient,  (3.  5) 

and  (3.8). 

R^,  R^  the  magnitude  of  the  "theoretical"  plane  earth  reflection  coef¬ 
ficient  for  horizontal  and  vertical  polarization,  respectively, 
(3.  12a)  and  (3.  12b) . 

estimate  of  an  effective  reflection  coefficient,  (3.  5). 

S  ,  S,  terms  defined  by  (3.42)  that  are  used  in  estimating  the  for- 
b  o 

ward  scatter  attenuation  A  ,  (3.43). 

s 

v  subscript  referring  to  vertical  polarization. 

v  ,  v^  ,  parameters  used  to  compute  the  double  knife-edge  at- 

v  ,  v  tenuation,  (3.25)  thr  ough  ( 3. 27) . 

1.42.4 

w  weighting  factor,  determined  empirscally  as  a  function  of 

radio  frequency  and  terrain  parameters,  (3.23). 
w  parameters  used  in  computing  rounded  earth  attenuations, 

C 

defined  by  (  3.  3  5a) . 
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Wy  estimates  of  w  corresponding  to  d  =  d^  .  d  =  d^  ,  (3.25). 

x  a  parameter  used  in  computing  the  theoretical  plane  earth 

reflection  c oefficient,  (3.9a). 

2  3  4  distances  defined  by  (3.  30)  through  (3.  31b). 
y i »  y^  parameters  used  in  computing  the  theoretical  plane  earth 
reflection  coefficient,  (3.  14). 

Ah  an  asymptotic  value  of  Ah(d)  which  is  used  to  characterize 
the  terrain,  table  1,  and  (3). 

Ah(d)  interdecile  range  of  terrain  heights  above  and  below  a  straight 
line  fitted  to  elevations  above  sea  level,  (3). 

Ar  the  difference  in  path  length  of  the  direct  and  reflected  ray, 

(3.  2)  and  (3.  3). 

i  the  permittivity  or  relative  dielectric  constant  of  the  ground, 

(3.9)  through  (3.  15). 

0  angular  distance  for  a  transhorizon  path,  (7). 

0^  the  sum  of  tv  '  elevation  angles  0  ^  and  0  ,  (6b). 

0,0  the  angles  by  which  the  horizon  rays  are  elevated  or  de- 
el  e  2 

pressed  relative  to  the  horizontal  at  each  antenna,  (6a)  and 
fig.  2. 

6  ,  0  angular  distances  corresponding  to  d_,  d«'  d - »  <1  in  (3.  26), 

3  ‘i  3  4  5  6 

0C.  0.  (3.29),  (3.40),  (3.  41)  and  (3.  42). 

o  0 

\  radio  wave  length,  used  for  example  in  (3.  23). 

c  the  conductivity  of  the  earth's  surface,  (3.9)  and  following. 

the  rms  deviation  of  terrain  and  terrain  clutter  within  the 
h 

limits  of  the  first  Fresnel  zone  in  the  dominant  reflecting 
plane,  (3.6). 

4j  the  grazing  angle  of  a  ray  reflected  from  a  point  on  the  sur¬ 

face  of  a  smooth  earth,  (3.  7). 
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3-5.  Computer  Program  Listing  and  Sample  Output 
A  computer  program  listing  and  a  sample  of  the  output  are  given  in 
this  section.  The  program  is  written  in  Fortran  IV  for  a  dibital  com¬ 
puter.*  A  list  that  relates  program  symbols  to  corresponding  term9 
and  equations  in  the  report  is  provided  as  well  as  a  brief  flow  chart  of 
the  program.  Attenuation  is  computed  at  fixed  distances  in  addition  to 

the  parameters  required  to  obtain  curves  of  A  versus  d. 

cr 

The  sample  output  shows  computations  for  paths  in  northeastern 

Ohio  and  in  the  Colorado  plains  and  mountains  that  correspond  to  those 

where  measurements  were  made.  The  curves  of  A  versus  d  shown 

c  r 

in  figures  3  through  8  were  plotted  from  this  output.  Calculations  were 
made  at  frequencies  of  100,  50  and  2u  MHz,  for  appropriate  antenna 
heights  above  ground.  The  terrain  of  the  area  Ohio  and  the  Colorado 
plains  is  characterized  by  Ah  =  90  while  for  the  mountain  paths 

Ah  =  65°  m.  Values  of  surface  refractivity  used  are  N  =  312  in  Ohio 

s 

and  =  290  in  Color.,  io.  For  the  longer  mountain  paths  a  somewhat 
lower  value  of  N  would  be  appropriate.  The  terrain  parameters  d  , 

S  Li  1 

u  ,  9  and  9  were  calculated  using  equations  (5)  and  (6),  and  for 

Lc  ei  a 

these  low,  randomly  located  antennas  we  assumed  h  =  h  .In 

ei.z  gi.z 

each  area  the  first  three  sets  at  100  MHz  are  for  vertical  polarization 
and  the  second  three  sets  are  for  horizontal  polarization.  At  frequen¬ 
cies  of  50  and  20  MHz  only  vertical  polarization  is  shown.  The  smooth- 


earth  horizon  distances  d,  ,  and  corresponding  attenuation  Ar  are 

Ls  Ls 

listed.  Similar ly.  the  distance  d  ,at  which  diffraction  and  scatter  at- 

x 

tenuation  are  equal,  and  the  corresponding  attenuation  A  are  given. 

dx 


*The  program  was  w-ritten  for  a  Control  Data  Corporation  CDC-3600 
computer  and  may  require  slight  modification  for  use  with  other  com¬ 
puters. 
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From  this  output  the  calculated  reference  attenuation  A  ^  may  be 


obtained  at  any  desired  distance  d 


F  or 

0  <  d  <  d. 

Ls 

F  or 

dT  <  d  s  d 
Ls 

F  or 

d  >d 

X 

A  =  A  +  k  d  +  k  log  d  dB. 
cr  e  1  2  10 


A 

=  A  J 

+ 

m  ,  d 

dB 

cr 

ed 

d 

A 

=  A 

+ 

m  d 

dB 

cr 

es 

8 
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Reference  List  of  Program  Symbols 


Program 

Report 

Equation 

Program 

Report 

Equation 

A 

a 

U) 

AOG 

A 

o 

ACR 

A 

cr 

(8) 

A  1 

A1 

(3.  18b) 

AD 

Ad 

(13) 

Al,  2 

ai,  2 

(3.29a) 

ADO 

Ado 

A3,  4 

*3,  4 

(3.29b) 

ADX 

Adx 

A3,  4 

A  3,  4 

(3.25) 

AED 

Aed 

(3.  38b) 

B 

b 

(3. 10) 

AES 

A 

e  s 

(18) 

Bl,  2 

Bl,2 

(3.  32) 

AFO 

Afo 

(J.  38c) 

B3,  4 

B3,4 

(3.  32) 

AG 

A 

(3.2) 

C 

c 

(3.  13)  to  (  3.  15) 

A  H5,  6 

A 

5,6 

(3.48) 

D 

d 

AH  50 

>> 

o 

(3.48) 

DEDO 

Aod 

(3. 17a) 

AK3,  4 

Akj,  4 

(3.27c) 

DED1 

Aid 

(3.  17b) 

ALS 

ALs 

(3. 17c) 

DH 

Ah 

(3) 

AR 3,  4 

A 

r  j,  * 

(3.28) 

DHD 

Ah(d) 

(3) 

AS 

A 

s 

(17) 

DHD3,  4 

'ih(d3,  4) 

(3) 

ASX 

A 

sx 

(3.48) 

DL 

dL 

(5d) 

A V  1  3,  23 

A(V1.3. 

2>3)(3.27) 

DLS 

dLs 

(5b) 

A V 14,  24 

Alvl.«. 

2  4)(3.27) 

DL  1,  2 

dLl,  2 

(5c) 

AXO 

A 

xo 

(3.47) 

DLS  1 ,  2 

dLs  1 ,  ^ 

(5a) 

AO 

A 

0 

(3.  18a) 

DX 

X 

(3.44) 
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Reference  List  of  Program  Symbols  (continued) 


Program 

Report 

Equation 

Program 

Report 

Equation 

DXO 

d 

xo 

(3.46) 

MS 

m 

s 

(18) 

DX 1,  2 

d 

x  1,  2 

(3.46) 

NS 

N 

s 

(2) 

DO 

d 

o 

(3.  16c) 

P 

P 

(3.9b) 

DOl,  2 

do,.2 

(3.  16b) 

POL 

polarization 

D1 

dl 

(3.  I6d) 

PSI 

4* 

(3.7) 

D3,  4 

d3,4 

(3.24) 

Q 

q 

(3.9a) 

D5,  6 

i  , 

j,  6 

(3.  39) 

RE 

R 

e 

(3.  5)  &  (3.8) 

E 

c 

(3.  9)  to ( 3. 

15)  S 

cr 

(3.9) 

F 

f  in  MHz 

SH 

"h 

(3.  6) 

FX1,  2 

F(*i. 2 

(3.  34) 

SHDLS 

(3.38) 

GX3,  4 

(3.37) 

SP 

sin 

(3.8a) 

HIE,  2E 

h 

ei  ,  z 

(4b) 

S5,  6 

S5,6 

(3.42) 

H1G, 2G 

h 

gi  ■  * 

figure  2 

TD 

ad 

(3.43) 

H5,  6 

h5.  6 

(3.41) 

TE 

0e 

(6b) 

Kl.  2 

S.2 

(3.  10)  to (3 

.22)  TE  1,  2 

0 

ei ,  z 

(6a) 

Kl,  2 

(3.  33) 

T3,  4 

93,4 

(3.  2  5c) 

K  3,  4 

K<i3,  4> 

(3.  33) 

T5,  6 

9  5,  6 

(3.40) 

M 

m 

(3.  11) 

V 1 3,  23 

V, 

1.3,  2.  3 

(3.26) 

MD 

md 

(3.  38b) 

V 14,  24 

V 

1.4,  2.4 

(3.  26) 

MDO 

nido 

(3.47) 

W 

W 

o 

(3. 18) 
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Reference  List  of  Program  Symbols  (continued) 


Program 

Report 

Equation. 

Wl,  2 

Wl,2 

(3.  35) 

W  3,  4 

W3,4 

(3.23) 

X 

X 

(3.9a) 

XI,  2 

*1.2 

(3.30) 

X3,  4 

X3.4 

(3.31) 

Yl,  2 

yl.2 

(3.14) 
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Line-oi’- Sight  Subroutine,  LOS 


KM  =  2 

Compute  | 

Call  DIFF 
KM  =  0 
Compute 

DO,  Dl,  DEDO 

'1  DEDl.DHD  . 

A  ED,  MD 

(entr y  2) 

Compute 


X,  Q,  P,  B,  M  for 
RH  or  RV 
SH  SP 
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Scatter  Subroutine,  SCATT 
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Computer  Program  I  1 3ting 


C 

C 


C 

c 

c 


5  6 


4  7 

58 

59 


PROGRAM  COMTE 

PROGRAM  TO  DETERMINE  PARAMETERS  AND  WRITE  OUTPUT 
COMMON  /M/F.D.NS.A.DH.DHS.S.E.POL.KM 

COMMON  /HP/  HIE .H2E.H1G.H2G.DLS1 .DLS2.D11 .DL2.DL.DLS.TE1.TE2.TE.KL 
COMMON  /MLDS/  AG.AD.AS.ACR.AED.MD.AH50.AH5.D5.MS.AES.DX.H5 
COMMON  /ml/  P0.D1.D01.D02 .AO.AI.k: .K2.AL.ALS.A0G 
DIMENSION  ANSI  3 ) ,DKM( 6 ) . D E LH I  6 ) » SD ( 6 ) . SA ( 6 ! 

REAL  NS.MD.MDO.MS.MSS.MDS .K1.K2.K3 .K4.LBF 

DATA  <ANS*290. .290.  ,312.  ) 

data  (DKM»5..10.,20..30..50..80.) 

DATA  (DELH-105.  .165.  .23A.  ,315.  .575  .  ) 

CALCULATION  OF  INPUt  PARAM:  £Rs 

S-.005  S  E*15. 

DO  500  I  X- 1 . 3 
NS*ANS I  I  X l 

A*6370./(1.-.04665»£XPF(.  00  5577*NS  )  ) 

WRITE  i  2.56)1 

'ormat  < r i ) 

IP  (IX  .EQ.  1 )  WRITE  (  2.57) 

IF  ( IX  .EQ.  2)  WRITE  (  2.58) 

IF  (IX  .EQ.  3)  WRITE  (  2.59) 

FORMAT  ( 2X .*C0L0PAD0  PLAINS  NS-290.*//) 

FORMAT  ( 2X .‘COLORADO  MOUNTAINS  NS=29p.*//) 

FORMAT  ( 2 X . *0h I  0  NS*  312.*//) 

DO  400  1*1.9 
KK  =  0 

DO  300  12*1.6 

IF  (IX  .EQ.  2  .AND.  12  .EQ.  61  GO  TO  300 

IF  (IX  .EQ.  3  .AND.  (12  .EC.  1  .OR.  12  .EQ.  6))  GO  TO  300 

D*DK“ (12) 

DH»pHS»90. 

IF  (lx  .EQ.  2!  DH*PH$*650. 

F  -  100. 


IF 

(I  .GT. 

6  ) 

F  *50  . 

IF 

(I  .EQ. 

9  ) 

F  *20  . 

POL 

*♦  1  • 

IF 

(I  .GT . 

3 

•  AND.  1 

[  ,LT.  7 

) 

POL •-  1  . 

1 1  C»  H  ;  E  *  4  . 

I  F 

(IX  . NE . 

3 

.  ANO. 

I  .EQ. 

9  ) 

H  1  G»  H  1  E  * 

1 

.  3 

IF 

(IX  .FQ. 

3 

.  AND  • 

I  .0T. 

6  > 

H1C.*H1E* 

4 

.24 

I  c 

(IX  .EQ. 

3 

.AND. 

;  .eg. 

9  ) 

h  1 G  *  H I  E  » 

** 

.6  8 

H2u*h2E*3. 

IF 

< 1  .EQ. 

2 

.OR.  I 

.EQ.  5) 

H 

2G*h2  E  -6  • 

;  f 

( I  .EQ. 

3 

.OR.  I 

•EQ.  6) 

M2G*HJ E  »  9 . 

IF 

(IX  . EO. 

3 

•  AND. 

I  .  EQ  • 

7  ) 

H2G-H2E* 

1 

• 

IF 

(IX  .NF . 

3 

.  A  N  P  . 

I  .EQ. 

7  ) 

H2G*»-'2E* 

• 

•s  S 

I  F 

(IX  .NE . 

3 

.AND. 

I  .EQ. 

8  ) 

H2G*  t  * 

1 

.  7 

IF 

(IX  . Nr . 

j 

•  AND. 

I  .EO. 

9  ) 

h  2  G  *  h  2  E  ■ 

1 

.  3 

DLS1-SQRTF ( .00 2* A* HIE ) 
DLS?»SQRTF ( ,OL2*A*h2E ) 
DlS*OLS1 »PLS2 
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Computer  Program  Listing  (continued) 


DL1=DLS1»EXPF( -  •  07  *SQR  T  F ( DH/HlE )  ) 

0L2*DLS2*EXPF (-.07*S0RTF( DH/H2E) ) 

DL«DL1+DL2 

TE1»(  .00065  /DLS1  )*  t  (  DLS1/DL  1-1.  )  *2  •‘-3 -077 *H1  E  ) 
T E 2  =  {  .0006  5/ 0LS2  )*  < (DLS2/l_: - i. ) *DH-3.077«H2E ) 
T  £xMAXlf (  ( TE1+TE2)  .(-DL/A)  ) 

XX*XX-Fl 

IF  (D  .GT.  OLS)  GO  TO  40 
CALL  LOS 
50 l XX  )  «D 
SA(XX)*ACR 
GO  TO  300 
40  CALL  DIFF 
SD  (  KtC. )  »0 

sa  un  =acr 


IF 

(  I  X 

.EQ. 

1  .AND. 

XX 

.NE. 

6  ) 

GO 

TO 

300 

I  r 

1  IX 

.EQ. 

2  .AND. 

XX 

•  NE. 

5  ) 

GO 

TO 

300 

IF 

(  IX 

.EG. 

3  .AND. 

XX 

•  NE . 

4  1 

GO 

TO 

300 

AE  : 

:  AOG- 

•X  1  *00 

-X  2*  AlOG 1 0 I 

00  ) 

c 

C  WRITE  OUTPUT 

WRITE  (  2.60!  F ,DH .H1G.H2G. Tf ,DX 

60  FORMAT  ( 4X .*Fs*F6. 1 .*  DH«»F6.2»*  H1G=»F6.2*»  H2G*»F6.2» 

C*  TE*»F10.6.»  0X»»F8.2) 

WRITE  (  2.61)  AE.X1.X2.0LS.ALS 

61  FORMAT  ( 4X ,»AE**F8.2 .*  X1**F10.5.*  <2=»F10.5.*  DL5'*F8.2. 

C  *  A|_S*»F8.2) 

aox>aed.md»ox 

WRITE  (  2.62)  AED.M0.AES.MS.40X 

62  FORMAT  { AX » *AED»*F  8 • 2  *  *  MP.»F10.5»*  AES»*F8.2.*  MS»*Fie.5. 

C*  A0X.WF8.2t 

WRITE  (  2.63)  C  S  0 , J  2 ) .  J2*1.XX) 

WRITE  l  2.64)  ( SA( J2 )  .  JZ»1  .XX ) 

63  FORMAT  ( 4X .»D»6F  10 .2  ) 

64  FORMAT  (  wx  .  »  A  »6F  10 .2  ) 

wR  I  rE  (  2.56  1 
300  CONTINUE 
400  CONTINUE 
600  CONTINUE 

call  Exit 

END 


SUBROUTINE  01 FP 

C  SUBROUTINE  TO  COMPUTE  DIFFRACTION  ATTENUATION 

C 

COMMON  /NR/  JZ.W.SW3  .SW4,S *3.5*4. SAFO 
COMMON  /MAR  14  /03  .04  .T 5 
COMMON  /M/F .D.NS.A .OH ,0HS.S.E, POl .XM 

COMMON  /MP /  HIE .M2E .M1G.H2G.PL  SI .CLS2 .OL 1 .DL2 .DL .OLS.TE 1 . TE2 . TE .XL 
COMMON  /MIPS/  AG.AD. AS. ACR. AEO.MO. AH60 . AH6 ,Q5 , MS. AES.OX .MS 
real  NS.MO .MPC.MS.MSS.MDS  -  XI  .X2  .X3  .X4 
FNA(C ) -6.02+9. 11*0-1 .2’*C*C 
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Computer  Program  Listing  (continued) 


FNB(C)«=12.93  3  +  20.*ALOG10(C> 

FNC(C)=416.4*F**. 3333333333*1 1 .607-0 

FNO (C )  =  (. 362 78/! C* F )**. 3333333333  > *  1 . / ( < E-l . ) **2  +  X*X ) ** . 25 
FNE ( C ) =C*SOR  TF(E*E+x*X) 

rdl-dl 

KK-0 

10  ICIt.ltKM 

D3«DL+. 5* ( A*A/F )**. 333333 3333 
IF  ( D3  .LT.  DLS )  D3-DLS 
D4-D3+1 A*A/F ) **.3333333333 
T3-TE+D3/A 
T4-TE+D4/A 
C 

c  CALCULATION  of  knife  edge  diffraction 

c 

V13«1.2915*T3*SORTF| F*DL 1  * ( D3-DL ) / (D3-DL2  )) 
V23*1.291F*T3*SQRTF(F*DL2*tD3-DL)/(D3-DLl)) 
V14-1,2915*T4*SQRTF(F*DL1*( D4-DL ) / (D4-DL2  )  ) 

V24«l .2915*T4*S0RTF ( F«DL2*( D4-DL ) / ( D4-DL1 ) ) 

AV13-FNAI V13  ) 

I F l V 1 3  .GT.  2.4)  AV13-FNB t V13  ) 

AV23-FNAI V23  ) 

IF(V23.GT.2.4  )  AV23*FNB( V23  ) 

AV 1 4  ®  FN  A ( V 1 4  ) 

IF  ( V 1 4  .GT.  2.4)  AV14*FNB(V14) 

A V 2 4 'FNA (  V24  > 

IF (V24.GT.2.4  !  AV24  =  FNB1V24) 

AK3» AV13+AV23 
AK4. AV14+AV24 


CALCULATION  OF  ROUNDED  EARTH  DIFFRACTION 

A  1  «D L 1 *DL 1 / ( . C02 *H  1  E  ) 

A2«DL2*DL2/( .002*H2E) 

A3- ( D3-DL !  /T  3 
A4 - ( D4-DL ) /T4 
X. 18000. *S/F 
K 1 -FNO (All 
K2-FND  f A? ) 

K3-FND!  A?) 

K4  -F.ND  !  A4  1 

IF  (POL  .EQ.  -l.i  GO  TO  15 
A.1-FNE  (Kl) 

K2-FNE < K2  ) 

K3-FNE ( K  3 1 
K4-FNE ( K  4 ) 

15  B 1 -FNC ( K  1  ) 

B2-FNC (  K?  ) 

83*FNC ( K  3  ) 

B4-FNC ( K4 ) 

X1-B1*DL1/A1*«. 66 66666666 
X2-B2*DL?/A2**.656646fc6E66 
X3*B3*(  D3-DL  )  / A3 ••  .fcfr6666E  ‘-♦X1+X2 
X4-B4*(D4-DL)/A*.**.666fc6666fr*Xl*X2 
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Computer  Program  Listing  (continued) 


p 


XL  1 *450 • /ABSF (ALOG101X1 )**3) 

XL2*450. /ABSF (ALOG 10 <X2)**3) 

IF1X1.GT.0..AN0.X1.LE.  200.  .AND. XI  .GE.  0..AND.  XI. LE.  .t  301) 

C16 . 17 

16  T«=40.*ALOG10(X1)-117, 

T1--117. 

T 2  =  M I N 1 F ( ( ABSF ( T ) > . ( ABSF ( T1  ))  ) 

FX  1  *T 

IF  ( T 2  .EQ.  ABSF  I  T  1  )  )  FX1-T1 

17  IF (X2.GT.0..AND.  X2.LE.200.  .AND.  A2  .GE.C..AND.  X2.LE.  .00001) 

C  18  *  19 

18  T*40.*ALOG10( X2 1-117. 

Tl*— 117. 

T2-M I N 1 F ( ( ABSF ( T ) ) . ( ABSF ( T 1  )  )  ) 

FX2-T 

IF  ( T 2  .EQ.  A8SFIT1  )  )  FX2=T1 

19  I F ( X 1  . GT . 0 .  .AND. XI  .LE.200.  .AND. XI  .GT.  .00001  .AND.  XI  .LT.  1. 
C)  21.22 

21  F  X 1 *4  0 . *ALOG 10(X1)— 117. 

IF  (XI  .LE.  X  L  1  )  FX  1  *20.  *  A  LOG  10  (XI  >-*-2.5*1.  E-5*X  1  *X  1  ft  1- 1  5  . 

22  IFIX2  . GT . 0 .  .AND.X2  .LE.200.  .AND.X2  .GT.  .00001  .AND.  XI  .LT.  1. 
C)  23.24 

23  FX2*40.»ALOG10(X2)-117. 

F  (  X  2  .LE.  X  L  2  )  F  X  2  *  2  0  .  *  *  .  1  '  X  2  )  •*•  2 . 5  *  1  .  E-t,*X2*X2/X2-15. 

24  W1«.0134*X1»EXPF(-.005*X1 )  $  W2 « .0 1 34 »X 2 *EX PF ( - . 00 5 »X 2  ) 
IF(X1.GT.200.  .AND.  XI. LE. 2000.) 

C  FX1  =  W1*(43.*ALOG10(X1 : - 1 1 7 . )  -M  1 . -Wl ) * ( .05 75 1 *X 1 - 1 0. * ALCG 1 0 ( X 1  )  ) 

I c I X2  .GT.  200.  .AND.  X 2  .LE.  2000.) 

CFX2*W2*(40.*ALOG10(X2)-117. )+( l.-W2>*l . 0 5 7 5 1 *X 2- 1 0 . * ALOG 10 t X 2 ) > 

I F ( X 1 .GT .  2000.)  FX1».05751*X1-10.*AL0G10(X1) 

I F ( X  2  .GT.  2000.)  FX2*.05751*X2-1p. *ALOG 1 0 1  X  2  ) 

GX3*. 05751 *X  3- 10 . • ALOG 1 0 ( X  3 ) 

GX4».05751*X4-10.* ALOG 101X4) 

A93sGX3-Fxl-FX2-20. 

AP4.GX-.-FX1-FX2-20. 

r 

C  COMBINATION  OF  ROUNDED  EARTH  ANC  XNIFE  EDGE  DIFFRACTION 

C 

28  DHD3-DH* 1 1 .-.8  »EXPF ( - . 0 2 • D3 ) ) 

DHD4»DW* ( 1 ,-.8  *EXPF ( -,C2*D4 ) ) 

P13*SGRTF( (H1E»H2E)/(H1G*H2G) ) ♦ ( A» TE +DL ) /D3 
p I  4  *  SQR  T  F 1  (H1E»H2E ) / tHlG»H2G>  ) ♦ ( A* Tt +  DL ) / D4 
DOL3*M!NlF(1000..(DHD3*F/299.7925)) 
DuL4«MlNlF(I00C.»(DHC4»F/299.7925i) 

63  W3*l./(l.*.l  •SQRTF(DOL3»P13) ) 

W4  *  1  .  /  {  1  .  ♦  .  1  *SQRTF(o'CL4»P14)  ) 

31  A3«( 1.-W3) *AX3>W3»AR3 

A  4  * i  1.-W4) »AX4  +  W4*AR4 
MD*( A4-A3) /(D4-D3I 
AED  *A4-MD«  D4 

DhDl  S*Dh» (  l.-.8«ExPF(-.02*DLS>) 

SHDLS*. 78*DHDLS*ExPF (-.5»(DHDLS»».25!) 

AFO»5.MLOG10(  l  . +H1G*H2G»,:»SmDLS*.  00001  ) 

AFO-M IN1F 1 Ar 0 . 1 5  .  ) 
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imputer  Program  Listing  (continued) 


91  AED* AED+AF  0 

IF  (KM  .EQ.  2  )  GO  TO  40 

IF  (KK.EO.  2)  GO  TO  20 
SAFO-AFO 

SW3*W3  $  SW4*W4  5  SA3-A3  $  SA4**4 

29  AD= AED+M0*D 

32  TO- ( TE+D/A)*D 

SDL1-DL1  5  SDL2-DL2  $  SDL*DL  $  STE1*TE1  S  STE2*TE2 
STE-TE 

CALCULATION  OF  SCATTER  ATTENUATION 
CALL  SCATT 

IF  ( H5  .LE.  10.)  AES-AH5  -MS*D5 
IF  <H5  .LE.  10.1  AS*AES+MS*D 
IFIH5.LE.10. )  GO  TO  30 
MDS-MD 
AEDS*AED 
dh*o» 

GO  TO  10 
2  C  ADC*AE0 
MDO*MD 
KD*MPS 
AED-AEDS 

0X1*1 AH50-MS*D5-AD0 ) / ( MDO-MS ) 

DX2*1RDL+.25*(A*A/F !*«.3333333333«ALOG10(F) ) 

DxO=DX 1* ( 3.-. 2*H5 )+OX2»(.2*H5-2.) 

AxO*AOO+M 00*0X0 
ASX  *  A  X0+ 1 AH5-AH5C) 

AES* ASX-MS*PXO 
AS*AE  S  +  MS*D 

30  0X*(AES-AE0) /(MO-MS) 

DxN«Dl+.25*( A*A/F ) **.3333 3333 33*AIOG10 (F ) 

IF  ( DXN  ,GT.  OX)  AES-AEO+ I Mp-MS ) *DX 
IF  (DXN  .GT.  DX  )  OX  *  DXN 
ACR* AD 

IF  (D  .GT.  DX)  ACR  *  A  S 

DL  1  *  SDL  1  *  DL2-SDL2  $  DL  *  SD  L  i  TE1-5TE1  %  1E2-STE2 

TE*STE 

DH»DHS 

41  CONTINUE 

40  RETURN 

End 


SUBROUTINE  SCATT 

C  SUBROUTINE  TO  COMPUTE  SCATTER  PARAMETERS 

C 

COMMON  /M/F.D.NS.A.DH.DHS.S.E .POL . KM 

COMMON  /MP/  HIE.H2E.H1G.H2G.DLS1.DLS2.DL1.DL2.DL.DLS.TE1.TE2.TE.KL 
COMMON  /MAR14/D3  .04  .  T5 

COMMON  /MIPS/  AG . AD. AS . AC R . AED.MD. AH50 . AH5 .05 .MS , AES .DX .H5 
REAL  NS. MD .MpO.MS, MSS .MPS .K 1 .K2.K3.K4 
KK  *0 
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Computer  Program  Listing  (continued) 


10  KK=KK+ 1 
D5=DL+  200. 

D6-0L+  400. 

11  T  5  =  TE  +  05 / A 
T6=TE+  06/A 

H5-MIN1FH (1./H1E  +  1 ./H2E >/ < T5#F*ABSF ! .007-.058«T5 ) ))  .  (15.)) 
H6=M1N1F((( 1./H1E+1./H2E ) /( T6*F*ABSF( .0o7-.058«T6 ) ) ) * ( 15. ) ) 

S5-H5  +  10.*ALOG10(F*T5**4)-.1*(NS-301. ) *EXPF ( -T 5 *05 /40 .  ) 
S6=H6+10.*ALOG10(F*T6»*4)-.1»(NS-3C1.) »E X PF ( -T6*D6 /40 . ) 

I F ( T5*D5  .LE.10.)  AH5 * S 5+ 1 03 . 4  +  . 33 2* T 5 *05 - 1 0 . *ALOG 1 0 ( T 5 *D5 ) 

I F ( T 6*06  .LE.10.)AH6»S6  +  1  03.4  +  .3  32  »T6*D6-10.*ALOG10( T6*D6 ) 

I F ( T5 #05  .GT.  10.  .AND.T5»D5.LE.  70.)  AH5 =  S5  +  97 . 1  +  .2 1 2 * T 5*D 5-2 . 5* 
CALCG10I T5*D5 ) 

I F ( T  6  #06  .GT.  10.  • ANO . T  6*06  .LE.  70.*  AH6  =  S6 +97 . 1  +  . 2 1 2  * T6 »D6-2 . 5 
C*ALOG 1 0 ( T6  *06 ) 

l F ( T  5*05  .GT.  70.)  Ah5 *S5  +86 . 8  +  .  1 5 7*T 5 #05  +  5 . *AL0G1 0 ( T5 *D5 ) 

I F ( T  6  #06  .GT.  70.)  AH6  =  S6 +86 • 8  +  . 1 5 7*T 6 *06  +  5 . * A  LOG  1 0 ( T 6*06 ) 

MS* ! AH6-AH5 ) / ' D6-D5  ) 

IF  UK  .  EQ.  2)  GO  TO  25 
IF  ( H 5  .LE.  10.)  GO  TO  30 
IF  UK  .EQ.  1  )  GO  TO  20 
25  MS=MSS 

AH5C= AH5 
AH5  * AH5  S 
05*05  S 
GO  TO  30 
20  DH  *  0 . 

DLl«DLSl*EXr'F(-.07»$QRTF(DH/HlE)  ) 

DL  2  =DL  S  2  *E  XPF ( -,07*SQRTF(0H/H2£ ) ) 

DL*0L1+DL2 

TE1*(  .00065/DLSl)*t  (0LS1/DL1-1. ) #DH-3.077*H1G) 

TE2  =  ( .00065/DLS2 )* ( (DLS2/DL2-1. ) #0H- 3 . 077 #H2G ) 
TE*MAX1F((TE1+TE2).( -DL/A  )  ) 

T-TT  ■  0/A 

AH5S*AH5 

mSS=mS 

0  5  S  =  . '  5 
GO  TO  10 
30  CON  T !  NUE 

RETURN 
END 


SUBROUTINE  LOS 

SUBROUTINE  TO  COMPUTE  LINE  OF  SIGHT  ATTENUATION 

COMMON  /M/F ,D .NS.A .OH.DHS .S.E ,  °OL . KM 
COMMON  /NR/  JZ »W ,SW3 .SW4, SA3 .SA4.S AFO 

COMMON  /MP/  HIE  »H2E .H1u.H2G.0LS 1 .DLS2.DL1 .DL2.DL .DLS.T  E 1 .TE2 .TE.KL 
COMMON  /MLOS/  AG.AD.AS.ACR.AED.MO.AH50.AH5.05.MS.AES.DX.H5 
COMMON  /ML/  DO.Dl.DOl .D02 . AO . A  1 .K1 ,K2 ,AL .ALS.AOG 
REAL  NS  »MD .MDC .MS.MSS.M0S.K1.K2.K3.K4.M 
KM*; 

CAH,  Oi;- 
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Computer  Program  Lifting  (continued) 


ICM*0 

CALCULATION  OF  TWO  RAY  THEORY 

DO.  «00004*H1E*H2E*F 
D02«MIN1F( (-AED/MD) * (DL-2. )  ) 

IF  (AED  .GE.  0.)  DO-MINlF!D0l*(.5»DL) ) 

IF  (AED  .LT.  0.)  D0».5*DL 

IF  (AED  .LT.  0.  .AND.  D02  .GE.  DO)  D0=D02 

D1*D0>.25* ( DL-DO ) 

IF  ( D 1  .LE.  DO)  D1=DO+.25*(DLS-DO> 

J*0  S  DS*D 
IF  (J.EO.  0)  D«D0 

D  I V* 1 .  *  PSI«ATANF< <HlE+H2E)/( 1000. *D) ) 

DHD-DH*( 1.-.8  *EXPF ( 02*D) ) 

SH-.78*DHD*EXPF(-.5*DHD**.25) 

L  *S I NF ( PS  I ) 

X«l8000.*5/F 

P2»  (SORT  F(  (E-COSF!  PS  I  )  *COSF  (  PS  I  )  )  *«  2  +  X  *X  ) -t-E-COSF  (  PS  l  )*COSF  (PS  I  )  ) /2 

C. 

P«SQRTF(P2) 

0«X/( 2 • *P ) 

IF  (POL  .EQ.  1. )B= (E*E+X«X) /(P2+Q*Q) 

IFIPOL  .EQ.  -l.)B-l./(P2+0«Q) 

IF  (POL  .EQ.  1. )M=2.*(P*E+Q*X)/(P2+Q*0) 

IF  (POL  .EQ.  -1. )M*2.*P/(P2+Q*Q) 

R2  *=  (  1  .  +  B*SP*SP-M*SP  )  /(  l.  +  B*$P*SP+M*SP  ) 

RE»$QRTF($P) 

SQEXF  =  SQRTF ( R2) »EXPF (-. 02 09564473* F »$H*SP ) «D I  V 
IF  (SQEXF  .GT.  .5  .AND.  SQEXF  . GT .  RE)  RE-SQExF 
C«ATANF(Q/(P*SP)  ) -AT ANF IQ/(P-SP  )  ) 

IF  (POL  .EQ.  -1 . )GO  TO  40 
Y 1 » ( X*SD+Q I / ( E*SP+P 1 
Y2*(X»SP-Q)/(E»SP-P) 

IF  (F*SP  .GE.  P)  C«ATANF(Y1)-ATANF ( Y 2 > ♦ 3 . 14 1 S926S4 

IF  ( E*SP  .LT.  P  .AND.  P*SP  .GT.  .5)  C* AT ANF ( Y 1 )♦ AT ANF < Y 2 > 

IF  (E*SP  .LT.  P  .AND.  P*SP  .LE.  .S)  C*  ATANF  »  i  1  ) -AT  ANF  (  Y  2  ) 

IF  ( J  .EQ.  0 ) 

CA0*-l0.»ALOG10( 1 .  +  RE*RE-2.»Rc«COSF ( .00  004 191 7*F *H l E »H2E /DO-C )  i 
IF  ( J  .EQ.  II  GO  TO  3 
D  =>  D 1  *  J*1 

GO  TO  22 
CONTINUE 
D-DS 

IF  ( J  .EQ.  1  ) 

CA1«-10.»ALO&10( l.+RE»RE-2.»RE»COSF ( .00004 19 17*F*H1E*H2E/D1-C ) ) 

COMBINATION  OF  TWO  RAY  THEORY  AND  DIFFRACTION 

AL  S* AED+MD*DL  S 
AL«AED+MD»OL 
DEDO- AED+MD»D0 
DED1«AED+MD*D1 
SAO-AO  *  S A  1  *  A  1 
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Computer  Program  Listing  (continued) 


W*(MINlF<HlE*H2E)/MAXlF<HlE#H2E)1/< 1  .  +  F *DH* . 000 1  > 

W  =  1 . / ( 1 .+F*DH*.000 1 ) 

A0=MIN1F< (W*A0+( l.-W)*DED0) ,0£D0) 

A 1  *M IN  IF ( I W*Al  +  ( l.-W)*DEDl  >  .DE01 ! 

10  <2  = ( ( ALS-AO! *1 01-00 )-( A1-A0 )*( DLS-DO) I / ( (Dl-D0)*ALOG10(DLS/D0)~ 

C ( DLS-DO ) *AL0G1 0 ( 01 /DO ) ) 

K2-MAX1C ( <2.0. ) 

Kl=< ( ALS-AO )-<2*ALOG 10 (OLS/DO) ) /(DLS-DO) 

IF  (<1  .GE.  0.)  GO  TO  50 

<1=0. 

K2*(ALS-A0)/( ALOG 1 0 ( DL S/DO )  ) 

5  0  AG  =  A0+<l*(D-D0)+<2*ALOG10( D/00  1 

IF  (AG  .IT.  0.)  AG=0. 

51  AOG«AO 

53  ACR=AG 

RETURN 
END 


3-40 


p 


Computer  Program  Output 


Colorado  plains  ns-290. 


rm  100*0  DM*  96.00  MiG"  A* 00 
AC«  28.09  Kl*  0;A9356  K2* 

ACO*  39.24  MO"  0.28151  AES" 
0  5.00  10.00  26.00 

A  35,24  39.72  *4.87 


H2G*  3.00  Tf"  0*004861  Dx»  133.88 
6.69916  DLS"  15.23  ALS"  43.53 

'69.68  MS*  0.05418_  AOX"  76.93 

30.00  50*06  80.00 

47.69  5302  61.76 


F"  100.0  OH"  9o.OO  MlO"  4.00 
AC"  24,93  Kl"  6»*0l59  K2* 

AED"  36*15  MO"  5.27636  AES" 
D  5,00  10*00  20.00 

A  31.91  36.05  41.68 


h28>  6.00  IE"  0.002464  OX"  130.11 

7.10267  DLS"  18.16  ALS"  41,17 

‘64,85  MS*  0*OS5?5^  ADX"  72.11 

30.00  50*06  80.66 

44.44  4979?  58*26 


F"  100.0  OH"  90.00,  MiG"  4.00 
AC"  22,62  Kl"  0i36l62  K2" 

AED"  33.81  MD"  0.27609  ~ES« 

D  5.0o  10.00  26.00 

A  29.47  33.44  39.23 


H2G"  9.00  TE"  0.001556  OX"  126.02 
7.20867  DLS"  20.41  alS"  39,44 

‘61.49  MS"  0*05641  ADX"  68.60 

30.00  50r06  80.00 

42,09  47.61  55.09 


F"  100.0  OH"  9o.OO  MJO"  4,00 
AC"  30,15  Kl"  6.32410  K2* 

ACD"  40*18  MO"  0.28417  AES" 
D  5.6o  10.00  26*00 

A  37.34  41.36  45.87 


H23"  3.00  Tf"  0.004861  OX"  135.25 

7.97226  DLS"  15.23  AlS-  44,51 

'71.29  MS"  0.054l8_  AOX"  78.62 

30.00  50*Oo  80.00 

48.71  54739  62.92 


F"  100*0  OH"  9o,00  Ml®*  4*00 
AC-  26.22  Kl-  0i24421  K2" 

AED"  36*16  MO"  0.27885  AES" 
0  5.00  10.00  26.00 

A  33,30  37.05  41,73 


m2G"  6.00  TC»  0.002*64  Dx»  129.21 
8.39353  DLS"  18.16  ALS"  41.22 

"64.98  MS*  0.05575  ADX-  72.19 

30.00  50>On  80.00 

44.52  50*1 0  58?46 


F«  100* 0  OH"  96.00  Ml®*  4.00 
AC"  23. TO  Kl-  6*23397  K2" 

ACO"  33*88  MO"  6.27858  AES" 
D  5.00  10.00  26.00 

A  30.79  34.51  39.40 


H2G-  9.00  TE"  0.001556  OX-  125.44 
8.47052  DLS"  20.41  AlS"  39.57 

61.75  MS"  0.05641  AOX"  68.83 

30.00  50»Oo  80. go 

42.24  47.8]  56.17 


F"  50.0  OH«  96.00  Ml°*  4.00 
Ac-  31.31  Kl"  1*01536  KZ* 
AED"  47.89  MD"  0.2967ft  ACS" 
0  5.00  10. 06  20.00 

A  42.87  49.60  53.82 


M2G"  6.55  Tf"  0.029474  DX"  121.70 
8.1428ft  DLS"  11.19  AlS"  51.21 

‘71.36  MS"  0.04634  ADX"  04.00 

30.00  50.06  80.00 

56.79  62.7?  71*63 


F"  So . 0  OH"  96.00  Hl®*  *.00 
AC"  29.15  Kl"  6*47  486  KZ" 
ACO"  42.92  MD"  5.2492o  AES" 
D  5,00  10.00  26.00 

A  38.16  43.39  47.90 


h2G"  1.7o  Tf«  8.608505  Ox-  136.42 
9.49046  DLS"  13.48  AlS"  46.28 

‘69.78  MS"  6.05227  ADX"  7ft. 91 

30.00  5o»0n  80.66 

50.40  55.3*  62.86 


F- 


AC" 

ACO" 

0 

A 


20.0  OH"  96.00  Ml®*  3.30 
30.25  Kl"  0i3U08  K2- 

43.48  MO«  0.2007*  AES- 

5.60  10.00  2o » 00 

39,49  44.36  47.50 


H2C"  1.36  Tf" 
10.97428  OLS" 
‘64.35  MS* 
30.00  50>0n 

49.51  53*5? 


0.011970  OX"  139.21 
12.07  ALf  45.91 
0.05084  AOX"  71.43 
80.66 
59.54 
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Computer  Program  Output  (continued) 


COLORADO  MOUNTAINS  NS-290* 


fm  100.0  DH-650.00  HlG-  *.00  H2G-  3.00  TE-  0.180463  DX»  131.08 

AE-  52, *2  Kl-  0«68S66  K2-  0*66909  DLS-  15.23  AlS-  63,70 

AEO-  59.29  MO-  0.28955  AES-  '91.24  MS-  0.0*574  aDX-  97.24 

0  5.00  10*00  20*00  30.00  50*0n 

A  56.33  59.97  65.08  67,97  73;76 

F-  100.0  OH-650. 00  HlG-  4.00  H2G-  6.00  TE-  0.117T12  Dx-  16*. 84 

AE-  47.62  Ki-  6*57113  K2-  0*69038  DLS-  18.16  AlS-  58.86 

AEO-  53*99  M0-  0.26795  AES-  "94.89  MS-  0*01984  AOX-  98.16 

D  5.00  10*00  20*00  30.00  50*06 

A  50.95  54.02  59.35  62.03  67.39 

E-  100.0  DH-650.00  Hl°"  4.00  H2G-  9 . 00  Tf-  0*100178  Ox-  164.32 

AE-  45.85  K 1  a  0*52997  K2-  0*63153  0L5-  20**1  *LS-  57.49 

AEO-  52*17  MD-  5.26061  AES-  '87.84  MS»  0*0*351  AOX-  94.99 

D  5.00  10*00  20.00  30.00  5o*On 

A  48.94  51.78  57.27  5«,99  65*26 

fm  100.0  OH-650. 00  H19-  *.00  H2G-  3.00  Tf»  0.180463  Ox-  125.95 

AE-  53,20  K 1  a  6*67804  K2-  0*80526  OLS-  15.23  ALS-  64.48 

AEO-  60*01  MO-  5.29374  AES-  91.24  MS-  _  0.04574  ADX-  97.00 

D  5.00  10*00  20.00  30.00  50;06 

A  57.16  60.79  65.88  68.82  74*76 

F-  100.0  OH-656. 00  M18-  4.00  H2G-  6.00  TE-  0.117712  DX-  161.69 

AE-  48.13  K 1  a  0*55941  K2-  0*83653  OLS-  18.16  AlS-  59.34 

AEO-  54.43  MO-  5.27074  AES-  "95,00  MS-  0.0198*  AOX-  98.20 

D  5.00  10*00  20.00  30.00  SOjOp 

A  51.51  54.56  59.84  62.55  67.96 

F-  100*0  OH-65Q.OO  h1«-  **00  H2Q-  9.00  TE-  0.100l"78  Dx-  161.52 

AE-  46,33  Kl-  6*51989  K2-  0*78175  DL5-  20**1  AL5-  57.96 

AEO-  52.59  MO-  0.26327  AES-  88.08  MS-  0*0*351  AOX-  95.11 

0  5.00  10*00  20*00  30.00  50;0p 

A  49,47  52.31  57.74  60.49  65*7? 

F-  50.0  DH-650.00  *10-  **00  H2G-  6*55  Tf-  1.482328  Ox-  163.97 

AE-  60.73  Kl-  1.87J81  K2-  1.59156  DLS-  11.19  AlS-  91.37 

AEO-  86.21  MO-  5. *6115  AES-  120.32  MS-  0*25309  aOX-  161.82 

0  5,00  1 6 • 0  6  26*00  30.00  5  0 • 0  n 

A  79.21  89.06  95. *3  lOO.O*  109*27 

F-  50*0  Om-650.00  Mi<J-  **00  H2G-  1.76  Tr-  0*305643  OX-  128.35 

AE-  51,08  Ki-  1.01365  K2-  2.02007  OLS-  13.48  AlS-  67.03 

AEO-  62*71  MO-  6.32004  AES-  95.06  MS-  __  0.06799  AOX-  103.79 

0  5.60  10.00  20.00  30.00  50;0ft 

A  57.56  63.23  69.11  72*31  78. 7i 

F-  20.0  OH-650. 00  HlG-  3.30  H2G-  1,36  Tr-  0. *37338  Dx-  146.37 

AE-  44,25  Kl-  1.523*6  K2-  4.18953  OLS-  12.07  AlS-  67.17 

AEO-  63*15  MO-  6.33315  AES-  '98.92  MS-  _  0-08873  AOX-  111,91 

0  5.00  10.00  26.00  30.00  50*0n 

A  54.80  63.68  69.81  73*14  79. 8o 
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Computer  Program  Output  (continued) 


OHIO  NS-312. 


F*  100.0  DH*  96.00  H1Q“  4,00 
AE*  28.13  Kl*  6:47405  K2» 

AED"  39.26  MD*  6.27418  AES* 
0  10.00  20.00  30.00 

A  39.60  44.74  47.48 


H2(3»  3.00  TE-  0.004762  Dx*  140.04 

6.72692  OLS*  15.55  AlS*  43.52 

‘69.81  MS"  0.05598  ADX*  77.65 

50.00 
52.97 


f*  loo.o  dh-  9o.oo  hio*  4.00 

AE*  24.98  Kl*  0*38606  K2* 

AED"  36*19  M0«  6.26912  AES* 

0  10. 00  20*00  30*00 

A  35.97  41.57  44.26 


H23«  6.00  TE"  0.002414  Dx*  135.88 

7.13023  DLS*  18.53  ALS*  41.18 

"64.94  MS"  6.65754  ADX*  72,76 

50.00 
49.65 


F*  100.0  OH*  96.00  Ml<3"  4.00 

AE*  22.67  Kl*  0*3*787  K2* 

AED"  33.86  MO*  6.26883  AES* 
D  1 0  •  0  0  20*00  36*00 

A  33.39  39.04  *1.92 


H20*  9.00  TE"  0.00152*  Ox*  131.49 

7.23686  OLS*  20*83  ALS"  39.46 

'61.55  MS"  O.OSBil  ADX*  69.21 

50.00 
*7.30 


F«  100.0  OH*  96.00  HlO"  4.00 
AE*  36.18  Kl*  6*30841  K2* 

AED"  40*19  MO*  6.27675  AES* 
D  10.00  20*00  30*00 

A  41.26  45.73  48.50 


H2G*  3.00  Tf"  0*00*762  Dx*  141.46 
7.99234  OLS*  15.55  ALS*  44,50 

'71.43  MS"  0*05598  ADX*  79.34 

50.00 
54.03 


F*  100.0  OH*  96,00  Hl9"  4.00 
AE*  26,26  Kl*  0:23228  K2* 

AEO"  36.20  MO*  6.27151  AES* 
0  10.00  20*00  36*00 

A  36.99  41.63  44.34 


H2G*  6.00  TE"  0.00241*  Ox*  134.94 
8*41271  OLS*  18.53  ALS-  41,23 

65.07  MS"  0*05754  ADX*  77.83 

50.00 
49.77 


F*  100.0  DH*  96.00  *10"  4,00 

AE"  23.74  Kl*  0:22313  K2* 

AEO"  33.94  MD*  0,27121  AES* 
D  10.00  20 * 00 '  30*00 

A  34.46  39.2$  42.07 


h2Q*  9.00  TE"  0.00152*  DX*  130.88 
8*49031  DlS*  20*83  ALS*  39.59 

"61.81  MS"  0*05821  ADX*  69.43 

50.00 
47,50 


F*  56.6  dh*  96.00  H10"  4,24 

AE*  36.00  Kl*  0’61595  K2* 

AED"  44.45  MO*  6.25734  AES* 
D  10.00  20*00  30.00 

A  45.10  49.60  52.17 


H20*  1.00  TE"  0.014366  Dx*  138.42 

8.93826  DLS*  12.74  ALS*  47.73 

"72.89  MS"  0.05189  ADX*  80.07 

50. 00 
57.32 


F*  50.6  DH*  96.00  M1G"  4.24 
AE*  27.49  Kl*  6:39565  K2* 

AED*  41.53  HO*  0.23414  AES* 
0  10.00  20*00  30.00 

A  41.04  46.22  48.56 


H2G*  3,00  TE*  0.004564  Dx*  144.92 
9.59105  DLS*  15.79  ALS*  45.23 

'67.33  MS"  0.05610  ADX*  75.46 

50.00 
53.24 


F* 


AE" 

AEO" 

0 

A 


20.0  DH*  96.00  Mla"  3.68 
27,77  Kl*  0*l’l76  K2* 

41.83  MD*  6.18600  AES* 
10.00  20*00  30.00 

41,50  45.55  47.41 


M20*  3.00  Te»  0.005664  Dx*  148.33 

11.81777  DLS*  15.21  AlS*  44.65 

"61.14  MS"  0.05580  ADX*  69.41 

50.00 
51.13 
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